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Zusammenfassung 

In dieser Arbeit werden Koharenzeffekte in der tiefinelastischen Streuung (DIS) und 
im Drell-Yan (DY) Prozess an Kernen untersucht, insbesondere der Shadowing Ef- 
fekt. Es wird im Ruhesystem des Targets und in der Farbdipol Formulierung gear- 
beitet. Die Glauber-Gribov Theorie fiir Mehrfachstreuung im Kern wird so modi- 
fiziert, dass der Formfaktor des Kernes in alien Streutermen beriicksichtigt ist. Ferner 
wird die mittlere Koharenzlange fiir einen Fockzustand definiert. Damit ist es moglich 
abzuschatzen, dass das Gluon- Shadowing fiir xbj > 0.01 vernaclilassigbar ist. Pa- 
rameter freie Rechnungen werden mit Daten von NMC und E665 fiir DIS und mit 
E772 Daten fiir DY verglichen. In beiden Fallen wird gute Ubereinstimmung fest- 
gestellt. Der von HERMES beobachtete Effekt kann jedoch nicht reproduziert wer- 
den. Fiir DY-Dileptonen aus Proton-Kern Kollisionen bei RHIC Energien wird fiir den 
gesamten xp Bereich deutliches Shadowing vorrausgesagt. Der Einfluss des Kerns auf 
die Transversalimpuls-Verteilung der DY Paare wird ebenfalls untersucht. Des weit- 
eren wird eine neue Parametrisierung des Dipol-Wirkungsquerschnittes prasentiert. 



Abstract 

In this work, coherence effects in deep inelastic scattering (DIS) and in the Drell-Yan 
(DY) process off nuclei are investigated, in particular nuclear shadowing. The target 
rest frame and the color dipole formulation are employed. Multiple scatterings are 
treated in Glauber-Gribov theory, which is modified to include the nuclear form factor 
to all orders. Based on the mean coherence length, which is defined in this work, it is 
estimated that gluon shadowing is negligible at xbj > 0.01. Parameter free calculations 
are compared to NMC and E665 data for DIS and to E772 data for DY. In both cases, 
good agreement is found. It is however not possible to reproduce the effect observed 
by HERMES. For dileptons in proton- nucleus collisions at RHIC energies, considerable 
shadowing for the whole xp range is predicted. The influence of the nucleus on the DY 
transverse momentum distribution is also studied. Furthermore, a new parametrization 
of the dipole cross section is presented. 
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...all exact science is dominated by the idea of approximation. 

Bertrand Russel 
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1 INTRODUCTION 2 

1 Introduction 

The use of nuclei instead of protons in high energy scattering experiments, hke deep 
inelastic scattering, provides unique possibilities to study the space-time development 
of strongly interacting systems. In experiments with proton targets the products of 
the scattering process can only be observed in a detector which is separated from the 
reaction point by a macroscopic distance. In contrast to this, the nuclear medium 
can serve as a detector located directly at the place where the microscopic interaction 
happens. As a consequence, with nuclei one can study coherence effects in QCD which 
are not accessible in DIS off protons nor in proton-proton scattering. An important 
question that can only be answered with the help of nuclei is for instance, how quarks 
and gluons evolve from the early stages of a collision to the hadrons which are finally 
observed in the detector. The large extension of the nuclear medium makes it possible 
to investigate, by which time scales this hadronization process is governed. Note that 
the radius of a heavy nucleus like lead is approximately eight times as large (~ 6.8 fm 
in the nuclear rest frame) as the radius of a proton (?» 0.86 fm). 

This work is mostly concerned with the theoretical analysis of coherence effects in 
DIS off nuclei and in Drell-Yan (DY) dilepton production in proton-nucleus scattering, 
in particular with the phenomenon of nuclear shadowing. Before we turn to nuclear 
targets, we shortly review the physics of a proton target. 

In DIS, a lepton is scattered off the target. This lepton radiates a virtual photon, 
which can resolve the microscopic substructure of the proton. Therefore, the huge 
colliders like HERA, where such experiments are performed, are big microscopes that 
allow to investigate, how the proton is made up of quarks. It is well known today that 
the proton contains three quarks which are called valence quarks. These valence quarks 
carry the quantum numbers of the proton and are the analog of the valence electrons 
which are responsible for the chemical properties of an atom. In addition to the quarks, 
there are also gluons inside the proton. Gluons mediate the color forces between the 
quarks. Since gluons have no electromagnetic charge they are not directly observable in 
DIS. One can however conclude that there must be neutral partons in the proton since 
quarks carry only about half of the momentum of the proton. The missing momentum 
is believed to be carried by gluons. It is natural to ask, whether quarks and gluons 
are elementary or whether they have a substructure themselves. In order to find an 
answer to this question, one has to increase the resolution of the microscope, i.e. build 
colliders with higher energies which can measure at higher momentum transfer. Today, 
the highest energies are reached at HERA where positrons collide with protons at a 
center of mass (cm.) energy VS ~ 300 GeV. As the resolution is increased, one finds 
that a quark, carrying the longitudinal momentum fraction Bjorken-x^j of the proton 
consists of a quark and a gluon which carry both together the momentum of the parent 
quark and therefore each a smaller momentum fraction of the proton. In addition, also 
gluons can split into quark-antiquark [qq) pairs. This is confirmed experimentally by 
the increase of the quark density at low values of xbj < 0.1. At very low xbj -C 0.1, 
the partonic content of the proton is dominated by gluons and the photon sees only 
the gg-pairs which originate from gluon splitting. However, one cannot distinguish in 
DIS whether the virtual photon couples to a quark or to an antiquark. Complementary 
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Figure 1: At low xbj and in 
the target rest frame, the vir- 
tual photon converts into a qq- 
pair long before the target. The 
quark carries momentum frac- 
tion a of the 7*, the antiquark 
1 — a. The transverse separation 
between the particles is denoted 
by p. The curly line represents a 
gluon. 



information about the antiquark density is provided by the Drell-Yan (DY) process. In 
the DY process in proton-proton collisions, a quark from the projectile can annihilate 
with an antiquark from the target and produce a massive photon. This photon decays 
into a lepton pair which can be detected. 

How does a nucleus look like at high energies, i.e. at low x^j? The answer depends 
on the reference frame. In a frame where the nucleus is fast moving, the so-called 
infinite momentum frame, the nucleus is strongly Lorentz contracted. However, the 
localization of gluons which carry only a very small momentum fraction of the nucleus 
is determined by the uncertainty principle. In the infinite momentum frame, the cloud 
of these Iow-xbj gluons extends over the whole nucleus and the nucleons are able to 
communicate with each other. The situation looks different in the rest frame of the 
nucleus. In the nuclear rest frame, the nucleons are well separated from each other by 
a distance of ~ 2 fm. How can these two pictures, infinite momentum frame and rest 
frame, be reconciled with each other? Of course, all observables have to be Lorentz 
invariant. 

Note that not only the partonic structure of the nucleus is frame depended, but 
also the partonic interpretation of the scattering process. At high energies, nuclear 
scattering is governed by coherence effects which are most easily understood in the 
target rest frame. In the rest frame, DIS looks like pair creation from a virtual photon, 
see fig. |l[ Long before the target, the virtual photon splits up into a gg-pair. The 
lifetime Ic of the fluctuation can be estimated with help of the uncertainty relation 
to be of order ~ I/itinXbj (cf. section |27^ ) where niN ~ 1 GeV is the mass of a 
nucleon. The coherence length can become much greater than the nuclear radius at 
low Xbj- On a nuclear target, the pair will experience multiple scatterings off different 
nucleons within the coherence length. This corresponds to the overlap of gluon clouds 
from different nucleons in the infinite momentum frame. The long lifetime of the qq 
fluctuation, which extends over the whole nucleus, leads to the pronounced coherence 
effects observed in experiment. The measurable cross section is independent of the 
reference frame, but our physical picture depends on it. The target rest frame is 
especially well suited for the study of coherence effects. 

The most prominent example for a coherent interaction of more than one nucleon is 
the phenomenon of nuclear shadowing. Naively one would expect that the cross section 



1 INTRODUCTION 4 

for scattering a lepton off a nucleus with mass number A is A times as large as the 
cross section for scattering the lepton off a proton. In experiment it is however seen 
that the nuclear cross section is significantly smaller. Shadowing in low xbj DIS and 
at high photon virtualities was first observed by the EM collaboration [|T|. The same 
reduction of the cross section was measured by E772 for the Drell-Yan (DY) process 
at low X2 0. 

What is the mechanism behind this suppression? If the coherence length is very 
long, as indicated in fig. |I|, the gg-dipole undergoes multiple scatterings inside the 
nucleus. The physics of shadowing in DIS is most easily understood in a representation, 
in which the pair has a definite transverse size p. As a result of color transparency 
0, ^, H, small size pairs interact with a small cross section agg{p), while large pairs 
interact with a large cross section. The term " shadowing" can be taken literally in the 
target rest frame. The large pairs are absorbed by the nucleons at the surface which cast 
a shadow on the inner nucleons. The small pairs are not shadowed. They can propagate 
through the whole nucleus. From these simple arguments, one can already understand 
the two most important conditions for shadowing. First, the hadronic fluctuation of 
the virtual photon has to interact with a large cross section and second, the coherence 
length has to be long enough to allow for multiple scattering. At xbj -^ the coherence 
length becomes infinite and shadowing in DIS can be calculated from a formula very 
similar to the Glauber-Gribov eikonal formula f^, 01) if one knows the cross section 
(Tqq{p) for scattering a gg-pair with transverse size p off a nucleon. The application of 
the eikonal formula is possible because at infinitely high energies gg-dipoles with fixed 
separation in impact parameter space are eigenstates of the interaction 0, ||, H . The 
dipole cross section is the main nonperturbative ingredient to all calculations in this 
work. 

Nuclear shadowing has received much attention within the last two decades. A 
detailed review of different approaches to nuclear shadowing can be found in [10|, see 
also the more recent review [P^. 

Note however that shadowing has not yet been measured at RHIC {y/S ~ 200 
GeV) nor LHC (Vd ~ 5.5 TeV) energies. Furthermore, it is only known for quarks. 
Shadowing is also expected for the nuclear gluon distribution. This gluon shadowing 
could have a strong influence on the expected quark gluon plasma formation at RHIC 
and LHC and is only poorly understood at present. Only quark shadowing is calculated 
in this work, but the light-cone approach can also be extended to gluon shadowing |12| . 
Further insight into the physics of nuclear parton distributions could be provided by 
DIS experiments off nuclei at eRHIC or at HERA. Nuclei could also be included into the 
THERA option at TESLA (y/S ~ 1.4 TeV). The advantage of using nuclear targets is 
that parton densities which could be probed only at LHC energies with proton targets 
are already accessible at HERA energies. Such experiments are not only important in 
view of the quark gluon plasma. They can also clarify, up to which energies the QCD 
improved parton model can be applied. 

While the target rest frame picture of DIS is very popular, the light-cone approach 
of Kopeliovich |jT3|, |I4|, ^, which describes the DY process in the target rest frame, is 



less known. In the light-cone approach, DY dilepton production in the rest frame of 
the target appears as bremsstrahlung, see fig. 0. A quark from the projectile scatters 
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Figure 2: In the target rest 
frame, DY dilepton produc- 
tion looks like bremsstrahlung. 
A quark or an antiquark in- 
side the projectile hadron scat- 
ters off the target color field 
and radiates a massive pho- 
ton, which subsequently de- 
cays into the lepton pair. The 
photon can also be radiated 
before the quark hits the tar- 
get. 



off the target and radiates a virtual photon. This photon decays into a lepton pair. 
Remarkably, the DY cross section can be expressed in terms of the same dipole cross 
section that appears in DIS, cf. section |2.6| . This is a result of the QCD factorization 
theorem [|14|. On a nuclear target, the quark will of course scatter several times. The 
effect of multiple scattering on bremsstrahlung is well known in QED as the Landau- 
Pomeranchuk-Migdal (LPM) effect [||, p, |lll' see also ]19|]. The LPM effect leads to 



a reduction of the cross section due to destructive interferences in multiple scattering 
within an amorphous medium. An electron incident on a dense target with n scattering 
centers, will radiate bremsstrahlung after the first scattering. If the electron energy 
(in QED) is high enough, the longitudinal momentum transfer in a single interaction 
is small and according to the uncertainty relation, the electron needs a long time to 
recreate it's electromagnetic field. Within this coherence time, the electron can travel 
macroscopic distances. This is in complete analogy to the case of DIS. If the electron 
is hit several times within this length, it cannot radiate again, because it has not yet 
recreated it's field. Therefore, the overall cross section is less than n times the single 
scattering cross section. It is remarkable that a microscopic process, where all momenta 
are large, can exhibit a length scale of order millimeters (in QED). 

Note that both, shadowing in DIS and DY, can be calculated from the Glauber 
eikonal formula at infinitely high energies, where the transverse motion of the particles 
in the pair, fig. [^, or the incident quark, fig. Q, can be neglected. In this work, a Green 
function approach is developed that allows to do calculations at energies which can be 
reached in experiment. 

Note that shadowing in DIS can also be regarded as LPM effect for pair pro- 
duction. The physical interpretation is however much less intuitive. Both processes, 
bremsstrahlung and pair production, were studied in |TB[. After the experimental dis- 
covery of the LPM effect for bremsstrahlung in QED 
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J, more than forty 

years after its theoretical prediction, this effect again received much attention. In 
particular the QCD variant of the LPM effect was extensively studied with respect 
to suppression of gluon radiation and energy loss of particles propagating in a quark 
gluon plasma. The two most important approaches are the path integral approach 



by Zakharov |2^-[^ and the diagrammatic technique of the BDMPS collaboration 
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|SH]. It is demonstrated in 



work on the LPM effect in heavy ion colhsions is done in [37, 



present theoretical understanding of the LPM effect can be found in |^ 



that the two formulations are equivalent. Further 

l9| . A review of the 
The formula- 



tion of nuclear shadowing presented in this work is the one of Zakharov since our both 
approaches are extensions of Glauber-Gribov theory. The treatment of shadowing in 
DIS as LPM effect was first suggested in |^^ and elaborated in |^, ^ . 

This work mainly consists of two parts. In the first part, section |^, DIS and diffrac- 
tion on a proton target and DY dilepton production in proton-proton scattering are 
considered. The deep inelastic structure functions and the definitions of kinematical 
variables are introduced in section |2.1| . In section 2^, the parton model and the QCD 
improved parton model are shortly reviewed. The first two sections can be skipped by 
a reader familiar with these concepts. The space-time picture of low xbj DIS in the 
color dipole formulation is discussed in some detail in section p.3| . The basic concepts of 
the light-cone approach and the techniques for treating nonperturbative effects, which 
are employed throughout this work, are presented in this section, too. Section P]^ is 
a continuation of the preceding section and contains a discussion of diffraction in the 
color dipole picture. The sections ^ and TE treat the DY process. While in section 
275| the well known parton model of the DY process is shortly explained, section |2^ in- 



troduces the color dipole formulation of the DY process. This formulation is employed 
in section |^ where nuclear targets are considered. Finally, a new parametrization of 
the dipole cross section is presented in section |2.7| . All results of numerical calculations 
for proton targets and comparison with experimental data are shown in this section. 

The light-cone approach is extended to nuclear targets in section |^. This section 
contains the main results of this work. The connection between nuclear shadowing in 
DIS and diffraction is discussed in section B.ll. In addition, the Glauber-Gribov multi- 



ple scattering theory is shortly introduced and the physical conditions which have to 
be fulfilled in order to observe shadowing are explained. 



In the sections |3.2| and p.3| , it 
is explained how the nuclear form factor can be introduced in all multiple scattering 
While section |3.2| is intended to explain the physics of our approach, a more 



terms 

formal derivation and a discussion of approximations which are employed is given in 
section |3.3| . The mean coherence length is defined in section |3.4| . This quantity is in- 

In section ^31 an estimate is given, on 



tended to be a tool for qualitative considerations 
how many dilepton pairs in the low mass region are produced via the bremsstrahlungs- 
mechanism. Proton-nucleus and nucleus-nucleus collisions are considered in this sec- 
tion, but no shadowing is taken into account. Shadowing for DY dilepton production 



is discussed in section p.6| , where also predictions for RHIC are presented. Section |3] 
contains the comparison of our calculation with data for shadowing in DIS and DY. 
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2 Soft contributions to hard QCD reactions 

Deep inelastic scattering (DIS) and the Drell-Yan (DY) process are the two classical 
examples for the application of perturbative QCD (pQCD). The high virtuality of the 
photon or the large mass of the dilepton pair, respectively, provides the hard scale, 
which is necessary for perturbation theory. Due to confinement, the fundamental QCD 
degrees of freedom, quarks and gluons are not directly observable but occur only in 
bound states. Therefore all QCD reactions also involve a soft scale, given e. g. by the 
radius of the hadron. 

Such two-scale processes have been extensively studied for more than 20 years and 
their theoretical treatment has reached a high level of sophistication. Basic to the 



standard approach is the operator product expansion (OPE), invented by Wilson P3 |. 
Combined with asymptotic freedom, the OPE is a systematic way to separate hard 
and soft scales in DIS, DY and other hard processes. The observable cross sections can 
be written in factorized form, namely as a convolution of a hard partonic cross section 
and of the soft parton distribution of the hadrons (or of hadronic matrix elements, 
more generally speaking). Contributions which do not factorize vanish in the limit 
Q"^ ^ oo. While the partonic cross section is governed by the hard scale and can be 
calculated perturbatively, the parton distributions are of completely nonperturbative 
origin and no way is known to calculate them from QCD. Even lattice calculations are 
not applicable in the kinematical region of DIS and DY. 

Instead, the parton distributions have to be extracted from experimental data. 
Their most important property is, that they are independent of the particular hard 
process under consideration and depend only on the hadron. This universality is cru- 
cial to make the theory predictive. Today, several collaborations exist which provide 
parametrizations of these distributions [Q, ^ ^. The parton distributions are typ- 
ically given at a semihard input scale Qq. The evolution of the parton distributions 
to higher virtuality Q^ is then correctly described by the DGLAP evolution equations 
^-fS^]. This is one of the great successes of QCD and an important argument, that 



QCD is the correct field theory of the strong interaction. 

One should however bear in mind, that the soft physics is completely described in 
terms of fit functions and nothing is known about the nonperturbative mechanisms, 
leading to the specific shape of the parton distributions. In particular, all coherence 
effects, which are so important in high energy nuclear physics, are hidden in these 
parametrizations. The aim of this work is to develop an approach which makes the 
physics underlying these coherence effects transparent and intuitively understandable. 

The purpose of this sections is to introduce the basics of the light-cone approach 
to DIS and DY before nuclear targets are considered in section |^. The fundamental 
input to all calculations is the cross section for scattering a quark-antiquark dipole off a 
proton. A parametrization of the dipole cross section is presented, which describes well 
all F2 data from HI and ZEUS at Q^ < 35 GeV^ and also reproduces total hadronic 
cross sections. We do not aim to achieve the same level of sophistication as the OPE 
approach in this work. We rather intend to develop a good physical understanding 
of high energy DIS and DY in the framework of the light-cone approach, which is 
especially well suited to describe coherence effects in multiple scattering off nuclei. 
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This section also serves as a reminder of the definition of structure functions and 
kinematical variables. 

2.1 Description of DIS in terms of structure functions 

The natural language to describe hard QCD processes is based on structure functions 
which contain all information about the substructure of the target. These functions are 
defined by relating them to the experimentally observable lepton-target cross section 
and can be expressed in terms of hadronic matrix elements. They acquire a physical 
interpretation within a parton model of the structure of the target. The definition of 
structure functions in DIS is explained in more detail in standard textbooks, see e. g. 



5T| , |52| , |53| , we will only give a short overlook. Only the scattering of charged leptons 



is discussed in this work. Deep inelastic neutrino scattering, as well as Zq exchange are 
not considered. 




X 



Figure 3: Inclusive deep-inel- 
astic lepton-proton scattering. 
A charged lepton with four- 
momentum k^ = {E, k) scat- 
ters off a proton with four-mo- 
mentum Pf" = {Ep,P). In DIS 
a spacelike photon with mo- 
mentum g^ = —Q^ < is ex- 
changed. Only the scattered 
lepton with momentum k^' is 
observed in inclusive measure- 
ments. 



The kinematics for inclusive DIS is explained in fig. |^. A lepton scatters off a proton 
or a nucleus. Due to the impact of the virtual photon, the target breaks up into an 
unobserved final state X. Only the scattered lepton is observed in the final state. We 
work with the standard kinematical variables. The virtuality of the photon is denoted 
by 

Q^ = -q^ = {k-k')^>0 (1) 

and the mass of a proton by 

ml = P^^l GeVl (2) 

Furthermore the lepton-nucleon center of mass ( cm) energy squared is defined as 

S = {k + Pf (3) 

and the same for the 7*-nucleon system 

s = {q + PY. (4) 
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For convenience, we have chosen our convention shghtly different from the standard 
convention, where s is denoted by W^. The Bjorken variable is given by 

= -91- = -91- ~ -91- (^^ 

^""^ 2P-q 2mMi^^Q^ + s ^' 

and the relative energy loss of the lepton is 

P ■ q Q' + s 



We will also frequently use the Lorentz invariant variable 

P-q 

V = . 

friN 



(7) 



In one photon exchange approximation, the differential cross section for inclusive 
scattering reads 

da{eP -. e'X) = ,^j^^^^, E(2^)''^'(^ + k - Px - k')\A\\ (8) 

The sum over final states is understood to include also the integration over phase space. 
The absolute square of the matrix element reads 

\A' = -QPUP\J'Ho)\x){x\rmP) (9) 

after summation of the polarizations of the virtual photon and in one photon exchange 
approximation. Here, J'^ is the electromagnetic current. The leptonic tensor is known 
as 

l^, = 2ik'^K + klk^-g^,k'-k). (10) 

We have averaged over spins of the initial lepton and summed over the final spins. All 
information about the target resides now in the hadronic tensor defined by 

W>^^ = Y.{P\rmX){X\.rmP)i2iTYS\P + q-px) (11) 

X 

- /<''--(^l^^(^)i"(0)|P). (12) 

In order to obtain (0), one has to make use of the integral representation of the 6- 
function and of the completeness of the final hadronic states X. The hadronic tensor is 
then related to the discontinuity of the forward virtual Compton scattering amplitude 
T^"" via 

p^M- = ^DiscT'^'^ = hm -l-iT^^'iqo + is) - T^^qo - ie)), (13) 

ZTT £^0 47ri 



with 



T-'^i/.%e-(P|T(.-(.).^(0))|P)^ (14) 
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- Im 




Figure 4: The hadronic tensor W'^'^ is related to the imaginary part of the forward 
{P = P' , q = q'), virtual (Q^ > 0) Compton scattering amplitude via the optical 
theorem. 



Here, T is the time ordering operator. Note that the discontinuity of T equals the 
imaginary part of T since the current operators J^ are hermitian. As illustrated in fig. ^, 
the relation between the hadronic tensor and the Compton amplitude is a manifestation 
of the optical theorem. 

Making use of Lorentz covariance, gauge invariance and parity conservation, one 
finds that the most general structure of W^"^ is 



W^" 



-9 



flU 



q^q"" 



Wi{xBj,Qn+\P^ + 



qh 



2x 



Bj 



pv 



q" 



2x 



Bj 



W2{XB,.Q^). 



The dimensionless invariant structure functions for DIS are commonly defined as 



Fi{xBj.Q^) 
F2ixBj,Q^) 



WiixBj,Q^) 
l^W^ixBj^Q' 



(15) 

(16) 
(17) 



and with this notation, the cross section for ep — > e'X reads 



dV 



Airai 



dxBjdQ^ Q^ 



y 



XBjy^N\ F2{XBJ,Q' 



Q' 



Xbj 



y^F^ixBj.Q' 



Since the only purpose of the lepton is to radiate the virtual photon, it is convenient 
to think about DIS as 7*p-scattering and to define corresponding cross sections for 
transverse and longitudinal photons. The cross section for a virtual photon with helicity 
A can be defined as 



o-A 



K 



'-e,{\)el{X)\mT^r 



(19) 



Note that the flux factor of a virtual photon is not well defined. We employ the 
convention of Hand Hi, i. e. 



K = 2P-q-Q'^ 

With this convention one obtains 

47r^ap 



2s. 



Cx 



(^L 



(5^(1 - XBj 

Hin Otem 
(52(1 -XBj 



-2xBjFi{x,Q' 



1 + ^ j F2{xbj,Q^) - 2xBjF^{xBj, Q' 



(20) 



(21) 



f22) 
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Here, ar is the cross section for transverse photons and ai for longitudinal photons. 
One also often introduces the helicity structure functions, 

FrixB^Q') = ^^^f^cr^, (23) 

'±71 Clem 

F,(.B„Q') = ^^-L. (24) 

The invariant structure functions Fi and F2 can be expressed in terms of the helicity 
structure functions Ft and Fl and vice versa. In particular, one obtains in the Bjorken 
limit, Q^ — > 00, xbj fixed, 

F2{xBj,Q'') = Ft{xbj,Q^) + Fl{xbj,Q^) (25) 

F,(x.„Q^) = ^^^^g^. (26) 

There are only two independent structure functions in DIS off unpolarized targets. 

2.2 The parton model of DIS 

It came as a big surprise, when the first DIS measurements at SLAG [^ showed that 
the structure function F2{xbj,Q'^) is nearly constant as function of Q^ at fixed xbj- 
An explanation of this phenomenon was given by Bjorken [Q and by Feynman |]57| . 



Feynman's intuitive explanation is depicted in fig. |^a. He proposed that the proton is 
made up of pointlike charged constituents, so called partons, and the total 7*p cross 
section is the incoherent sum of photon-parton cross sections. The transverse momenta 
of the partons are neglected, which is well justified in a frame where the proton is fast 



moving, e. g. in the Breit frame or in the infinite momentum frame [^. In such a 
frame the longitudinal momenta of the partons are much larger than the transverse 
ones. 

The cross sections for transverse and longitudinal photons scattering off spin-1/2 
partons, i. e. quarks (fig. ^), are given by 

An'^a^m.Zj f XBj 



aV = \ 5(l-^ (27) 

Q^{1 - XBj) V ^q J 



^l' = 0, (28) 

where Xq is the momentum fraction of the proton carried by the struck quark and 
Zf is the flavor charge in units of the elementary charge. The (5-function arises from 
momentum conservation and gives a physical meaning to the Bjorken variable. In the 
Breit frame, xbj is the momentum fraction of the proton carried by the struck quark. 
For massless quarks, the longitudinal cross section is zero due to helicity conservation. 
Introducing the density Qjixsj) of quarks of flavor /, inside the proton, one obtains a 
simple partonic interpretation of the structure functions, 

F2{xBj) = XBj^ZJ{qf{xBj) + qf{xBj)), (29) 

/ 
Fl = 0. (30) 
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^Bj 




X 



(a) 




(b) 

Figure 5: The parton model for DIS (a) and its QCD improvement (b). In the 
naive parton model, the virtual photon hits a quark (q) inside the target. In the QCD 
improved parton model, the quark can e. g. be created from a gluon, denoted by a 
curly line, by the splitting process G -^ qq. The Xi are the longitudinal momentum 
fractions of the particles with respect to the target. Therefore xbj < Xi < X2 < ■ ■ ■■ 



In the naive parton model, the structure functions depend only on xbj and not on Q^, 
the longitudinal structure function vanishes and one obtains the Callan-Gross relation 



F2 - 2xBjFi = 



(31) 



These equations are the basic results of the parton model and they are approximately 
confirmed by experiment. 

One of the great achievements of QCD is the successful description of the deviations 
from the naive parton model seen in experiment. In particular at low xbj deviations 
from Bjorken scaling become quite pronounced, see fig. ^. In the QCD improved 
parton model, perturbation theory is applied to calculate corrections to the parton 
model predictions. An example of such a correction, which is important at low xbj, 
is shown in fig. ^d. The quark seen by the photon is generated perturbatively from a 
gluon by the splitting process G — *> qq. Higher order corrections allow also to take into 
account further gluons radiated from the gluon that splits into the gg-pair. 

The evolution of the structure function with Q^ is described by the DGLAP equa- 
tions. There are several ways to motivate these equations. The DGLAP evolution 
equations resum logarithms in Q^ originating from ladder diagrams like the one in 
fig. ^D. Historically, they originate from a time even before QCD. Gribov and Li- 
patov E^, Esi were the first who found that these ladder diagrams receive dominant 
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Figure 6: Deviations from Bjor- 
ken scaling can be explained by 
pQCD corrections to the naive 
parton model. While Bjorken 
scaling is approximately fulfilled 
at high XBj, scaling violations 
become stronger and stronger as 
XBj decreases. At very low xbj, 
practically no scaling can be ob- 
served any more. The figure is 
taken from 

from ZEUS M, CTl, E665 
and NMC 



I[ and the data are 

m 0. 



contributions from configurations with strong ordering in the transverse momenta, i.e. 



pIi > Pt2 > 



(32) 



They obtained integrals of the form dp^/p\ which could be resummed to all orders. 
However, their calculations were done in an abelian theory. A very intuitive derivation 
in non-covariant Hamiltonian perturbation theory was given later (in QCD) by Altarelli 
and Parisi ^^ . The illustrative interpretation of the DGLAP equations is due to them. 
The photon acts as a microscope with a resolution determined by Q^. As one increases 
Q^, the photon resolves more and more of the target substructure. The probability 
to find a quark and an antiquark of fiavor / inside a gluon, fig. ^3, is then described 
by a splitting function Pfc- Similarly, there exists a splitting function which describes 
the probability to find a quark and another gluon inside the parent gluon. Finally, 
Dokshitzer [^ derived the same equation again with the method proposed by Gribov 
and Lipatov, but now in QCD. 

From the process where the quark seen by the photon is generated by the splitting 
of a gluon, one obtains in perturbation theory the contribution 

F2{XBJ,Q^) 



XBj 



+ 





Qf (xBj) 






1 
as f 
2tt J 

XBj 


|^G(.,{P. 


fxBj^ 


ln(«>...}l 



(33) 



to the structure function. Here, G{xi) ist the gluon density of the target and xbj/xi < 
1 is the momentum fraction of the gluon carried away by the quark in the ladder. 
The logarithmic divergence in (|33D is regularized by the cutoff k. The dots denote 
renormalization scheme dependent functions which can be calculated and contain no 
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such divergence. These functions are of no interest to our more quahtative discussion 



and we refer to |52] for more details. The sphtting function is given by 



Pfciz) = ^ [.^ + (1 - zf] . (34) 

The divergence in ( |55D is a coUinear divergence, originating from configurations with 
Pt — ^ in fig. Pd. This hmit corresponds to a long range part of the strong interaction 
which cannot be calculated in perturbation theory. Therefore, the quark (g/) and 
gluon (G) densities in (^) are unmeasurable, bare parton distributions, g/o and Gq. 
The renormalized quark density may be defined as 



qf{xs,,l^F) = qfoixB,) + g £ f^^o(xi) {P/G (^) In (§)+...} . (35) 

The coUinear singularity is absorbed into the bare quark density at the factorization 
scale fip- With this prescription, the partonic interpretation of F2 (pOf ) is maintained, 
but the structure function is no longer independent of Q^. 

The distribution qf{xBj, I^f) cannot be calculated in pQCD, since it is dominated 
by soft physics. It has to be extracted from measurements of F2. However one can 
calculate in pQCD, how the parton distribution runs with the factorization scale. Since 
^p is not a physical quantity, observables cannot depend on it. Therefore, 

dF2{xBj,Q^,f^F) ^ Q ,2g^ 

dfip 



g,JjAE^ = ^ r ^P.c f^ ) G(.„Q% (37) 



and 

dQ'^ 271 J,^^^ xi \xi 

This is one of the DGLAP equations which describes the evolution of the quark density. 
Note that also the gluon density depends on Q^, because of pQCD corrections. 
Altogether one obtains A'^^ + 1 coupled equations describing the Q^ evolution of the 
singlet parton densities, 

d (^A-B,.Q^)\ a. }dxj ""ffi"^) Pf-i"^)] 



Q 



dQ' \ G{XB,.Q') ) 2^4 ^1 



\Pcf(^^) Paa(^^) J 



g/(xi,Q2) 

G'(xi,g2) 

(38) 

How can one calculate Pgj and Pgg- The photon does not couple directly to gluons. 
The most straightforward way is to consider a longitudinal photon in fig. ^3. In this 
case, the a^ correction (|33|) does not contain a divergence and therefore does not 
contribute to the running of the quark densities. It is easy to understand intuitively, 
why there is no divergence for longitudinal photons. The divergence originates from 
configurations with p^ — ^ 0. This case however is exactly the naive parton model, 
in which o"/^ = due to helicity conservation. In order to obtain Pcf and Pgg-, one 
has to take into account at least one rung in the ladder. The splitting function is 
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the coefficient of the logarithmic divergence occuring at pxi — > 0. For this reason, 
longitudinal photons are often called " gluonometer" |50[ . 



At present the splitting functions are known to next to leading order in pQCD. 
They can e. g. be found in |51]]. The function Pgg{xbj/xi) ~ 6xi/ xbj (at low xbj) 
exhibits a pole at xbj = 0. It is widely believed that this pole is responsible for the 
steep rise of F2 and the high gluon density at low xbj, which behaves approximately 
like 

xbM^Bj, Q') oc exp ( 2v/ ^ Hl/XB,) HQ^Ql) J (39) 



for fixed strong coupling constant xbj- There are however also different opinions ||64 



65 



DGLAP is a special kind of a renormalization group equation. One of its most 
important properties is that it separates scales. All the soft physics is contained in 
the parton distributions. These are of entirely nonperturbative origin and have to 
be parametrized at some input scale Qq. Such parametrizations are provided by the 
three collaborations GRV, MRST and CTEQ ||4^, ^, ^ in leading and in next to 



leading order. With the parton distributions as input, one can then calculate F2 at 
a higher value of Q^, because the splitting functions can be calculated in pQCD. The 
other essential property is, that the parton distributions are universal, i. e. they do not 
depend on the process under consideration, but only on the hadron state. 

A more rigorous technique to separate scales in QCD is Wilson's operator product 



expansion |^. Using the OPE, the matrix of anomalous dimensions, i.e. the Mellin 



transforms of the splitting functions, was derived by Georgi and Politzer |66] and by 



Gross and Wilcek |]B^. Explaining the OPE is however not within the scope of this work 



and we rather refer to standard texts |5^. We only mention that in this method the 



time ordered product of currents in the virtual Compton amplitude (|T^ is expanded 
near the light cone in a series of local operators with coefficient functions which are 
singular on the light cone. These coefficient functions can be calculated in pQCD and 
are related to the splitting functions by Mellin transform. It is interesting to see within 
this framework, that (i) DIS is a long distance process receiving large contributions 
from the region near the lightcone, (ii) the partonic interpretation of F2 is valid only 
in the infinite momentum frame and (iii) additional terms occur, which cannot be 
interpreted as parton densities. These are the so called higher twist terms which are 
formally expressed as matrix elements of multiple field correlators. A typical example 
is a double scattering process. These matrix elements are, like the parton densities, not 
accessible by means of perturbation theory and have to be adjusted to experimental 
data. The matrix elements are however universal, since they depend only on the state 
in which the field correlator is evaluated. Higher twists are typically suppressed by 
additional powers of Q^, but can be enhanced by a factor A^^^ in processes involving 
nuclei, where A is the nuclear mass number. The A enhancement was argued by Luo, 
Qiu and Sterman [^, |69| . The calculation of higher twist terms is rather involved and 
they are theoretically not well under control. 
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2.3 The color dipole picture of low-x DIS 

In this section we have a closer look at the space time picture of DIS at low xbj- 
The two possible time orderings for the photon-quark interaction are depicted in fig. 
|7]. In [7|a, the virtual photon penetrates the target and hits a quark. In |7]b however, 
the 7* splits up into a qq pair that subsequently scatters off the target color field. 
In covariant Feynman-perturbation theory, both contributions are automatically taken 
into account. It is however instructive to consider DIS in the target rest frame and in 
a Hamiltonian, i.e. not manifestly covariant, picture, because this will give us a clearer 
physical understanding of the dominant mechanism and in turn we can make full use 
of our physical intuition. 





target 



target 



(a) 



time 



(b) 



Figure 7: The two time orderings for the interaction of a virtual photon with a target: 
(a) The photon hits a quark inside the target, (b) the photon splits up into a gg-pair 
which subsequently scatters off the target color field. In the rest frame of the target. 



the latter contribution is the dominant one. The figure is taken from ||TT 



The relative importance of the two contributions in fig. |^ may be estimated from 
the ratio of the energy denominators as it was done in [^ . For large photon energies, 
I' ^ rri]^, these energy denominators were found to be 



AEa = Ea{t2)-Ea{U) 

AEh = Eb{t2) - Ehih) 



-{plY" 



2v 






Q' 



2v 



(40) 
(41) 



where M^q is the invariant mass of the gg-dipole and (pi) is the mean square of the quark 
momentum inside the target. In order to estimate the ratio of the energy denominators, 
one usually sets M^^ —* Q^, since Q^ is the only large dimensionful scale available. This 
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was also argued in |7D[ with the result 



AE, 



AE,, 



ITiM 

2^B,(^y7^- (42) 



This expression makes very clear, that in the target rest frame, the process in fig. |^ 
is suppressed by a factor of order xbj <^ 1- For our physical intuition it is therefore 
sufficient to think of fig. |^. Indeed, we can estimate the lifetime of the qq fiuctuation 
with help of the uncertainty relation, 

1 1 

/^ = ^ ^ , (43) 

IXEh zrriNXBj 

where we again replaced M^^ by Q^. One recognizes that the lifetime or coherence time 
can become very long at low xbj-, e. g. 10^ fm at the lowest values of xsj accessible 
at HERA and 10 fm at NMC energies [xbj ~ 10"^). This is illustrated in fig. |l|, 
which corresponds to the Feynman diagram in fig. Pd, when only the gluon [xi^pTi) is 
taken into account. It is interesting to see, that although the 7*-proton cross section 
is a Lorentz invariant quantity, the partonic interpretation of the scattering process 
depends on the reference frame. 

The coherence length Ic is one of the key quantities in low-x DIS and we will present a 
more precise calculation than the prescription M^^ — »■ Q^ later in this work. Everything 
which happens to the photon within the length Zc, is governed by coherence effects. It is 
therefore essential to achieve a good physical understanding of these effects, especially 
in view of processes involving nuclei [|lT|, where multiple scatterings will occur. It is 



the main purpose of this work to provide insight into the mechanisms underlying these 
coherence effects. 

The cross sections for transverse and longitudinal photons are most conveniently 
written in a mixed representation. The two transverse directions are treated in coordi- 
nate space, while the longitudinal direction is described in momentum representation. 
Let p be the two dimensional vector pointing from the quark to the antiquark in the 
transverse plane and a the fraction of the photon energy v carried by the quark. The 
momentum fraction of the antiquark is then 1 — a, see fig. |I|. The cross section reads 

^T,L= da d^p '^^^^{a,p) agg{p), (44) 

where the "^g^ {a, p) are the light-cone (LC) wavefunctions for the transition 7* — > qq. 
The LC wavefunctions can be calculated in perturbation theory and read in first order 
in the fine structure constant aem 

ON ^^' 

l*I>'P)r = -^^ZJ{[l-2ail-a)]e^Kl{ep)+m}Kl{ep)}, (45) 



8Nra, 



N 



f 



l^>'P)r = ^^9l5?EW«'(i-«)'K^M' (46) 
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where Kq i are the modified Bessel functions of the second kind. These functions are 



e^ = a{l - a)Q^ + m), (47) 



also called MacDonald functions |7^. Furthermore we have introduced the extension 
parameter 

which depends on the fiavor mass rraj. The cross section for scattering a gg-dipole off 
the proton is denoted by (Tgg{p)- In Born approximation (two gluon exchange) it is 
independent of energy and related to the two-quark form factor of the proton via 

rr (n-) ^^"^ / ^2„ [1 " {P\ exp(ipri ■ (fj - r2))|P)] [1 - exp(ipTi ■ p)] , , ^ 

3 J Pti 

where pti is the transverse momentum of the exchanged gluon, fig. ^b. In Regge 
phenomenology, two gluon exchange is a simple model of pomeron exchange, cf. section 
|2.4] . Indeed, the dipole cross section takes only the pomeron contribution to the total 
cross section (0) into account. Therefore ( ^i)) can be applied only at high energies. 
Note the color screening factor [1 — exp(ip2"i ■ p)] in (^), which makes the dipole cross 
section vanish oc p^ at p -^ 0. This salient property of the dipole cross section is 
the heart of the color transparency phenomenon 0, ^, |^. In Born approximation, the 
two-quark form factor is related to the differential gluon density of the target by 

a(.,G(.,,.ph)) ^ AaAAA 11 _ (PI e.p(ip„ . (ft - f.))|P)| ^ (49) 

The energy dependence of the gluon density is caused by higher order QCD corrections, 
i.e. one has to take the gluon rungs in fig. ^d into account. The dipole cross section 
reads then |7^ 



— 2 (^| 2^ f ^Vl [^ " exp(ipTl " p}] djXBjGjXBj , P ti)) 
^P0.M/P)J ^^^ ^2^^^2 51og(p2^J 



(Tgg{XBj,p) = —p ".(A/p ) / —^ -^-^ L^„(J ^ ' (^O) 



with A = 2.25. For small distances p -^ one can expand the color screening factor in 
(0), which leads to 

71-2 

agg{xBj,p) = —p^as{\/p^) XBjG{xBj, A„/p^). (51) 



It was found in WM that A^r ^ 10. The dipole cross section still vanishes quadratically 



at small p, up to logarithms which originate from the gluon density. Note that the 



result of Nikolaev and Zakharov (|51|) coincides with the result from [|75|, ^ ^2. 



aqq{xBj,p) = yp as{X'/p ) XbjG{xbj,X' /p ), (52) 

except for the scale, at which the strong coupling constant enters. For a detailed 
derivation of (|5^ ) see |7^. The gluon density in ( ^2]) is tested at the same scale as the 



one in (0), i.e. A ^ 10 [0. 

For convenience, we do not write out the energy dependence of the dipole cross 



section explicitly, until section \2.1\ . The dipole cross section is largely unknown and 
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has to be parametrized. Several fits already exist in the literature |7^, ^, ^, ^], we 



will present our own parametrization in section \2.1\ . As already pointed out above, the 
space time picture of DIS depends on the reference frame. In the target rest frame 
the virtual photon splits up into a gg-pair which then scatters off the target. In the 
infinite momentum frame, the photon couples to the electric charge of the quarks in 
the target. Note, that there are valence and sea quarks in the target. Since the dipole 
cross section is proportional to the gluon density of the target, only quarks which are 
generated from gluon splitting are taken into account by (^^. In other words, the 
valence quark contribution to the cross section is neglected and therefore (|4^) is only 
applicable when sea quarks dominate, i.e. at low xbj- It is worth noting that it is 
impossible to decide whether a sea quark belongs to the target and is generated by 
gluon splitting or whether it is part of the hadronic structure of the virtual photon. 
This depends on the reference frame. 

It is important to note that while (^51) and (^) are valid only in perturbation 
theory, our formula (0) for the total cross section is completely general and does not 
rely on the applicability of pQCD. In the mixed p — a-representation, the scattering 
matrix is diagonal. Using the standard notation for the scattering matrix, see e.g. 



Sf, = 6f, + i{2nY6(''^ (y2Pf-Y.p)^. 



_ ^ fu (53) 

/ 

where the index / denotes the final state and i the initial state. One also often intro- 
duces the T-matrix element 



T,, = i2nr5^'^('£pf-J2p)^. 



_ , f^ (54) 

/ 

such that Sfi = 1 + iTjj. The unitarity of the S*- matrix leads to the optical theorem as 
a special case of the Cutkosky rules |8^ , 

2lmAii = Katou (55) 

where K k. 2s is the flux factor. The eigenstates of interaction are defined by 

5^lV'fc)(V'fc| = 1, (56) 

k 

T\iJk) = hliJk). (57) 

Here, tk is the probability that the eigenstate \ipk) scatters off the target. The photon 
can now be expanded in this basis [H, ^, 

k 

with 

(7*|7*) = 1. (59) 
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For simplicity, we do not consider wavefunction renormalization and assume that alss 
states are normalized to unity. Thus one obtains 



E 



4 ]% = sa-- '. (60) 



Comparing this expression with (0), we identify the eigenstates of the interaction 
with color dipoles with fixed transverse separation and fixed longitudinal momentum 
fractions. The index k labels the (p, a) of the dipoles. Furthermore, we can identify 
tk/s with the cross section of a dipole with a given p off the target, i. e. with aqq{p), 

^M-(P) = -J- (61) 

The identification of interaction eigenstates with color dipoles of fixed transverse sepa- 
ration in QCD was first made in 0]. There are of course also higher Fock-states in the 
photon, for example the qqG state will play an important role for diffraction in DIS. 
This will be discussed in section |2.4] . 

It is a special virtue of the color dipole description of DIS to work with the correct 
degrees of freedom, which are the eigenstates of the interaction. This will be of great 
use, when we consider nuclear targets in the section ^. 

The p-integral in ( ^4|) extends from small distances, where pQCD is applicable, up 
to very large distances, which are governed by infrared properties of QCD that are 



not understood at present. The extension parameter (|47|) controls which distances p 
in (|4^ ) contribute to the integral. Note that the asymptotic behavior of the Bessel 
functions is 

Ko{z) = -\og{z) for ^^0, (62) 

Ki{z) = - for z^O, (63) 

z 

Ko(z) = KUz) = J — e-' for z -^ oo. (64) 

V 2z 

Both functions, Kq and Ki, decay exponentially for large arguments and are divergent 
at 2 — >■ 0. While Kq has only an integrable, logarithmic divergence, Ki diverges much 
stronger and the wavefunction is not square integrable. Numerically, the Ki-part is the 
dominant term. 



Following [13], we will now investigate, how these hard and soft contributions inter- 
play. Special attention is paid to the question, which regions of integration give scaling 
contributions to F2 and which are higher twist, i.e. are suppressed by additional powers 

ofg2. 

First consider transverse photons and estimate is the probability w, that the photon 
splits into a small size pair p^ ~ l/<5^- This probability is of course practically 1, 
because the LC wavefunction (|i5| ) decays exponentially at large ep. The largest dipoles 
that can contribute have a size of order 1/e. A small size dipole interacts with a 
cross section roughly proportional to its size, because of color transparency. As result, 
one finds that the contribution of small size dipoles from transverse photons give a 
contribution to a'^ ^ that decays as 1/Q^ and thus a leading twist contribution to F2. 
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hard 
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Table 1: The Q^ behavior of hard and soft contributions to the total photon-proton 
cross section for transverse photons. The probability that the photon converts into a 
small (hard) or large (soft) dipole is denoted by w, a is the order of magnitude of the 
dipole cross section and the last column contains the Q^-behavior of the total photon- 
proton cross section. (Jh is a typical hadronic cross section, i.e. independent of Q^ . For 
transverse photons both, small and large dipoles, give leading twist contributions to 

F2. 



How can the virtual photon convert into a large size dipole? The only possibility 

is, that the extension parameter in the Bessel function (^5]) e becomes small, i.e. one 

needs a — »• or a — > 1. These limits are known as Bjorken's aligned jet configurations 

8^. The phase space for these highly asymetric configurations is however very small. 



of order m'i/Q'^, and thus the probability of finding a large dipole is also small. This is 
however compensated by the large interaction cross section, which is of the order of a 
typical hadronic cross section, ah- One therefore obtains a leading twist contribution 
to F2 from large dipoles that cannot be treated perturbatively. These considerations 
are summarized in tab. |l|. 

For longitudinal photons and for small dipoles, the situation looks like in the trans- 
verse case and the hard part gives a leading twist contribution to F2. This is however 
different in the case of large dipoles. Because of the term a^(l — a)^ in the longitudinal 
LC wavefunction (0), aligned jet configurations are suppressed by an extra power of Q^ 
compared to the transverse case. There are no large dipoles from longitudinal photons. 
Also the large cross section an cannot compensate for this and the soft contribution to 
F2 is higher twist. This is summarized in tab. |[ 

Note that the integral for the transverse cross section in ( PI) is logarithmically 
divergent (for a quadratically rising dipole cross section), if one would not include the 
quark mass mj. This divergence originates from the end points of the integration over 
a and is the analog of the logarithmic divergence in (0). Again the reason is, that the 
integral extends over domains, where pQCD is not applicable. In contrast to this, the 
longitudinal cross section is finite, even with mj = 0. No collinear divergence occurs 
in this case, because the endpoints of the a-integration correspond to the naive parton 
model in which aj ^ = 0. Therefore, no divergence can come from these configurations. 
We point out, that there is a one-to-one correspondence between the divergences in 
the standard pQCD in momentum space and the color dipole picture of DIS, which is 
formulated in the mixed representation. Factorization properties in impact parameter 
space are also considered in M . 
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Table 2: The same as tab. but now for longitudinal photons. In the longitudinal 
case, only small (hard) dipoles give a leading twist contributions to F2. There are no 
large dipoles from longitudinal photons. 



For a photon, the LC wave functions, i.e. the coefficients c^ in (pBD, can be cal- 
culated in perturbation theory. For an introduction to LC wave functions and how 
to use them, see |^. Although they can be calculated in a manifest covariant way, 
the situation is more transparent in a Hamiltonian framework. One has in first order 
pQED 

iZf,/a^ u{pg)eT,L ■ av{pq) 

_ 1/1/ I im I 

qq 



^!l'J'{a,pT) 



47ra(l — a)u 



Pt- 



iO. 



(65) 



2a{l-a)i/ 

where u{pg) and v{pg) are Dirac spinors of the quark and the antiquark, respectively 
and d is the three-vector of Dirac matrices, not to be confused with the longitudinal 
momentum fraction a. Furthermore, eTis the polarization vector of the photon. The 
energy-denominator in (p5| ) depends on mass Mgq of the pair. In the mixed represen- 
tation, one can express the LC wave function in terms of the Green function for the 
propagation of the pair in vacuum, integrated over time, or over the z-coordinate which 
is the same. 



an IP2, '^2\P\, Zl) 



{27: f 



exp I 1 



pl + e^ 
2a(l — a)v 



iO, 



(22-^1) . (66) 



Integrating this expression over the 2;-difference, Az = Z2 — zi, one obviously recovers 
the energy-denominator. 



/ 



° C (Pvr Q-^PT-ih-pi) 

dAzG-{p2,z2\puZi) = -ly 7^ ,.^.. 

2a(l~a)i/ ^ 



-i2a(l - a)uKo{6{p2 - pi))- 

(67) 



In the last step the relation 



Koiep) 



(fpT- 



-ipT-p 



2tx J " '"" pl + e"^ 



(6J 



was used. Finally, the LC wavefunctions can be written as 



,T,L 



^qh ("'P) 



47ra(l — Ci)v 



Xq^qq Xq 



rfAzG™'=(p,^2|0,^i) 



(69) 
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where the x ^^re two component spinors and the operators Og^ are defined by 

d^g = mja ■ e + i{l -2a){a-n){e-Vp) + {a xe)-Vp (70) 

dj;^ = 2Qa{l - a)a ■ n. (71) 



The two dimensional gradient Vp acts on the transverse coordinate p, n is the unit 
vector parallel to the photon momentum and a is the three vector of the Pauli spin- 
matrices. For a Green function G""^^ without interaction, one obtains of course the 



perturbative LC wavefunctions on page |T^. To see this, one simply notes that 

d 



K,iz) 



dz 



Ko(^) 



(72) 



Introducing a quark mass which acts like a cutoff and prevents the integrals from 
diverging, is obviously a poor model of the nonperturbative effects at large p. A 



better solution was suggested in |jT2[, where an interaction between the quark and the 
antiquark was explicitly introduced. For this purpose, the Green function in (^) was 
modified by introducing a harmonic oscillator potential. 



K?(«' P) 



a] 



2ua{l 
The Green function including the nonperturbative interaction reads 



(73) 



Gr(P2,^2;Pi,;^i) 



a? ia] 



•qq 



27rsinh (cj Az) 



exp 



2vail 



a] 



X exp 



where 



a^ ia) 



UJ 



(Pi + Pt) coth {uj Az) : 



2pi ■ p2 



sinh {u Az) 



a a 



,(74) 



(75) 



V a[l — a) 

is the oscillator frequency. It was found in ||12| that data for diffraction dissociation 
and the photoabsorption are well reproduced with the parametrization 



a^ ia] 



,1.15 



(112MeV)^ 



,1.15 



(165MeV)^a(l 



a] 



(76) 



The first term in this ansatz prevents that the transverse distance becomes arbitrarily 
large at the endpoint a — »■ 0, 1. Note however that the more conventional ansatz, which 



results from a relativistic approach to the gg-bound state problem p6[ is setting the 
parameter f = 0. Numerical results are surprisingly insensitive to the value of the 
parameter w, which could not be fixed in [|12|. We will always use the value v = 0.5. 

The nonperturbative interaction changes of course the LC wavefunctions, which 
read now 



|*7*( 



a,p)\ 



Nf 

2Ncaem sr^ 



i2ny 



z}{[^ 



2a{l-a)]^l{e,\p) 



m)^l{e,\p)] 



i^r^p)! 






Nf 



5^ZjgV(l-a)2$g(e,A,p), 



(77) 
(78) 
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with 



<^o{e,\,p) 



du 



P 



$1 (e, A,p) = — I duexp 



X 




sinh (An) 




r AeV' 




4 



2^2 



Xe'p 



coth (An) 



n 



coth (Am) 



M 



(79) 



^0) 



The strength of the nonperturbative interaction is described by the dimensionless pa- 
rameter 

A = ^^. (81) 



Compared to the expression for $i in |T2[, we have integrated by parts over the pa- 
rameter u. This considerably simphfies the expression. 

It is worthwhile to investigate the two limiting cases of vanishing interaction, A — > 0, 
and of strong interaction, A — > oo. With the useful relations [^ 



Ko(ep) 



epKi{ep) 



oo 

1 r du 



2 J u 



oo 



exp 



4u 



u 



duexp 



4n 



u 



f82) 



(83) 



one recognizes that in the limit A ^ 0, where the nonperturbative interaction becomes 
negligible, the perturbative LC wavef unctions, (^5]) and (^6]), are recovered. In the 
strong interaction limit, A -^ oo, the functions $o,i again acquire simple forms [O, 



%{e,X,p) 
$i(e,A,p) 



A— »oo 



A— ►oo 



47r " I 2 



p I a'^(a)p'^ 

-exp' 



271 p- 



B5) 



This limit is appropriate in the case of real photons. The nonperturbative interaction 
confines even massless quarks. However, we will include current quark masses in all 
our calculations. 

We point out that our approach to describe the soft contributions to DIS is quite 
different from parametrizing all the unknown physics into parton distributions. Instead 
we develop a model, that smoothly interpolates between the hard and the soft part of 
the interaction. It will turn out later that this way of treating the soft physics allows 
us to develop a better physical understanding of coherence effects in QCD. 



2.4 Diffraction 

Diffraction in hadron-hadron collisions has been known for a long time. A process is 
called diffractive, if only the quantum numbers of the vacuum are exchanged in the 
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t-channel pB|, 53 . The best example for such a process is elastic scattering, A + B - 



A + B. In the language of Regge phenomenology, the trajectory with the quantum 
numbers of the vacuum is called Pomeron trajectory, because it was first proposed 
by Pomeranchuk [^. Elastic scattering is not the only type of diffractive processes. 



There is also inelastic diffraction, where one or both colliding particles are excited into 
states with the same quantum numbers as the incoming particles. Since no quantum 
numbers are exchanged, the observation of a large rapidity gap between the outgoing 
particles is characteristic for a diffractive event and serves as experimental criterion of 
diffraction. 

An intuitive picture of diffraction was proposed by Good and Walker [^, intro- 
ducing an analogy between diffraction and wave optics. The beam particle can be 
decomposed into a coherent sum of interaction eigenstates. These eigenstates scatter 
with different amplitudes off the target and the coherence is destroyed. New particles 
are produced in the same way, as white light is decomposed into different colors, by 
sending it through a prism. The term " diffraction" originates from this analogy. 

Typical properties of total hadronic cross sections are among others, (i) a slow rise 
of the total and the elastic cross section with energy, (ii) a large imaginary part of the 
forward scattering amplitude, compared to the real part and (iii) a small elastic cross 
section, compared to the total cross section. One can argue from these observations 
that the theory of strong interaction must be nonabelian. Indeed, if only one particle 
exchange is considered in an abelian theory one would have elastic scattering in the 
first order of the coupling constant and therefore a large elastic cross section and a 
predominantly real forward scattering amplitude ||90 . 



The Pomeranchuk theorem ||9^, |93[ states that any scattering process which involves 
charge exchange, must vanish at asymptotically high energies. The converse is also true. 
Foldy and Peierls MM have proven that, if a cross section does not vanish asymptotically. 



the process must be dominated by the exchange of vacuum quantum numbers. A 



simple model of the Pomeron has been suggested by Low ||^ and by Nussinov [^, |^ . 
The basic idea is that elastic hadron-hadron scattering is mediated by exchange of two 
gluons in a color singlet state. This model qualitatively explains the basic experimental 
observations. One obtains a purely imaginary forward scattering amplitude and the 
total cross section is independent of energy. It is argued that the observed increase 
with energy, atot oc s°°^, can be attributed to higher order corrections. 

It came as a big surprise, when large rapidity gap events were discovered in DIS at 
HERA. Naively one would expect that the virtual photon hits a quark inside the proton 
and produces a jet. Due to the strong color forces, hadronic activity is expected in the 
whole rapidity region between the jet and the proton remnants. However, about 10% 
of all DIS events at xbj < 0.1 show a large rapidity gap between an almost elastically 
scattered proton and a diffractively excited state. The amount of these events shows 
only a weak Q"^ dependence, suggesting that diffraction in DIS is a leading twist. 
The energy dependence is however much stronger than in the case of hadron-hadron 
scattering, approximately a"'*^ oc a;^°'^, which lead to the assumption that there are 
two Pomerons, see e. g. ||6^, ^. The so called soft Pomeron shows only a weak energy 
dependence and is responsible for the rise of total hadronic cross sections with energy, 
while the hard Pomeron is observed in the stronger energy dependence in DIS. 
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Figure 8: Diffraction in DIS seen 
in the target rest frame, cf. fig. |l|. 
Only tlie quantum numbers of the 
vacuum are exchanged in the t- 
channeL This can be modeled by 
the exchange of two gluons in a 
color singlet state. The proton re- 
mains almost intact and a large ra- 
pidity gap is observed between the 
outgoing proton and the diffrac- 
tively excited state qq. 



A popular way to think of diffraction is to imagine that the virtual photon scatters 
off a preformed color neutral cluster inside the proton, the Pomeron. This picture is 
known as Ingelman-Schlein model ^^. The space-time picture of a diffractive event 
is shown in fig. ^ where the Pomeron is represented by two gluons. The hadronic 
fluctuation of the photon is developed long before the target. 

Diffraction in DIS requires the introduction of additional kinematical variables. The 
quantity 



X-p 



Ml 



Q' 



6) 



may be regarded as the momentum fraction of the proton carried by the Pomeron. Here, 
Mx is the mass of the diffractively excited state. One also often uses the variable 



/5 



Q' 



Q^ + Ml 



XBl 

x-p 



(87) 



instead, which can be interpreted as the momentum fraction of the struck quark relative 
to the Pomeron. Furthermore, if the proton is observed in the final state, one also has 
to deal with the four momentum transfer squared at the proton vertex, 

t = {P- P'f < 0, (88) 

where P and P' are the four momenta of the incoming and outgoing proton, respec- 
tively. 

Thinking about diffraction as DIS off the Pomeron allows to introduce parton densi- 
ties for the Pomeron and to evolve these densities with the DGLAP evolution equations. 
A fit to these parton densities can be found in [ p.0(j| , p.01[| . It is however only a postulate 
that diffractive parton distributions evolve according to DGLAP. It is not even clear, 
whether these distributions have to fulfill a momentum sumrule. Indeed, in ||100|| the 
momentum of all partons together is significantly larger than 100% and the parton 
densities become even negative at small (3. 

The color dipole picture of diffraction does however not rely on the assumption that 
the 7* scatters off the pomeron. The decomposition of the virtual photon in interaction 
eigenstates (E3) is completely independent of this picture. Furthermore, we can not 
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only decompose the photon, but any hadron in such a series. However, in the photon 

* 

case the LC wavef unctions, i. e. the coefficients c^ , can be calculated perturbatively 
at least at small p. This is not the case for hadron wavefunctions. Proceeding like on 



pp. |T9|, we obtain for the mass integrated diffraction dissociation cross section 



dM^-—^- 
^ dMldt 



t=o 



IGtts^ dt \t=o 


1 


E i-r X - (y: let %] 
k \ k ) 



(89) 
(90) 



where we have subtracted the elastic part. This part is however of order Og^ and will 
be dropped from now on. Keeping only the qq Fock-component of the photon, one 
obtains 

= Y^J daJd^p\^,gia,p,Q^)\^(Tl^{p), (91) 



dMl-—^- 



^x^ 



where \Ef denotes either the perturbative ( ^ , |i6| ) or the nonperturbative (^, |78|) LC 
wavefunctions of the virtual photon. We omit the indices T and L from now on and 
assume summation over all polarizations. 

It is instructive to consider also off-diagonal diffraction, where the photon goes into 
a final hadronic state h. Decomposing the hadron into interaction eigenstates. 



\h) = Y.^km^ 



(92) 



one obtains 



da^^iY -^ h) 



dt 



t=o 



1 
1 



167r 



{h\T\Y) 



E(-9*-rt. 



d^p da'^lia, p)^gg-(a, pWqJp) 



(93) 
(94) 

(95) 



where "^h is the LC wavefunction of the hadron. It is convenient to express this result 
in terms of the diffractive amplitude 



/(7* 
da^^iY -^ h) 



dt 



qq) =i^qq{a,p)aqq{p), 

' if{h-^qq)\f{Y^qq)). 



t=o 



IGir 



(96) 
(97) 



We see from (p^ ) that the photon cannot go into an orthogonal state, if all its eigencom- 
ponents would scatter off the target with the same amplitude tj. = s (Tgq{p). Offdiagonal 
diffiaction occurs, because the eigencomponents scatter with different amplitudes and 



thus the coherence between them is disturbed. This is exactly the picture of |95 



We mention that the photon can of course also be represented as a superposition of 
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T 


w 


a^ 


iO 


hard 


1 


1 


^Q^ 


soft 




-I 


^Q^ 



Table 3: The Q^ behavior of hard and soft contributions to the diffractive cross section 
in DIS for transverse photons, cf. tab. |1]. The probabihty that the photon converts into 
a small (hard) or large (soft) dipole is denoted by w, a^ is the order of magnitude 
of the dipole cross section squared and the last column contains the Q^-behavior of 
the diffractive cross section, ah is a typical hadronic cross section, i.e. independent 
of Q^. The probabilities for finding small or large dipoles is the same as in DIS. The 
contribution from small dipoles is higher twist in diffraction, the leading twist originates 
from large, nonperturbative dipoles. 



hadronic states with the same quantum numbers, i. e. vector mesons, rather than in- 
teraction eigenstates. This is done in the vector dominance model (VDM). A review 
article can be found in 11102 . 



We will now examine, like in the preceding section, from which values of p the 
leading twist contribution to diffraction originates. This was also done in ^^. Since 



the LC wavefunctions are the same as in DIS, the probabilities w to find small or large 
dipoles does not change, when one considers diffraction. Note however we have now 
<Jgq{p) in the integrand (pl|). As a consequence, the contribution of small size dipoles 
is suppressed by an extra power of Q^, see tabs. H and §. For transverse photons, the 
leading twist originates purely from soft contributions, i. e. from large dipole sizes, 
making diffraction a phenomenon dominated by nonperturbative physics, in spite of 
the large Q^. Furthermore, we have seen that in the case of longitudinal photons, the 
only scaling contribution in DIS comes from small size dipoles, see tab. 0. This is no 
longer the case in diffraction and the diffractive cross section for longitudinal photons 
seems to fall away as oc 1/Q^. This is however not the case, it will be shown later, that 
the qqG Fock state gives a contribution oc l/Q^- 

Consider the M|- dependence of diffraction for large M|-. First one notices that a 
gg-pair with fixed transverse separation does not have a well defined mass, because the 
transverse momenta are completely undetermined in coordinate space representation. 
In momentum space representation one has 



Ml 



p^ + rni 
a(l — a) 



(9J 



We have seen above that the diffractive cross section for transverse photons is domi- 
nated by contributions from the endpoints, a ^ and a -^ 1. These are large size 
fluctuations, which have small transverse momenta. The mass is therefore approxi- 
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a^ 


K) 


hard 
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^Q^ 


soft 




-I 


^Q^ 



Table 4: The same as tab. |^ but for longitudinal photons. In the longitudinal case 
both, hard and soft, contributions to the diffraction dissociation cross section are higher 
twist. The leading twist contribution conies from the qqG Fock state of the photon, 
see fig. p. 



mately Mj- ^ mi /a and thus 



dMl 
-Ml 



da 



a 



In the large mass limit, one can write |72 



dM^dt 



f 
-o°^Ml 



rfVlC(«,p)l\\(p). 



(99) 



(100) 



Here we consider the limit a — > 0. The other endpoint gives the same contribution. 
Since large dipoles interact with a typical hadronic cross section ah, one sees that the 
diffr active mass spectrum decays as oc 1/Mx- This behavior is a result of the small 
phase space for the aligned jet configurations and is at variance with the experimental 
observation that the mass spectrum decays only as ex l/M^-. 

The solution to this problem is that one of the quarks can radiate a gluon before 
the impact on the target, fig. ^ Then the multi parton configuration, frozen in impact 
parameter space scatters off the target. Although the gluon is radiated on the cost of 
an additional power in the strong coupling constant as, the spin-1 nature of the gluon 
leads to a much weaker decay of the diffractive mass spectrum and the longitudinal 
cross section for diffraction gets a leading twist contribution. The LC-wavef unction for 
the transition 7* — >■ qqG reads in the small recoil approximation, ac ^ 0, [|12| 



"^qqcia, "G, Pi, P2) = ^^(a, Pi - P2] 



V a 



m 



qG 



Qg 
1 — a' 



P2 



(101) 



where ^ qc is the LC wavefunction for the transition q -^ qG. Here, a is the longitudinal 
momentum fraction of the photon carried by the quark and ac is the momentum 
fraction of the gluon with respect to the photon. Furthermore, pi is the distance 
between the quark and the gluon and p2 the distance between the antiquark and the 
gluon. The first term in the square brackets corresponds to radiation of the gluon 
from the quark, while in the second term the gluon is radiated from the antiquark. In 
analogy to (|69|) , ^qc is defined as 



%G(tt,P) 



27rQ;(l — a)u 



Xq^Xl, 



d{^z)Gll?{p,Z2\^,z,) 



(102) 
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Figure 9: The same as fig. || but 
now for tlie qqG Fock-component 
of the virtual photon. For longitu- 
dinal photons, a leading twist con- 
tribution to diffraction arises from 
configurations, where the gluon is 
soft and far away from the gg-pair. 



where the operator 



O, 



qG 



mifa^e ■ {n x a) — i(2 — a){e- Vp) + ae- {a x V^ 



(103) 



comes again from the spinor structure of the vertex. Without a potential between the 
quark and the gluon, the integral over the Green function reads 



d/:^zGl^{p, Z2|0, zi) = -i2ua{l - a)Ko(rp), 



where 



2^2 



T = a m 



f 



(104) 



(105) 



plays now the role of the extension parameter. The LC wavefunction reads explicitly 



^gG(a,Pi)^5G(a,P2) 



IQa. 



3 



[l + (l-af] 



21 2P1 ■ P2 



P1P2 



Ki(rpi)Ki(rp2) 



+ m^ra^Ko(rpi)Ko(rp2) 



(106) 



The contribution of the qqG Fock state to the diffractive cross section is given by 






dt 



167r 



da 



dac 



2 ,2 



d pid p2\'^qqGia,aG,pi,p2)\ agg-G(pi,P2). (107) 



t=o -^" JO J "G 

The cross section for scattering the qqG state off the proton can be expressed in terms 



of the dipole cross section ||103| , |104 | 



9 

C^W-g(P1, P2) = - \^qq{p\) + (yqq{p2) " (^ qq{p\ " P2)] 



(108) 
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The solution to the problems mentioned above lies in the integral over dac/o^G- The 
singular behavior at ac — ^ is characteristic for the radiation of vector bosons. Note 
the phase space in ac is no longer suppressed as ac — > 0. Therefore, the gluon can go 
far away from the quark and the qqG system scatters with a large cross section. This 
makes diffraction for longitudinal photons a leading twist effect and also leads to the 
asymptotic behavior 

dMldt ,^0 ^ Ml ^ ^ 

for Ml ^ Q^, as observed in experiment. In the language of Regge theory [|105|| , the 

qqG component corresponds to the triple Pomeron vertex, which dominates the large 

mass behavior. In the language of pQCD, the additional gluon gives us the evolution 

of the gluon density. We explained in section |2.2| why it is necessary in DIS to take at 

least one gluon-rung in the ladder in fig. ^ into account to obtain the gluon splitting 

function Pgg- 

There are two interesting limiting cases of (|108|) . First, the separation of the qq 

dipole is very small, of order l/Q^. This situation is typical for longitudinal photons. 

Then we obtain 

9 

(7qqG{pl = P2) = ^^gg-(p) = C^Gg(p), (HO) 

the cross section reduces to the cross section for an octet-octet, or gluon-gluon, dipole. 
Indeed, after radiation of the gluon, the qq dipole is with a high probability in an octet 
state, as one can see from color algebra. If the size of this dipole is small, it appears 
like a gluon and we have an effective gluon-gluon dipole. The factor 9/4 is the ratio of 
the two Casimir operators of the adjoint and the fundamental representation of SU(3). 
In the opposite limit, where the gluon stays close to one of the quarks and the qq dipole 
has a large separation, one obtains 

(^qqG (Pl ,p2 = 0)=aqq{pi). (Ill) 

This limit corresponds to low virtualities Q^. 

In the last section, it was demonstrated that the quark mass acts like an effective 
cutoff, rendering all integrals finite. Also in the case of the qqG Fock state, the integrals 
would diverge, as the gluon gets infinitely far from the quark it has been radiated from. 
We do not want to introduce a gluon mass, instead we follow [jl2| , where an interaction 



between the quark and the gluon was proposed, which is again described by a harmonic 
oscillator potential. 

The ansatz for the parameter that describes the strength of the interaction is chosen 
like in the qq case, 

b'^{a)=bl + 4:bla{l-a). (113) 

However, since the dominant contributions come from the region ac -^ 0, only 

6(0) = bo = 650MeV, (114) 
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was fixed in [|I2] with help of a triple pomeron analysis of CDF data |106| |. For this 
reason, we will always use the limit 



*J(p) ^ Jim *-(a..p1 ^ lyi^^exp (-%!) (115) 



in our calculations. We point out, that the qualitative corrections due to the additional 
gluon are not a result of the ansatz for the nonperturbative interaction, but are purely 
due to the spin 1 nature of the gluon. 

It is worth noting that the relatively large value of 60 does not allow that the distance 
between the gluon and the quark becomes as large as 1 fm. The maximum distance is 
rather about 0.3 fm, considerably smaller than a typical confinement radius. There are 
several phenomenological indications that this must be the case. The small size of the 
gluonic fluctuation is related to the high octet-octet string tension, K[8] fa 4 GeV/fm, 
compared to the triplet string tension, ^[3] fa 1 GeV/fm. These string tensions are 
related to the slopes of the pomeron and reggeon trajectories respectively ||105|| . If 



one thinks of diffraction in proton-proton scattering as proton-pomeron scattering, 
one can extract the proton-pomeron total cross section from CDF data. This cross 
section is about one order of magnitude smaller than the total pp cross section. The 
only explanation is that the pomeron is a small size object compared to the proton 
|107|| . The conjecture that the pomeron has a rather small size is also supported by the 



approximately fulfilled quark counting rules [|108|| . Indeed, the pomeron seems to couple 
to the single valence quarks inside a hadron. A large pomeron would couple to the whole 
hadron at once. Finally, we mention that a rather short gluon correlation length, about 



0.3 fm, also results from lattice calculations ||109|| . These lattice calculations are used 
to determine the values of the free parameters in the stochastic vacuum model of Dosch 
and Simonov | 110 , |1 1 1|| . Calculations within the SVM for soft diffraction are in good 



agreement with experimental data. 

2.5 The parton model of the DY process 

The concepts of the parton model, originally invented for DIS, can also be applied to 
certain processes in hadron-hadron collisions. The most prominent example for this is 
the Drell-Yan process ||112|| , where vector bosons are created in hadronic collisions. We 
will consider only photons. The model of Drell and Yan is depicted in fig. [1^. Two 
hadrons collide and a quark from one hadron annihilates with an antiquark from the 
other hadron into a timelike photon. This photon decays into a lepton pair that can 
be detected. 

We denote the mass of the spacelike photon by 

M^ = q^>0, (116) 

where q'^ is the four momentum of the virtual photon. The square of the center of mass 
energy of the colliding hadrons is 

^=(Pi + P2)', (117) 
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Figure 10: The DY process to leading order. Two hadrons collide and a quark from 
one hadron annihilates with an antiquark from the other hadron into a timelike photon 
with mass M^ > 0. This photon decays into a lepton pair that can be detected. 

where Pi and P2 are the four momenta of hadron 1 and hadron 2, respectively. A 
convenient variable to work with is Feynman xp, 

Xf = -^ ^Xi-X2. (118) 

Here, pf^ is the longitudinal momentum of the dilepton in the hadron-hadron center 
of mass frame and xi and X2 are given by 

2P2 ■ g 2Pi-q , , 

xi = —g— ' 2:2 = ^ ■ (119) 

These variables have the meaning of the longitudinal momentum fractions of the quarks 
participating in the hard process. The quark q in fig. |l^ has longitudinal momentum 
xiPi and the antiquark X2-P2- It holds 

xia;2 = — , (120) 

where the transverse momentum of the virtual photon is zero in the naive parton model 
and has been neglected. Another frequently used variable is 

-^. (121) 

Similar structure functions as in DIS can be introduced for the DY process, starting 
directly from the hadronic tensor, 

W"^ = I d^xe''i%PiP2\J^'{x)r{Q)\PiP2), (122) 
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as we did in section pTTI. Since there are two hadrons involved in DY, there are more 



possible Lorentz structures than in DIS. One can define four independent structure 
functions for DY, rather than two in DIS. We do not pursue this approach in detail 



and refer only to the literature [ p.l3|| . 



The partonic cross section for fig. 10 reads 



da _^-^lmZ^^,^,^s-M-). (123) 



The factor A'c, number of colors, appears in the denominator, because quark and 
antiquark must have the same color in order to annihilate. Embedding the partonic 
cross section into the hadronic environment yields 

-^ = / dx^dx^ Y. {'?/(^i)^/(^2) + (1 -- 2)} -^, (124) 

Jo ^ 

where qf{xi) is the probability to find a quark of fiavor / with longitudinal momentum 
fraction xi in hadron a and g/ is the analog for antiquarks. The (5-function in (|1 23|) 
allows to write the cross section in the scaling form 

M'^ = ^^ f dx, J2 Zj {qf{x,)qf{T/x,) + (1^2)}. (125) 

The r. h. s. of (|125|) depends only on r and not separately on M^ and S. 



The observation of this scaling property in experiment, see e. g. |p.l4| , |115|| , confirms 
that the mechanism illustrated in fig. ^ is correct. There are however features of 
dilepton production which cannot be understood in the lowest order picture. 



• Cross sections calculated straightforward from ( 125 ) are too small by a factor of 
2 — 3 compared to the measured value. This discrepancy is usually treated by 
introducing a so called i^-factor. The K-iactor is approximately independent of 
M^, but it is process dependent. 

• Large transverse photon momenta, of order few GeV, are observed in experi- 
ment. There are however no transverse momenta in the naive parton model. 
Phenomenologically, one can introduce a primordial momentum distribution of 
the quarks. One usually assumes a gaussian shape for this distribution, but the 
width necessary to describe what is observed in experiment is much larger than 
what one would expect from Fermi motion. 

These problems can be overcome by taking into account the first order QCD cor- 



rection, shown in fig. |TT| The first row contains the virtual corrections to the quark 
propagator and to the vertex. The second row shows the so called annihilation process, 
where the quark or the antiquark radiates a gluon before it annihilates with a parton 
from the other hadron. Due to the radiation of the gluon, the quark acquires a trans- 
verse momentum. In this way, the pQCD correction provides the missing mechanism 
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Figure 11: Higher order QCD 
corrections to the DY process. 
The upper row contains virtual 
corrections. The diagrams for the 
annihilation process are shown in 
the middle, where the quark or the 
antiquark can radiate a gluon be- 
fore the annihilation process. The 
QCD Compton process is depicted 
in the last row. Here, the pro- 
jectile quark scatters off a gluon 
from the target and radiates a 
massive photon. The Compton 
process is the dominant contribu- 
tion. These higher order correc- 
tions account for most of the K- 
factor and explain the occurance 
of large transverse momenta. The 
figure is taken from |[116| . 



for the production of lepton pairs with large transverse momentum px- This explana- 
tion was suggested first in ||117| , |118| , |119|| . The last row displays the diagrams for the 
QCD Compton process, where a quark in one hadron picks up a gluon from the other 



hadron and radiates a photon. This mechanism is dominant at large px |52 



The graphs in fig. |Tl| contain of course divergencies. While the infrared divergencies 
cancel in the sum of virtual and real corrections, the coUinear divergencies, which occur 
e. g. when the transverse momentum p^ of the intermediate quark in the QCD Compton 
process goes to 0, turn out to be identical to the divergencies observed in DIS, cf. section 
Thus they can be absorbed into a redefinition of parton densities. For example. 



the QCD Compton process gives the additional contribution [p.l8|| 



M^ 



d^ 



c 



1^/ ^/ ^E^HW-)+./M)Gfe)H-(l«2)} 



X 9(a;iX2 — r) 



27r' 



fG 
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a;iX2 






K^ 



(126) 



where G is the step function and Pfc is given by (Q). It was first noticed in 

that the whole DY cross section, including (|126| ) can be recasted in the form (|125|) by 

redefining the quark density. 



qf{x) 



qf[x) 



2-K 



In 



M2 



Uy 



P 



fG 



X 



G{y). 



(127) 



y \y. 

This is exactly what was done for DIS in section |2.2| . Thus the redefined parton 
densities in DY obey the same renormalization group equation as the DIS parton 
densities, namely the DGLAP equations (^). 
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The first order correction solves most of tlie problems of tlie naive parton model. 
It explains how large pT-dileptons are produced and account for almost all of the 
ii'-factor ||116|| . However, not all problems are solved by this correction. Since it is 
numerically large, one has to investigate how much higher order corrections change 
the result. Furthermore, the transverse momentum spectrum is even qualitatively not 
well described. The theoretical result agrees only at p^ ~ M^ with the data and even 
diverges at p^ — ^ 0, 

da QisiPT) 



dpT Pt 



128 



while the experimental result is of course finite. The reason for this behavior is that 
large logarithms ln(M^/p|.) occur in higher order corrections and one has to resum all 
these terms. This is possible within pQCD ||12CI| , p.21| , p.22|| , by a resummation of soft 



gluons radiated from the quark or the antiquark, respectively. The result indicates that 
one needs almost no intrinsic transverse momentum and practically all pt is generated 
perturbatively ||123 . 



2.6 Light-cone approach to DY 

Although cross sections are Lorentz invariant, the partonic interpretation of the mi- 
croscopic process depends on the reference frame. We have seen this in section ^]3| for 
the case of DIS and similar considerations can be made for DY. It was pointed out 
by Kopeliovich ||13| that in the target rest frame, DY dilepton production looks like 



bremsstrahlung, rather than parton annihilation. The space-time picture of the DY 
process in the target rest frame is illustrated in fig. ^. A quark (or an antiquark) from 
the projectile hadron radiates a virtual photon on impact on the target. The radiation 
can occur before or after the quark scatters off the target. Only the latter case is shown 
in fig. |. 

The cross section for radiation of a virtual photon from a quark after scattering on 
a proton, can be written in factorized light-cone form [|13|, |14|, |1^ , 



j^ = d p\'$^lg{a,p)\ agq{ap), (129) 

similar to the DIS case (0). Here a is the longitudinal momentum fraction of the 
quark, carried away by the photon. The LC wavefunctions for DY can be written in 
the same form as in DIS (|69|) , 



<4(«' P) = ^ (x.d^'U^) Koivp). (130) 



Summation over quark helicities and polarizations of the photon is understood in (|69D. 
Here, the extension parameter 

r]^ = mja^ + M^ {1 - a) . (131) 



is the anolog of the extension parameter e (|47|) in DIS. As usual, the index T stands 



for transverse photons and L for longitudinal. The x are again two component spinors. 
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but the operators O^lg read now |T^ 

d^.g = imja^e- {fix a) -i{2- a){e-Vp) +ae- (a xWp), (132) 

a^.^ = 2M(l-a). (133) 

— * 

The two dimensional gradient Vp acts on the transverse coordinate p, n is the unit 
vector parallel to the momentum of the projectile quark and a is the three vector of the 
Pauli spin-matrices. The LC wavefunctions for the transition q —>■ 7*g read explicitly 
for a given flavor of unit charge 

\%*,{a,p)\' = |^^.,(a,p)P + |vl/^.,(a,p)p, (134) 

l*^.>,P)r = ^{myKUrjp) + [l + il-af]v'Kl{rjp)}, (135) 

l*^*,(«,P)r = ^M'{l-afKlM. (136) 



vr^ 



Comparing (|135|) and ( p.36| ) with their DIS counterparts (^5[) and (^) shows that the 



factor Nc is no longer present in the q — > 7*g LC wavefunctions. This corresponds to 
the Nc in the denominator of ( p.23| ). Furthermore, we see that |\1/L (a, p)p has an extra 
factor of 2, because in the DY process, one has to sum over the transverse polarizations 
of the photon, rather than average like in DIS. 

For embedding the partonic cross section ( |129| ) into the hadronic environment, one 
has to note that the photon carries away the momentum fraction Xi = ( \/xp + At + 
xf)/2 from the projectile hadron. The hadronic cross section reads then 

da _ aem xi [^ ^Sr^ ^2 ( f^\ 4. pA) dajqp -» qYp) 

dM^dxp ~ 67rM2 xi + X2 X, a2 2^ -^ V^ \ a J ^^\ a J ) d\na 



a^ 



(137) 
-\- X2 Jxi Ci \ a J d\na 



67rM2 xi + 

where the factor aem/(67rM^) accounts for the decay of the photon into the lepton 
pair. Remarkably, the parton density of the projectile enters just in the combination 
F2 (P^), which is the structure function of the proton. Therefore we did not include the 



fractional quark charge Zf in the DY wavefunctions (135, 136 ). The structure function 
F2 is needed at large values of xbj- We will employ the parametrization from | |124[ in 
our calculations, unless mentioned differently. 

An interesting feature of the light-cone approach is the appearance of the dipole 
cross section in ( |129| ), although there is no physical gg-dipole in fig. |^. The explanation 
is illustrated in fig. |12[ The antiquark enters, when one takes the complex conjugate of 
the amplitude. The dipole cross section appears, because the quark is displaced in the 
impact parameter plane after radiation of the photon. If p is the transverse separation 
between the quark and the photon, the 7*g fluctuation has a center of gravity in the 
transverse plane which coincides with the impact parameter of the parent quark. The 
transverse separation between the photon and the center of gravity is (1 — a)p and 
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(a) 




(b) 



Figure 12: Two out of four con- 
tributions for the DY cross sec- 
tion. In the first row, the photon 
is radiated on both sides of the cut 
after the quark gluon vertex. The 
antiquark enters through the com- 
plex conjugate of the amphtude on 
the r.h.s. of the cut. The second 
row shows the interference term. 
The graphs in (a) and (b) have a 
relative phase factor exp(iap- p^), 
where px is the transverse momen- 
tum of the photon. This gives 
rise to the color screening factor 
in the dipole cross section (^) as 
explained in the text. 



the distance between the quark and the center of gravity is correspondingly ap. A 
displacement in coordinate space corresponds to a phase factor in momentum space. 
The two graphs (a) and (b) in fig. |1^ have the relative phase factor exp(iap-pT), which 
produces the color screening factor [1 — exp(iap ■ pr)] in the dipole cross section (^81). 
Therefore, the argument of the dipole cross section is ap. 

The physical interpretation of ( |129| ) is the same as for the DIS case, cf. section |2.3| . 
The projectile quark is expanded in the interaction eigenstates (p8|) . We keep only the 
first eigenstate. 



\q) 



Z2\qbare) +cL|g7*) + ... 



(139) 



where, Z2 is the wavefunction renormalization constant for fermions. Different eigen- 
states scatter with different amplitudes off the target. This disturbs the coherence 
between the eigencomponents and the photon is freed. Since there is no real dipole in 
DY, we will not include a nonperturbative interaction, as it was done for DIS in section 
but always use the perturbative LC wavefunctions 



135 



The light-cone approach to the DY process is applicable only at low X2, because the 
dipole cross section takes only sea quarks generated from gluon splitting into account. 
Therefore, the dipole cross section is proportional to the target gluon density. The 
statement, whether a sea quark belongs to the target or to the projectile, is however 
frame dependent. If the projectile quark or antiquark becomes slow in the limit a — > 1, 
in the infinite momentum frame, it can be interpreted in the infinite momentum frame 
as anti-seaquark or seaquark of the target which annihilates with the projectile parton. 
No annihilation with valence quarks from the target is taken into account. Note however 
that valence as well as sea parton distributions of the projectile are taken into account 
in the light-cone approach, if one employs a parametrization of the projectile structure 
function in ( |137| ). Therefore, the formulation of the DY process presented in this 
section is not fully symmetric between projectile and target. 
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It is also interesting to see that with the substitution 



1 + a a 

the transverse DY cross section (|137|) can be rewritten as 



« = ^ , P=-, (140) 






Jo ^ 



dM^dxF Gn^M^ X1+X2 Jn ^ f (1 + a 



X 



Jd'p{[l + 2a{l + a)]e'Kl{ep) + m}Kliep)}a,,{p}, (141) 

and similarly for the longitudinal case, where e = 5(1 + 5)M^ + ml is the analog of 
the DIS extension parameter (^Tf). One immediately recognizes the similarity with the 
qq wavefunction (|45|) in this equation. Indeed, the exact LC wavefunctions for DIS can 
be obtained by analytic continuation from timelike to spacelike region M^ -^ —Q^ and 

5 — > —a. 



In the rewritten form (|141|) , our expression for the DY cross section agrees with the 



result of Brodsky, Hebecker and Quack |T^ , where factorization in impact parameter 
space was considered. Note that the discussion of hard and soft contributions to DIS in 
section p.3| immediately carries over to DY. Configurations with 5 — > are the analog 
of the aligned jet configurations in DIS. Since 5 = p'g^/q^ is the energy of the quark in 
the final state divided by the photon energy, large distances become important, when 
5 — s> and a — > 1, respectively. It appears however only one limiting case in DY, 
whereas in DIS (^^ , two extremely asymmetric cases allow for large distances. 

For zero quark mass, the divergencies arising from the endpoints in the a- 
integration, correspond to the logarithmic divergencies occuring in the standard ap- 
proach, cf. sections p^ and |2.5|, which were absorbed into the nonperturbative parton 



densities. A similar procedure can be performed in impact parameter space |jTj]. In 
the parton model, the differential DY cross section can be written as 



J2 Z} {qf{x,)qf{x2) + g/(x2)g/(xi)} . (142) 



dxpdM^ 9M^ xi + X2 

The region 5 — »■ is the DY analog of Bjorkens aligned jet configurations and gives a 
leading twist contribution to the transverse DY cross section. As explained above, a soft 
outgoing quark in the target rest frame can be regarded as part of the nonperturbative 
sea-antiquark distribution of the target in the infinite momentum frame. Neglecting 
the quark mass m/, the divergent part of the integral ( p.41| ) can be recasted into the 



form (|142|) and one can read of the quark and antiquark distributions of the target at 



low X2. One find [l^ 



3 M^ r P 

X2qfJ{x2) = -^ y d^J d'paM^Kl{ep)ag,{p), (143) 

where A <^ 1 is a cutoff and the approximation e — > 5M^ for 5 — > was used. On the 
other hand, in the parton model of DIS, one has 



^TP 



^T '^ = 5^ ^J Ufi^Bj) + qfixBj)} , (144) 



4vr2a^,„ ^ 
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and one can again extract the divergent part of the integral (H3) for the transverse 



cross section and read of the target parton distributions [jl4 

XBjqffixBj) = ^ / da J d''paQ^Klisp)a,g{p). (145) 

Note, that in DIS a slow quark (a -^ 0) seen in the target rest frame, can be inter- 
preted as part of the sea-antiquark distribution of target in the infinite momentum 
frame. Obviously, the equations ( |143| ) and ( |145D are identical and therefore the DIS 



parton distribution is identical to the DY parton distribution. This manifestation of 



factorization in impact parameter space was first pointed out in |1J . 

The transverse momentum distribution of DY pairs can also be calculated in the 
light cone approach |jl5|. The differential cross section is given by the four- fold Fourier 
integral 



da{qp — > g7*p) 



d\n.a(PpT (27r)' 



/ d^pid^p2 exp[ipT ■ (Pl - P2)]^^.q(a,Pl)^7*g(«,P2) 



X - {(yqq{api) + (Tqq{ap2) " (T qq{a{pi - Pa))} • (146) 

after integrating this expression over pxi one obviously recovers ( |129| ). The negative 
term in the curly brackets arises from the interference between the two graphs for 
bremsstrahlung. The expressions for the LC wavefunctions needed here are 

*7^g(«,Pi)*7*g(«,P2) = ^ <^ mJo^Ko (r/pi) Kq (r/pa) 

+ [l + (l-«)^]r/2^i:^Ki(r/pOKi(r/p2)|, (147) 

P1P2 J 

vl>;^^(«,p-'i)vl>^.^(a,p-'2) = ^M2(l-a)^Ko(r/pi)Ko(r/p2). (148) 

The hadronic cross section is then given by 

da aem 1 f^ da^p /Xi\ dajqp -^ q^p) ^^^gs 



dM'^dxpd'^PT QttM'^ Xi + X2 Jxi ct \aJ dhiad^pT 

in analogy to ( |137] ). Three of the four integrations in ( p.46| ) can be performed analyti- 
cally for arbitrary a^q. This will be done in appendix |A.1| . 

2.7 The dipole cross section 

In this section we present a parametrization of the dipole cross section which is adjusted 
to DIS data and to total hadronic cross sections. In the preceding sections we have 
expressed the cross section for DIS and the DY-process at low x in terms of the dipole 
cross section. Thus the dipole cross section is the crucial missing input one needs to 
perform a calculation that can be compared to experimental data. The dipole cross 
section is of interest not only for DIS and DY. Total hadronic cross sections, like for 
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pion-proton scattering can also be written in terms of the dipole cross section. In this 
case, the pion wave-function in coordinate representation \l/7r(p), gives the probabihty 
to find a quark- ant iquark pair with a given transverse separation inside the pion. If 
only the lowest Fock state is taken into account, the total cross section at cm. energy 
s can then be written as 



^TTp\ 



ld'p\^^ip)W,is,p), (150) 



in complete analogy to the DIS case. This is remarkable since hadronic cross sections 
are governed by nonperturbative long distance physics in contrast to DIS. Still, (|150|) is 
valid, because the decomposition into interaction eigenstates (|58|) is completely general 
and does not rely on perturbation theory. 

The dipole cross section is largely unknown, only at small distances p it can be 
expressed in terms of the gluon density (0). However, several parametrizations exist 
in the literature, describing the whole function aqq{p), without explicitly taking into 
account the QCD evolution of the gluon density. A very economical parametrization 
is provided by the saturation model of Golec-Biernat and Wiisthoff []7U|, P0| , 



(Tqq{x, p) = (To 



, p'Ql 

1 — exp ' 



A{x/xoy 



(151) 



where Qo = 1 GeV and the three fitted parameters are ctq = 23.03 mb, Xq = 0.0003, 
and A = 0.288. This dipole cross section vanishes oc p^ at small distances, as implied 
by color transparency and levels off exponetially at large separations, which reminds 
of eikonalization. The authors are able to fit all available HERA data with a quite low 
X^ O and can furthermore also describe diffractive HERA data |^^. It is argued that 
saturation is necessary to obtain a nearly Q^ independent ratio of the diffractive to the 
total cross section, which cannot be reproduced with a dipole cross section that rises 
quadratically up to infinity. 

There is however no xbj in hadron-hadron scattering and the appropriate variable 
to describe the energy dependence of the dipole cross section is s. This is different in 
DIS where the virtual photon scatters off pointlike constituents inside the proton and 
Xbj is the correct variable to use. Parametrizing the dipole cross section as function 
of s will automatically lead to a violation of Bjorken scaling. This is easy to see, since 
s f» Q'^/xbj- If the cross section is supposed to grow like a power of s, it will grow 
with approximately the same power in Q^. 

A two pomeron model of the dipole cross section as function of s, motivated from 



Regge phenomenology, was given in |^T|. It was argued in |82| that the authors of |81 
had to parametrize the dipole cross section as a function of p^s, to get approximate 
Bjorken scaling. The ansatz of [^ can describe deep inelastic data up to Q^ = 60 
GeV^, with 10 parameters. 

However, as can be seen from fig. |^, Bjorken scaling is quite strongly violated at 
low Xbj- Since we have in mind to describe hadron-hadron scattering and DIS with 
the same dipole cross section, we write the dipole cross section as function of s. We 
are aware of the fact that our ansatz will fail at very large values of Q^. Nevertheless, 
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practically all physics which is interesting to us happens at quite small values oiQ^ < 10 
GeV2. _ _ 

An ansatz inspired by the saturation model of [[ 
reads 



80|| , was presented in [|12|. It 



0", 



QQ 



(s,p) =ao{s) 



exp 



P' 



with ro (s) = 0.88(s/so)-°-i^ fm, sq = 1000 GeV^ and 



cro(s) 



a 



tot\ 



H 



Vch/n 



(152) 



(153) 



where cr^^(s) = 23.6(s/so)'^'°^ mb and {rl/^)^ = 0.44 fm is the pion charge radius. This 
cross section is proportional to p^ for p ^ 0, but flattens off at large p. The energy 
dependence correlates with p. At small p the dipole cross section rises with a hard 
pomeron intercept, 0.36, and at large separations it still depends with a soft pomeron 
intercept, 0.08, on energy. This is different in ( |151| ), where the maximum value of 
the dipole cross section is independent of energy. We also mention that a value of 
(To = 23.03 mb is too small to reproduce the pion-proton cross section. Note that 
the energy behavior of our parametrization is also conceptually different from Regge 
phenomenology, where the pomeron intercept has to be independent of the hadron. We 
think however, that the growth of the cross section with s should depend on the typical 
distances involved and our ansatz smoothly interpolates between cross sections in hard 
QCD processes and soft hadronic scattering. The saturation scale ro(s) decreases with 
increasing s. This means that at very high s, only very small dipoles exhibit a hard 
pomeron behavior. At separations p ^ rQ{s), the dipole cross section grows only slowly 
with energy. The ansatz ( |152|) ensures that the dipole cross section is a monotonous 
function of s and p. This is not the case for the parametrization of |^, which can 
decrease with increasing p. 

The energy dependent saturation value of the parametrization ( |152| ) is chosen in 



such a way, that the pomeron part of the fit by Donnachie and Landshoff 
automatically reproduced, if one employs a gaussian pion wavefunction, 

3 / 3p2 



125 



I^ApW 



exp 



(154) 



To understand the factor of 3/8, note that it has to hold 

because p is only a two dimensional vector (2/3) and it has the meaning of a diameter, 
rather than a radius (2^). 



d'pp'\^Ap)\' 



(155) 



The dipole cross section (|152|) is however adjusted only by eye and no x^ ^^is 



calculated. Furthermore, it tends to overestimate F2 at Q^ > 10 GeV . We improve 
the parametrization (|152| ) in several ways. First, we use the more recent fit by Cudell 
et al. 11261 , 

-0.357 / Q \ -0.56 / Q \ 0.0933 



a^ 



26.2 mb 



1 GeV^ 



-7.63 mb 



1 GeV^ 



12.08 mb 



1 GeV^ 



(156) 
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p- / fm^ 



Figure 13: The dipole cross section ( |157| ) for -y/i = 200 GeV (HERA), solid line, and 
for ^/s = 1.4 TeV (THERA), dashed line. The left figure shows the dipole cross section 
vs. p, while the figure on the right displays aqq{s,p) vs. p^. At large separations, the 
dipole cross section becomes constant in p, but still grows slowly with energy. At small 
separations the dipole cross section depends quadratically on p. 



for ^/s > 9 GeV^. which has a slightly larger pomeron intercept and agrees better with 
experimental data than the old version by Donnachie and Landshoff ||125|| . A pomeron 
intercept larger than 0.08 is also suggested by recent SELEX measurements [|1271| , fig. 



21 



To improve the agreement with data, we multiplied (|152|) by a function that leads 
to a suppression at small value of p and becomes 1 at large p. We have tried several 
functions. Best agreement with HERA data is obtained with the parametrization 



(^qq{s, P) 



o-o(-s) 



1 — exp 



P' 



rl{s) 



where 



ao{s) = 12.08 mb 



0.0933 



1-0.9 exp - 



3rg 



P' 



[s] = (2.19 fm) 



1 GeV2 

0.44 fm^ 



Vchl-K 



n 



1 GeV2 
0.27 fm. 



-0.34 



(157) 

(158) 
(159) 

(160) 
(161) 



Our result is depicted in fig. |T3[ The energy dependence of ( [L57| ) is qualitatively 
the same as discussed for (|152|) , because ri does not depend on s. In the left figure, 
we plotted the dipole cross section vs. p. One can see how the cross section increases 
from HERA to THERA energies. The figure on the right shows aqq{s,p) vs. p^. At 
small separations, the dipole cross section increases quadratically with p. 

We determined the 3 parameters, ri and the two numbers in (|160|) , using Hl- 
data up to Q^ = 65 GeV^. We have tried to achieve good agreement with the same 
parameters for perturbative and nonperturbative wavefunctions. A comparison with 



data is presented in figs. p^4Hr7| . One observes that with perturbative LC wavefunctions. 
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(^,^), the very low Q^ data is not well reproduced. The agreement can be improved by 
introducing an unphysically small quark mass. This would however allow for arbitrarily 
large separations of the pair, at variance with confinement. We therefore fix the mass 
for u- and rf-quarks at a typical inverse confinement radius, 

mu,d = 200 MeV. (162) 

We also include strange and charm quarks in our calculation, 

nis = 350 MeV , m, = 1.5 GeV. (163) 

With these parameters, the lowest x^ that can be obtained with perturbative wave- 
functions is 

Xpi(Hl, g^ < 65 GeV^) = 523/233 d.o.f. (164) 

Statistical and systematic errors are added in quadrature. However, the dashed curves 
still make a reasonable impression at not too high Q^, when one looks at figs. |T^|T7. 



We will therefore perform all our calculations with perturbative wavef unctions, too. 
Better agreement is obtained with the nonperturbative wavef unctions, ( |77| , |78D , which 
describe the data especially well at low Q^. We obtain 

x'pi(Hl,g2< lOGeV^) = 86.4/100 d.o.f, (165) 

X^pi(Hl,g2<30GeV2) = 214/182 d.o.f. (166) 

For the nonperturbative wavef unctions, we fix the quark masses at rUu = mii = 10 MeV, 
nis = 150 MeV and rric = 1.5 GeV. We have also plotted the result of our calculation 
for very high virtualities Q^ = 2000 GeV^, see figs. ^^-[16|. The parametrization (|157|) 
overshoots the data at large Q^. 
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Figure 14: The points show HI data | 129 |-| 132 | for the proton structure function 
F2 vs. s. Only statistical errors are shown. The solid curves are calculated with 
the parametrization ( p.57| ) of the dipole cross section and with the nonperturbative 



LC wavefunctions (|77| , [78|) . The dashed curves are the same but calculated with the 
perturbative LC wavefunctions (^,^). The parameters in (|157|) are adjusted to this 
set of data. See text for more details. 
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Figure 15: The same as fig. |Tj. The points shown here are not used to determine the 
parameters in 
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Figure 16: The same as fig. |l^ but the points are from ZEUS ||60|, |6l], |133|| - |[T36[| . These 
points are not used to determine the parameters in ( |157| ). 
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Figure 17: The same as fig. |l^ but the points are from ZEUS ||60|, |6l], |133|| - |[T36[| . These 
points are not used to determine the parameters in ( |157| ). 
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Figure 18: The LC wavefunctions integrated over a, (|167D , at ^/s = 200 GeV. The 
sohd curves are calculated with the nonperturbative wavefunctions, (|77| , |T8D , while the 
dashed curves are calculated with the perturbative ones, (^ , ^6[) . The lower two curves 
(at large p) are for Q^ = 10 GeV^ and the upper ones for Q^ = 0.5 GeV^. 



To see in more detail, how the nonperturbative interaction modifies the LC wave- 
functions and which distances contribute to F2, we calculate the following three quanti- 
ties. First, we integrate the LC wavefunctions over the longitudinal momentum fraction 
a and over the angle, 

27rp [^ , |,T,T,L/ 



Ao{p) 



ap 



da\^qk (a,P)l' 



(167) 



The result is shown in fig. |T8[ The results for our calculation of the wavefunctions 



integrated over a and angle weighted with the dipole cross section, 

-4i(p) = agg{s,p) rfa|^Ji^(a,p)p, 



ap 



:i68) 



can be seen in fig. |T9|. This quantity is important for the structure function F2. Finally, 
we weight the integral with agg{s,p) 



A2{p) 



2lTp 9 

— -(J 



<?<? 



[s,p) 



C^«l*r^^(«>P)P, 



(169) 



which is important for diffraction. The result is depicted in fig. ^. In all three cases, 
the factor 27rp arises from the area element. It has to be included, if one wants to get 
a correct impression of the distances contributing to F2. For convenience, we divide 
AQ-A2 by the fine structure constant Ogm- AH calculations are done at ^/s = 200 GeV^ 
and for two values of Q'^, namely for Q^ = 0.5 GeV^ and for Q'^ = 10 GeV^. 

It is hard to see in fig. |T^, how the nonperturbative interaction modifies the trans- 
verse LC wavefunctions. Even at Q^ = 0.5 GeV^, the curves are almost identical. The 
transverse part diverges oc 1/p at p ^ 0. The wavefunction is not normalizable. The 
curves rapidly decay as p becomes larger, cf. (|62D. The situation is different in the lon- 
gitudinal case. Here, Aq goes to at p ^ 0. The longitudinal wavefunction is square 
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Figure 19: The LC wavefunctions integrated over a and weighted with <Jqq{s, p), (|168D 
at ^/s = 200 GeV. See fig. for more explanations. 



integrable. At high Q^ the integral Aq has a peak at small p ~ 1/Q and then decreases 
exponentially. At low virtuality, this peak is much less pronounced and Aq has quite 
a long tail into the large distance region. Somewhat surprisingly, the longitudinal LC 
wavefunctions are stronger affected by the nonperturbative interaction at low Q^ than 
the transverse. At Q^ = 10 GeV^, no effect of the nonperturbative interaction is visible 
any more. 

The quantity Ai displayed in fig. |19| integrated over p is nothing but the structure 
function F2, up to a factor Q^/(47r^). The region of small p is suppressed compared to 
fig. 18, because of an extra factor p^ from the dipole cross section. Therefore, larger 
distances as one would expect contribute. Both, for low and high Q^, the curves have 
a long tail into the large distance domain. Again, there is practically no infiuence 
of the nonperturbative interaction at high virtuality, but at low Q^ and large p, the 
perturbative part is slightly suppressed compared to the nonperturbative. This means 
that our choice of the quark mass, squeezes the wavepacket somewhat more than the 
interaction. This can also be seen in figs. p!^p!7| , where the perturbatively calculated 
curve is always below the nonperturbative. At very high Q^, the curves become in- 
distinguishable. We see again, that the longitudinal part is more strongly infiuenced 
by the interaction, than the transverse. Note also, that the longitudinal part is much 
smaller than the transverse, namely of order 20%. 

We also plot the «-integrated LC wavefunction, weighted with o"gg(s,p). In fig. [2^, 
we see that the small distance region is suppressed by an additional factor p^, compared 
. Small distances play numerically no role, if one would integrate A2 over 

This illustrates our discussion in section |2]J, where we argued that diffraction is 



to fig. 

P 

dominated by soft physics, despite the large virtuality. 

While the original parametrization of the dipole cross section from [|12| reproduced 
identically the pomeron part of the Donnachie-Landshoff fit, the improved parametriza- 
tion (|157|) does not have this property. We therefore check, how well we can reproduce 
the total iT~p cross section with the wavefunction ( |154|) . Of course we have to add the 
reggeon exchange contribution from (|156| ), since the dipole cross section describes only 



2 SOFT CONTRIBUTIONS TO HARD QCD REACTIONS 



51 



nonperturbative 
perturbative 




250 

S 200 

'> 
c» 
O 

:: 150 

c 

■5 

•I 100 

c 

_o 

< 50 



nonperturbative 
perturbative 




p/fm 



p/fm 



Figure 20: The LC wavefunctions integrated over a and weighted with o"g^(s, p), (|169D 
at i/i = 200 GeV. See fig. for more explanations. 



the pomeron part. The result is shown in fig. ^ We see that the cross section is not 
identically reproduced, as could have been expected. Our result is slightly smaller than 
the fit ( |156D , because we multiplied the old parametrization ( |152| ) with a function that 
is always smaller than 1. Nevertheless, the pion-proton cross section is reasonably well 
reproduced. The cross section calculated from our parametrization ( [L57| ) grows with 
the same power of s at high energy. One gets the impression from fig. ^ that the data 
grow even faster. At present, there is however no data available at higher energies than 
v^ = 35. 

We can now proceed and investigate how well our parametrization describes DY 
data. The hght-cone approach is applicable only at low X2, because it neglects valence 
quarks in the target. At present, there are however not many data for DY cross sections 
at low X2 = i^/xp + At — xp)/2 or at large xp, respectively. Only the data at the two 
largest values of Xp from the E772 collaboration, correspond to values of 0.01 < X2 < 
0.1. The result of our calculation, using (|137|) , is shown in fig. |2[ Although the data 
scatter a lot, one observes that the M dependence is only approximately reproduced. 
We believe that this can be improved by modifying the dipole cross section. Remarkable 
is also that the absolute magnitude of the cross section is quite well reproduced, without 
introducing a i^-factor. Indeed, we believe that it is not legitimate to use a K factor 
in our approach, since the dipole cross section is supposed to parametrize all higher 
order corrections. 

Furthermore, we also calculate the transverse momentum distribution of DY dilep- 
ton pairs from ( |149| ). The result is depicted in fig. |23|. The DY cross section is finite 
at pt = 0, in contrast to the result of first order pQCD. We do not compare to data 
in fig. ^, because all available data are integrated over xp and are therefore contami- 
nated by valence quark contributions. An extraction of the low-x part of the transverse 
momentum distribution is in progress [ |138|| . Note however, that the differential cross 



section is usually parametrized as |139, 115 1 



da 
dpp 



oc 



'i+{PT/Poyr 



(170) 
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Figure 21: The points show the world data ||128|| for the 7i~p total cross section. The 
point at ^/s ~ 35 GeV is from the SELEX collaboration ||127|| . Only statistical errors 
are shown. The solid curve is calculated with the dipole cross section ( |157| ), while the 
dashed curve is the fit ( |156| ) from [p.26|| . 



while our result decays much weaker and behaves approximately like 



1 



oc 



[1 + (pt/I GeV) 



212' 



(171) 



The behavior oc p^ is expected from pQCD at large p"^ > M^ and the finiteness at 
Pt ^ looks encouraging. It is at present not clear to us, whether a modified dipole 
cross section could reproduce the familiar shape ( |170| ), or whether the resummation of 
soft gluon radiation from the projectile quark requires additional modifications of the 
approach. 

To summarize, our parametrization reproduces well the proton structure function 
F2 from the lowest values of Q^ ^ 0.1 GeV^ up to Q"^ ^ 30 GeV^, although it does 
not respect Bjorken scaling. The low Q^ data is much better described if one employs 
the nonperturbative wavefunctions. Also at Q^ > 30 GeV^, the data is quite well 
reproduced if one looks at the plot only by eye. Since we have parametrized the dipole 
cross section as function of s, rather than xbj, we can also calculate total meson-proton 
cross sections, given an appropriate meson wavefunction. The total 7i~p cross section 
is reasonably well reproduced and also the few available DY data can be described 
without invoking a ii'-factor. Only the DY transverse momentum distribution looks 
different from what one would expect. At present, all data on the transverse momentum 
distribution is however integrated over xp, making a comparison with calculations 
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Figure 22: The points represent the measured DY cross section from ||137|| . Only 
statistical errors are shown. The curves are calculated with the dipole cross section 
( |157|) without any further fitting procedure. In order to make X2 small, xp has to be 
large. In these two figures 0.01 < X2 < 0.1. 



impossible. The differential cross section could give much more detailed information 
about the dipole cross section than the presently availabe total cross sections. Certainly 
( p.57| ) is not the final answer and has to be improved further and compared to more data. 
Especially diffractive DIS data which are more sensitive to the large distance region 
have not been taken into account so far. We will however compare our predictions for 
shadowing with data and nuclear shadowing is closely related to diffraction, as will be 
explained in the following section. 
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Figure 23: The transverse momentum distribution for DY pairs calculated from (|149|) 
with the dipole cross section ( |157D . 
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3 Nuclear shadowing in DIS and for DY 

The typical energy and momentum scales in nuclear physics are much lower than the 
virtualities reaches in DIS. One could therefore expect that using nuclei instead of 
protons as target would only increase the cross section. There are however very pro- 
nounced nuclear effects in DIS. At intermediate values oi xbj one observes a suppression 
of the structure function ratio, fig. ^, known as the EMC effect [|im . This effect was 
unexpectedly discovered at CERN. It is described theoretically in terms of new degrees 
of freedom in nucleus, i.e. others than nucleons, see review ||140|| . We do not discuss the 
EMC effect any further and turn our attention to lower values of xbj- After a small 
enhancement in the region of xbj ~ 0.1, which is called antishadowing, one reaches 
the shadowing region, where cr"^*"^ < Aa"' ^ {A is the nuclear mass number), i.e. a 
considerable depletion of the nuclear structure function is observed. During the last 
decade, nuclear shadowing has been measured at CERN (EMC, NMC), at Fermilab 
(E665) and recently also at DESY (HERMES). Values oixBj < 0.001 can be reached, 
but most of the data are at higher values. 



'•:'(-^' y J 




Figure 24: Schematic plot of nuclear effects in the ratio of the structure functions 
F^ of a nucleus and the proton structure function F^ vs. x = xbj- The suppression. 

The 



^2 



/AF2 < 1, at Xbj < 0.1 is called nuclear shadowing. See text for more details, 
figure is taken from | ]141 |. 



After the pioneering work of Gribov 0], where shadowing arises from the interaction 
of virtual photons at the surface of the nucleus, this effect was not unexpected. Nuclear 
shadowing is of high theoretical interest, because it can be understood in terms of the 
fundamental degrees of freedom of QCD, quarks and gluons. In the preceding sections 
the ground for an investigation of nuclear shadowing in terms of the color dipole picture 
was prepared. In this part of the work, an improved version of the Glauber-Gribov 
theory is developed, which includes the nuclear form factor in the usual Glauber series. 
The interplay of multiple scattering and coherence are discussed in detail. 
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3.1 Shadowing and Diffraction 

The partonic interpretation of shadowing depends on the reference frame. In the 



infinite momentum frame of the target, shadowing is due to parton fusion ||142|| - ||145|| , 
which leads to a reduction of the parton density at low xbj- This can be understood 



intuitively in the following way |[146|| : In DIS a picture of the microscopic substructure 
of the target is taken. However, this cannot be done in an infinitely short time. The 
uncertainty relation tells us that the time needed increases as xbj becomes small. This 
can be seen from the energy denominator in (|43|) on page |l^. More precisely, the role 



of time is played by the logarithm of xbj, as can be seen from the double leading log 
DGLAP equation, 

d^XBjGixBj^Q^ _ N,a 2 

d\n{l/xB,)d\nQ^ - -r'^B.GixB.^Q )• (172) 

For fixed strong coupling constant a^, the solution behaves asymptotically like (pQl). 
While a fast moving nucleus is contracted to a pancake shape, the cloud of partons with 
very small momentum fraction is much less contracted and parton clouds from different 
nucleons will overlap. As the density becomes higher and higher with decreasing xbj, 
the partons start to feel the neighborhood of other partons and fusion processes like 
GG —>■ G and GG — > qq will occur. These fusion processes reduce the parton density 
in a nucleus compared to a free nucleon. 

A very intuitive picture arises in the rest frame of the nucleus, where the same 
phenomenon looks like nuclear shadowing of hadronic fluctuations of the virtual photon 



||TT| , |70| , |7T|, ||147|| - ||153|| . Shadowing occurs, because of multiple scattering of the 
hadronic fluctuation inside the target. The nucleons at the surface of the nucleus 
cast a shadow on the inner nucleons, which then have less chances to interact with 
the photon. Multiple scattering at high energies is commonly described in the optical 
model of Glauber-Gribov theory. The basics of this approach have been formulated a 
long time ago by Glauber |^ within nonrelativistic quantum mechanics. The optical 
model is usually derived in eikonal approximation with the assumption that the phases 
obtained in different scatterings are additive. Of course, one would like to have a 
connection to quantum field theory and the Feynman diagram technique, which allows 
to include relativistic kinematics. This has been elaborated by Gribov [0]. However, 
the relation between Glauber theory and Feynman diagrams has been fully explored 
only for the double scattering term. The theory has been developed by many authors 



since then. In [ p.54|| , Gribovs approach is generalized to non forward scattering and 



an expression is obtained which also includes inelastic transitions. For a discussion of 



Glauber-Gribov theory in the context of Regge phenomenology see for instance |155 



We give a short summary of the Glauber-Gribov theory, following ||154||. For sim 



plicity, only forward scattering is considered, although the approach can be generalized 
to include also small momentum transfers. All correlations between the nucleons will 
be neglected in the following. The total amplitude for scattering off a nucleus is written 
as a sum over ra-fold scattering amplitudes. 

A 

F(.) = ^F(")(.), (173) 

n=l 
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where s is the photon-nucleon cm. energy squared and F^"'\s) is the amphtude for 
an n-fold scattering. It can be written as an integral over impact parameter (b) space 
and over the longitudinal coordinates Zi of the scatterings, 

oo oo oo 

F^^'^s) = l^j d% dZiPA{b,Zi) dZ2PA{b,Z2)... dZnPA{b,Zn) 

X J2 /r-/.ie-'^^-*-''i(^^-^^)/h,^h,e-^''-^-'^2(^^-^^) . . . A„_,^,*. (174) 
{h} 



Here, i^ Ri 2s is the flux factor. The expression ( |174| ) is illustrated in fig. ^. The 



nuclear density pa, which is normalized to A, enters after integration over the wavefunc- 
tion of the nucleus, neglecting correlations. It is assumed that the impact parameter b 
does not change in a collision. The limits of the integrations over the longitudinal coor- 
dinates Zi originate from the condition that the (z + l)th collision has to take place after 
the ith one. After the first collision, the virtual photon converts into a multiparticle 
hadronic state hi. The amplitude for this process is denoted by f^*^hi- Furthermore, 
the longitudinal momentum transfer qL^*^hi appears as an oscillating phase factor for 
each intermediate state. One has 

g2 + Mi 

qir^h, = — — -, (175) 

where M/j^ is the invariant mass of the hadronic state hi. Finally one has to sum 
over all states hi that can be reached by multiple scattering. According to the optical 
theorem, the total cross section is then given by 

a^*\s) = ^. (176) 

Straightforward calculations with ( |1 74| ) are however not possible, because the in- 
termediate states which can occur in an n-fold scattering are not known. Therefore, 
one is forced to introduce additional approximations. The idea to represent the virtual 
photon as superposition of hadrons with the same quantum numbers, namely vector 
mesons, leads to (generalized) vector dominance model (G)VDM |0, |102|| . The advan- 



tage of this approach is that the masses of vector mesons are well defined, in contrast 
to the mass of a qq dipole in the mixed representation. The masses of the intermediate 
states have to be known, since they enter the longitudinal momentum transfers qi- On 
the other hand it is not known, how to calculate scattering amplitudes, higher than the 
double scattering term, because vector mesons are no eigenstates of the interaction. A 
possible prescription would be to omit all off diagonal transitions which is called the 
diagonal approximation. Note that in the VDM, which is the diagonal approximation 
of the GVDM, one has to choose a vector meson-nucleon cross section that behaves 
as avp oc 1/My, where My is the mass of the vector meson, in order to reproduce 
the Bjorken scaling of the proton structure function. With a nearly constant avp, F2 
would increase linearly with Q^. However, with such a fast decreasing meson-nucleon 
cross section, it turns out that nuclear shadowing is a higher twist effect. It was indeed 
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Figure 25: Illustration of the n-fold scattering amplitude F^"'\s) ( |174|) with intermedi- 
ate hadronic states 1 ... n — 1. The upper part of the figure corresponds to the photon, 
which converts in different hadronic states due to multiple rescattering. The lower part 
corresponds the nucleus, where each line represents a nucleon. The amplitude for the 
transition i —>■ j is denoted by fj in this figure. 



believed for quite a long time that shadowing vanishes at large Q^, until it was shown 
experimentally that this is not the case [HI]. This phenomenon is sometimes called the 
Gribov paradox ||156|| . The solution is that one also has to take into account off-diagonal 
transitions, which correspond to diffraction dissociation. These off-diagonal transitions 
interfere destructively with the diagonal ones and one can reproduce the Bjorken scal- 
ing with an My independent avp ||157|| . The diagonal VDM is applicable only at low 



Q^ < 1 GeV^, where the proton structure function indeed rises approximately linearly, 
see fig. p. 

The close connection between shadowing and diffraction becomes most transparent 
in the formula derived by Karmanov and Kondratyuk ||158|| . In the double scattering 
approximation, the shadowing correction can be related to the diffraction dissociation 
spectrum, integrated over the mass. 



a 



.7M 



Aa'^'^-Sn / (fb / dzipA{b, Zi) / dz2PA{b,Z2 



21 



dM^ 



da{YN -^ XN) 



X 



dM'^dt 



cos — 1 



t^o 



J J dM^dt 



Q' + Ml 

2u ^'' - 


-z,)] (177) 


Flilc.b). 


(178) 



t->0 



Note that the kinematical limit for t is slightly larger than zero. This result is easily 
obtained from ( |174| ), when one neglects the small real part of the scattering amplitude 
F{s). Here 



FJ(/„ b) 



dzpA^b, z)e 



iz/lc 



(179) 
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Figure 26: Illustration of the Karmanov-Kondratyuk formula ( |177| ). The left figure 
shows the single scattering contribution to DIS off nuclei. The double scattering term, 
shown on the right, gives the driving contribution to nuclear shadowing. It scales like 

Al. 



is the formfactor of the nucleus. Again the coherence length 

2iy 



L 



Q^ + Ml 



;i80) 



emerges. This length is of course already present in ( |1 74[ ) . The Karmanov-Kondratyuk 
equation is illustrated in fig. ^. The first term in ( |177|) is the single scattering term, 
which is proportional to A. The second term represents the double scattering con- 
tribution and is responsible for shadowing. To model the influence of higher order 
rescattering terms, the shadowing correction is often multiplied by the term 



exp 



o. 



eff 



dzpA{b,z] 



:i8ii 



with the effective absorption cross section cTe// ~ 20 mb for quark shadowing. Note 
that at very high energies, the coherence length becomes infinite and the oscillating 
exponential in the nuclear form factor is practically unity which corresponds to maxi- 
mum shadowing. Of course, the shadowing term still grows with energy, because of the 
diffractive cross section. As the energy becomes lower, the exponential oscillates more 
and more rapidly and the shadowing correction vanishes. This observation makes the 
importance of Ic very clear. 

A simple estimate that demonstrates the relation between shadowing and diffraction 
is performed in | ]159| ] , see also [|T^] . In |159| , an infinite coherence length is assumed 



and the diffractive cross section is parametrized in the way. 



dM^ 



da{YN -^ XN) 



dM'^dt 



Ba, 



DD 

7*p- 



;i82) 



t->o 



The diffractive cross section o"£p ~ O.la^*^ makes about 10% of the total cross section 
and the slope parameter which describes the t-dependence of diffraction is adjusted to 
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B ^ 8 GeV~^. Although no higher order scattering terms are taken into account, the 



authors of ||159|| obtain good agreement with experimental data at very low xbj- This 



analysis confirms the physical picture that shadowing in DIS is governed by coherent 



interaction of diffractively produced states in the nucleus, fig. ^. 

From this discussion, two fundamental requirements emerge, which have to be ful- 
filled for shadowing. 

• The mean free path of the hadronic fluctuation is long enough to have multiple 
scatterings. More quantitatively it must hold 

^ < Ra, (183) 



PACTeff 

where Ra is the nuclear radius. 

• The coherence length l^ has to be bigger than the mean inter-nucleon distance, 
lc> 2 fm. This condition ensures that the virtual photon coherently scatters off 
the nucleons. 

One can also estimate, which masses Mx can contribute to shadowing. Note that 

since 

da{YN -^ XN) 



dWj^dt 



« TTT (184) 



t^o M' 



for large masses, one formally has to integrate the spectrum up to the kinematical limit. 
However, due to coherence length effects, large masses cannot contribute to shadowing. 
This can be seen in the following way. The coherence length can be written as 

Ic = -^- = -^— (185) 

rriNXBj niNX-p 

with 

" = (?&■ ''''* 

One finds from the second condition that only masses 

M^<0.1s (187) 

are relevant for shadowing. Note that this quantity can become much larger than Q^ 
at low Xbj- It is however small compared to s. 

In the VDM and for the Karmanov-Kondratyuk formula, the multiple scattering 
amplitude ( |174| ) is evaluated in a hadronic basis. The disadvantage of using this basis 
is that one cannot calculate higher order scattering terms in an unambiguous way. 
We can however make use of our knowledge of the interaction eigenstates, namely 
partonic configurations frozen in impact parameter space. No off-diagonal transitions 
occur in this basis. One encounters however another severe problem, the mass of the 
eigenstates, which one has to know for the coherence length, is undefined in the mixed 
representation and one cannot give an explicit expression for the scattering amplitude. 
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Only in the limit /c — ^ oo, it is possible to resum the whole multiple scattering series 



|154|| in an eikonal-formula 

7*A _ o /"^2> f , /"^2^|,T,/ , X _|2 A _^ ( ^qg-(g, p) 



a' 



2 I Shi da I Sp\^{a,p\^\^ Tl-exp (-'^d^T{h) 

+ 2 j d% j da j — - j d'^pid'^p2\^qqG{a,aG,Pi,P2)f 

x(l-exp(-^^^^^i^l^^r(6) 

Summation over photon polarizations is understood. The nuclear thickness function 
T{b) = f_^dz pA{b,z) is the integral of nuclear density over longitudinal coordinate 
z and depends on the impact parameter b. The light-cone wavef unctions are given by 
(^) and (^) on page 17 for the free case and by ( [77D and ( |75D on page ^ including 



the nonperturbative interaction. The LC wavefunction for the gluonic component is 
given by ( |115| ) on page ^. The cross section for the three body system qqG ( |108D can 



be expressed in terms of the dipole cross section. In the following more quantitative 
discussion, the ggG-term is dropped and for the gg-part the shorthand notation 

a^*^ = (2 I d% Tl -exp L^d^T{b)\\\ (189) 

is introduced. The condition l^ ^ Ra insures that the r does not vary during prop- 
agation through the nucleus (Lorentz time dilation) and one can apply the eikonal 
approximation. Note that the averaging of the whole exponential in ( |189| ) makes this 
expression different from the Glauber eikonal approximation where (Tqq{s, p) is averaged 
in the exponent, 



d26(^l-exp(^-^T(6)jj. (190) 

The difference is known as Gribov's inelastic corrections 0. In the case of DIS the 
Glauber approximation does not make sense and the whole cross section is due to the 
inelastic shadowing. Indeed, a''*^ is at most of order 100 /ib, for real photons, which 
corresponds to a mean free path in nuclear matter of // ~ 600 fm. This means that 
the bare photon does practically not interact and therefore no multiple scattering can 
occur. The photon interacts only via its hadronic components. 

The condition Ic ^ Ra is however not fulfilled in most of the experiments. Only 
for the lowest values of xbj in the E665 experiment, this condition is met. For the case 
Ic ~ Ra, one has to take into account the nuclear form factor, i.e. the variation of p 
during the propagation of the qq fluctuation through the nucleus. This can be done for 
the gg-component of the photon and in double scattering approximation ||152|, |153||, 



a^*^ ^ Aa^*^ - \ {a\s,p))ld%Fl{l,,b). (191) 

It is not well defined what one has to take as argument Ic of the formfactor ( |179| ). The 
usual prescription is to set M^^ = Q^ which leads to Ic = l/{2mMXBj), cf. discussion on 
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page |r^. In the limit Ic ^ oo, ( |191D can be obtained from expanding the exponential 
in ([T89D . 

Like ( p.77|) , (|191|) also resembles the close connection between shadowing and diffrac- 
tion. The connection between these two equations is given by ( pT]) on page ^ In 



section |2.4| it is discussed in detail that only large size pairs give a leading twist con- 
tribution to the diffractive cross section. Similarly, also shadowing is a soft reaction, 
although Q^ is large. The leading twist contribution to shadowing for longitudinal 
photons comes from the ggG-Fockstate. The qq contribution vanishes with a higher 
power of Q^. The fact that only large size configurations contribute to shadowing is 
the solution of the Gribov paradox in the dipole language. It follows from the first con- 
dition for shadowing that the pair has to scatter at least with a cross section of 10 mb 
to be shadowed in a large nucleus. This quantity can be several times larger for a light 
nucleus, like carbon. In the VDM all hadronic fluctuations of the photon interact with 
a large cross section, which leads to a structure function F2 that is proportional to Q^. 
However, VDM cannot be applied at large Q^, since then most of the qq fluctuations 
interact with small cross sections and this produces the scaling of F2. Only the very 
unlikely Bjorken aligned jet configurations make shadowing a leading twist effect. 
To summarize, two problems remain under discussion 

• How can the nuclear formfactor can be included in the higher order scattering 
terms which are of great importance for heavy nuclei? For instance, the shadow- 



ing term in ( |189| ) for lead is half as big as Aa'^ p, so the need of the higher order 



terms is obvious. 

Even for the double scattering term it is still unclear which argument should enter 
the formfactor. Indeed, the effective mass of the qq fluctuation needed for the 
coherence length in ( p.8CI| ) cannot be defined in a representation with a definite 
qq separation. On the other hand, ([ITT]) exhibits an explicit dependence on Mx 



and the longitudinal momentum transfer is known. However, unknown in this 
case is the effective cross section o"e//- 

A solution to these problems is proposed in the next subsection. 

3.2 Light-cone approach to nuclear shadowing 

The goal of this subsection is to give a physical explanation of the light-cone approach 
to nuclear shadowing which includes the nuclear form factor in all multiple scattering 
terms. A more formal derivation is given in section ^.3|. In order to study the difference 



between the correct quantum mechanical treatment of nuclear shadowing and known 
approximations, we restrict ourselves to the qq Fock component of the photon, neglect- 
ing gluons. The nuclear antishadowing effect is omitted as well, since we believe it is 
beyond the shadowing dynamics. It might e.g. be caused by bound nucleon swelling. 



Like in ( |189| ) the total cross section is represented in the form 

ctT*^ = Aa^'P - Aa, (192) 
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where Aa is the shadowing correction, 

oo oo 

Aa = -Re / (fb / dzipA{h,zi) / dz2PA{b, Z2) 

—00 zi 

X / da j d?p2^*qq{p2,a)(Tqq{s,p2)A{p2,Zi,Z2,a) 



with 



A{p2,Zi,Z2,a)= I d'^piW{p2,Z2\pi,zi)e '''"'"^^^ "^^ (Tqqis, pi)'^gg{pi,a) 



Here, 



Ql 



Q'^a{l — a) + m'j 



2ua(l 



a] 



(193) 



(194) 



(195) 



is the minimal longitudinal momentum transfer when the photon splits into the qq 
dipole. These equations were first suggested by Zakharov in the context of the LPM 
effect [^. Indeed, as already mentioned in the introduction, shadowing in DIS can be 
regarded in the target rest frame as the LPM effect for pair production. 
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Figure 27: Propagation of a qq- 
pair through a nucleus. Shown 
is the case of a finite coherence 
length, where the transverse mo- 
tion is described by the Green 
function W {p2,Z2\pi,zi). 



The shadowing term in ( |192| ) is illustrated in fig. ^. At the point Zi the photon 
diffractively produces the qq pair (7*A^ — > qqN) with a transverse separation pi. The 
pair propagates through the nucleus along arbitrarily curved trajectories, which are 
summed over, and arrives at the point Z2 with a separation p2. The initial and the final 



T,L, 



separations are controlled by the light-cone wavefunction '^^^ {p,Oi)- While passing 
the nucleus the qq pair interacts with bound nucleons via the cross section agg{s, p) 
which depends on the local separation p. The Green function W{p2, Z2\pi, zi) describes 
the propagation of the pair from zi to Z2, see ( |194|) . One recognizes the diffractive 
amplitude ( p6|) 



fil* -^ m) = i'^qq{pl,a)(^qq{s,pl] 



(196) 



At the position Z2, the result of the propagation is again projected onto the diffractive 
amplitude. The Green function includes that part of the phase shift between the 
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initial and the final photons which is due to transverse motion of the quarks, while the 
longitudinal motion is included in ( p.94| ) via the exponential. 

Thus, (|193|) does not suffer from either of the two problems of the approximations 
( p.77| ), ( |189| ) and ( |191| ). The longitudinal momentum transfer is known and all the 
multiple interactions are included. 

The Green function W{p2, Z2', pi, zi) in ( |194| ) satisfies the two dimensional Schro- 
dinger equation, 

.dW{p2,Z2\pl,Zi) A(p2) 



dZy, 



2ua(l 



a] 



W{p2,Z2\pl,Zi] 



- 2^^' 



P2) PA{b, Z2) W{p2, Z2\pi, Zi] 



(197) 



with the boundary condition W{p2, Zi\pi, Zi) = 5^'^\p2 — Pi)- The Laplacian A(p2) acts 
on the coordinate p2. The kinetic term A/[2i/a(l — a)] in this Schrodinger equation 
takes care of the varying effective mass of the qq pair and provides the proper phase 
shift. Indeed, the term z/a(l — a) can be regarded as the reduced mass of the pair. 
The role of time is played by the longitudinal coordinate Z2- The imaginary part of 
the optical potential describes the absorptive process. The nonrelativistic appearance 
of (|197|) comes from the fact that the energy in light cone coordinates can be written 
as 



P- 



P 



m 



P^ 



:i98) 



One recognizes that the two dimensional Laplacian originates from the Pj and that 
P'^ jl plays the role of the mass in this two dimensional Schrodinger equation. 

The Green function method contains the eikonal approximation (|189|) and the 



Karmanov-Kondratyuk formula (|1 77|) as limiting cases. 

In order to obtain the eikonal approximation, one has to take the limit v 



this case, the kinetic term in ( |197|) can be neglected and one obtains 



W{p2,Z2\puZx)\^_ 



5^'\P2 



Pi) exp 



o-(s,P2) 



dzpA{b,z) 



00. In 



(199) 



When this expression is substituted into (|194|) and with g™" ~^ one arrives after a 
short calculation at the result ( |189| ). 

One can also recover the Karmanov-Kondratyuk formula, when one neglects the 
absorption of the qq pair in the medium, i.e. the imaginary potential in ( |197| ) is omitted. 
Then W becomes the Green function of a free motion, 



W{p2,Z2\pl,Zi)\^^Q 



(fpT 

27r 



exp 



-ipT ■ (P2 - Pi) 



1Pt[Z2 - Zi, 

2ua{l — a) 



(200) 



where px is the transverse momentum of the quark. In this limit, the shadowing term in 
( p,92|) reproduces the second term in ( |1771) . To see this, note that the Fourier transform 
of the diffractive amplitude reads 



fooiPi 



^<^^(p,«)Mp)e'^^-^: 



(201) 
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In momentum space, the mass of the pair is well defined, 

= ^^. (202) 

a{l — a) 

Then one obtains ( |177| ) with this M]^. 

We calculate nuclear shadowing for calcium and lead with the Green function 
method and compare to the approximations (|177|) , and ( |191| ), where only the dou- 



ble scattering term is taken into account. As was mentioned before, only the valence 
gg-part of the virtual photon is taken into account, but the higher Fock components, 
containing gluons and sea quarks, are neglected as well as the effect of anti-shadowing. 



The dipole cross section is approximated by the form (Jqq{p) = C p\ C ~ 3 ||152| , |153 



which works remarkably well, even at large separations and is sufficient for our pur- 
pose. A uniform density pA = 0.16 fm~^ is used for all nuclei, and the quark masses 
are fixed at ttIu = rud = 300 MeV, rris = 450 MeV and rric = 1.5 GeV. Within these 
approximations it is possible to solve ( |197|) analytically. The solution is the harmonic 
oscillator Green function with a complex frequency [|160|| , 



^ ^P'^ "^' ^'^ ''^ = 2vrsinhr.;A.) '"^ i "I 



{pI + pI) cothiuAz) '^P^'P' 



sinh iuoAz) 



(203) 



where 



Az = Z2-Z1, (204) 



2 _ • CpA 



u' = i "" (205) 

ua (1 — a) 

a^ = -i CpAva (1 - a) . (206) 

This formal solution properly accounts for all multiple scatterings and the finite lifetime 
of the hadronic fiuctuations of the photon, as well as for fiuctuations of the transverse 
separation of the qq pair. 

Because of all these approximations, the calculation is not compared to data, but 
only to the standard approximations ( |177|) , ( |189|) and ( |191| ). The comparison with 
data will be presented below in section 



The results are shown in fig. ^. The dashed curves show the predictions of ( |191| ) 
which is called standard approach. The mean values of a^ and a are calculated using 



the perturbative LC wavefunctions (^) and (|^) given on page |T^. The intermediate 
state mass is fixed at M^ = Q"^. At low xbj < 0.01 shadowing saturates because 
q = 2miyXBj -C 1/Ra- The thin solid curve also corresponds to a double scattering 
approximation, but now only the imaginary potential in ( |197| ) was set to zero. The 
formfactor is treated properly, i.e. the kinetic term that describes the relative transverse 
motion of the qq pair, correctly reproduces the phase shift. The difference between 
the curves is substantial. The thin solid curve does not show full saturation even 
at X = 0.001. We conclude that the prescription M^^ = Q^ does not work well for 
shadowing in the transition region, where xbj — > 0.1 and formfactor effects are of 
crucial importance. This could have been expected from the estimate ( |187[ ) which 
shows that fluctuations with quite large masses participate in shadowing. 
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Figure 28: Nuclear shadowing for calcium and lead at different values of Q^. The 
dashed curve is calculated in the standard approach ( P76| ), i.e. with coherence length 
Ic = 1/2'mNXBj and in double scattering approximation. The thin solid curve corre- 
sponds to the double scattering approximation with the free Green function, ( |200| ), and 



the thick solid curve shows the full calculation, (|203|) . 
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The next step is to do the full calculation and to study the importance of the higher 
order rescattering terms in nuclear shadowing. The results are shown by the thick 
solid curves. Higher order scattering brings another substantial deviation (especially 
for lead) from the standard approach. At very low xbj the curves saturate at the level 



given by ( |189| ). Note that the NM collaboration [|161|| was indeed unable to fit the 
A dependence of the shadowing ratio at low xbj with an ansatz linear in the nuclear 
density, indicating the importance of higher order rescattering terms. 

However, calculations with ( |197| ) turn out to be numerically quite involved, if one 



wants to use a realistic parametrization of the DCS and a realistic nuclear density. 
The approximations introduced in this section significantly simplify the calculations, 
because one can find an analytic expression for the Green function W. We are forced 
to apply simplifying approximations also in our further calculations, but will use an 
effective C in the dipole cross section in order to account for the correlation between 



the energy and p dependence, as will be described in section p.7| . Then one can still 
make use of the solution ( ^U^ ) of the Schrodinger equation. This approach can be 
generalized to include also higher Fock components of the photon ||12||. 



3.3 Derivation of the formula for nuclear shadowing in DIS 

Before we start with the derivation, we briefly sketch the idea of the calculation. The 
applied approximations are summarized at the end of this subsection. For a nuclear 
target, the total cross section may be written in eikonal form ( |189| ), if the transverse 
separation p is frozen, i.e. at very small xbj- The eikonal formula can be obtained in 
exactly the same way as by Glauber |^ in nonrelativistic quantum mechanics, because 
of the nonrelativistic structure of the light-cone Hamiltonian. In the eikonal approxi- 
mation, one splits a fast oscillating phase factor from the wavefunction (f{r) = e'^^^(p{r) 
and obtains for the slowly varying part (p an equation of the form 

In the target rest frame, the interaction is given by a color-static potential V . However, 
anticipating that all dependence of the interaction will be absorbed into the dipole cross 
section aqg{p), we use an abelian potential. Strictly speaking, this is justified only in 
the case of an electron-positron pair propagating in a condensed medium, but since 
the influence of the potential will be expressed in terms of scattering amplitudes and 
because our final result reproduces well known approximations in limiting cases, we 
believe, that the results holds also for the case of a nonabelian potential. 

The Laplacian acts only on the transverse coordinates and is omitted in the eikonal 
approximation. Then, ( |207|) is easily integrated. The multiple scattering series is sum- 
med like in Glauber theory P). This is possible, because the typical distance between 
two nucleons inside a nucleus is roughly 2 fm, while the gluon correlation length is 
much smaller, presumably ~ 0.3 fm. After averaging over the medium, one obtains the 
eikonal expression ( |189| ). 

The idea is now to keep this Laplacian, because it describes the transverse motion 
of the particles in the pair. Taking into account this motion, one can correctly describe 
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the effective mass of tlie fluctuation in a coordinate space representation. The phase 
shift function has then to be replaced by the Green function for ( P07| ). We write this 



Green function as a path integral and averaging over all scattering centers yields an 
effective Green function W which fulfills the two dimensional Schrodinger equation 
( |197|) with the optical potential Vopt = — icgq (p) Pa ip, z) /2. 

In order to derive explicit expressions for the terms in (|192|), we start from the 



general expression for the cross section for production of a gg-pair, 

(27r) z/2p^2py 

/ em I M (27r)'4z/2A(l-A) ^ ' 

where Zj is the flavor charge, and aem = 1/137. Note that we prefer to call the 
longitudinal momentum fraction of the quark by A, because a is already used to denote 
many other quantities. The matrix element is given by 

Mfi = f d^r *5 (r) a ■ e'^g- (r) e''''. (210) 

The 2;-axis is chosen in direction of the propagation of the photon. The photon's 
momentum is denoted by k and it's energy by u, Pg is the momentum of the quark 
and Pq the momentum of the antiquark. In ( p09| ) we have introduced the energy 
fraction A = Pg/i^- In the ultrarelativistic case which is considered, pair production 
takes place predominantly in forward direction and therefore we distinguish between 
the longitudinal direction (^-direction) and the transverse directions. The longitudinal 
momenta are large compared to the flavor masses, Pz ^ mf, and to the perpendicular 
momenta, \p±\ ~ m'l. In our approximation, terms of order ruf/pz are kept only in 
exponentials and neglected otherwise. All higher order terms are omitted. 

Note, that the pair is created electromagnetically in a color singlet state, therefore 
the factor Zfaem appears in (|208| ), but the quarks interact with the nucleons as de- 



scribed by a phenomenological, scalar potential 0(r). The particles in the pair move 
in a potential U (r) that is a superposition of the potentials of all nucleons, 

A 

U{r) = Y,<Pir-f^). (211) 

i=i 

The vector fj runs over all positions of the nucleons. 

The wavefunction of the quark fulfills the Dirac equations and is an eigenstate 
with positive energy, while the antiquark is represented by an eigenstate with negative 
energy, 

(p;-f/(r)-m//3 + ia-V)^,(f) = 0, (212) 

(-pO-f/(r)-m//5 + ia-V)^g-(f) = 0. (213) 

No interaction between the quark and the antiquark is taken into account and there- 
fore, the two equations decouple. The wavefunction of the quark \E'q (r*) contains an 
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outgoing plane wave and an outgoing spherical wave in it's asymptotic form, while the 
wavefunction of the antiquark contains an incoming spherical wave and an incoming 
plane wave. These equations are transformed into second order equations by applying 
the operator (p° — U (r) + rrifP — id ■'V) on the first equation, and the corresponding 



operator on the second equation, as described in ||163|| . When the term quadratic in 
the potential is omitted, one obtains 

(A+|p;|'-2p;f/(f) + i«-(Vf/(f)))%(f) = 0, (214) 

{A + \pgf + 2pp{f) + ia-{VU{f)))^q{r} = 0. (215) 

The solutions may approximately be written as Furry ||164| , |165|| type wavef unctions. 



*. (f) = e'^^^-^ (^1 - ^ . v) F, (f) u (p„ A,) , (216) 

*, (r) = e-'^^-^^ (^1 + i^ . v) F, (f) V {p,, A,-) . (217) 

Here, u {pg, Xq) is the free spinor with positive energy and polarization A^ which fulfills 
the equation (|Sg— nif) u {pq, Xq) = and similarly {-^qi- rrif) v {pq, Xq) = 0. In Dirac 
representation the spinors read 



U [Pq, Xq) = \lp^q + Ul f | ^ ' Pq ^ I . ^ {Pq^ \) = \/ P^ + ^f I P? + ^f ^^ 

Pq + rnj " J \ Xq 

(218) 

The three Pauli spin matrices are denoted by a and the Pauli spin state referred to the 
rest frame of the particle is Xg, or Xq respectivly. This means explicitly s- aXq = XqXq 
and s ■ a Xq = —XqXq with a spin vector s normalized to unity and Xq, Xg = ±1. The 
functions Fq and Fg have no spinor structure any more. They contain all dependence 
of the potential and have to be calculated for a given U {"f). They fulfill the equations 



{A + 2ipg-V-2plU{r})Fg{r} = 0, (219) 

{A-2ipg-V + 2pp{r})qj{r} = 0, (220) 

with boundary conditions F — * 1 for the quark and for the antiquark as z -^ cxd. Note, 
that it is essential to take the correction proportional to a in ( |216|) and ( |217|) into 



account, although these terms seem to be suppressed by a factor of l/p°. It turns out, 
that when one calculates the matrix element of the current operator, [a ■ €) r-, between 
free spinors, the large part cancels and one is left with a contribution of the same order 
as produced by the correction term. Thus, the terms proportional to a may not be 
neglected in the matrix element, although they give only small corrections to the wave 
functions. 

In order to remove the dependence on the transverse momenta from the phase 
factors in eq. (216) and ( 217 ), we rewrite the solutions in the form 



^qirl = e^'^^^'^(l-^-V-|^-J?x,, + 0(|p;|-p.,,))^,(r)«(j>,,A,) 
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(221) 

(222) 

with 

V^Jf) = e'P^'«-^^^e-'('P^«l-P^-')"Fg(r), (223) 

^-(f) = e-'P^-^-''^e'^\P'i\-P^'i^'Fq{r}, (224) 

and \p\ = 4 /po — rnl for the quark and the antiquark respectively. In the following, we 

neglect the terms with a^ in eq. ( p21| ) and (|222|) , because they are of order O (l/PqP^). 
The functions ipq (r) and ipq (r) will play the role of effective wave functions for the 
quarks. 

The phase factors combine in the matrix element (|210|) to the minimal longitudinal 
momentum transfer, 

g— = k - \pq I - \pq \^^ + ^ + ^, (225) 

2z/ 2po,q 2po,g 

and one obtains 



X 



u' (pg, Ag) ( a ■ e +a — ^ a-e +a-ea — ^ J v [pg, Xgj 



X C (n) ^5 (^"5) _ • (226) 



r=ra=ra 



Here, the coherence length, 1^°-^ = l/q"^^^, enters as an oscillating phase factor. How- 
ever, [^"■^ does not depend on the transverse momenta. Their influence on the cross 



section is encoded in the rest of the wavefunctions, ( p21| ) and ( |222| ). The operator V (fg) 
acts only on the variable fg and the operator V (fg) only on fg. After the derivatives 
have been performed, the whole integrand has to be evaluated at r = -Tg = fg. 

With the representation ( pi8| ) one obtains after some algebra within the demanded 
accuracy 



00 



00 



Mf, = I dze<'"*^ / d'r^ 



VA(l-A) 
X xU T^f^ ■ er + i (1 - A) a ■ e-^ex ■ V (fj_,g) + iAa ■ e^eV ■ V (r_L,g) 

+ (1 - A) (e, X er) ■ V (f^,g-) - A (e, x e^) ■ V (f^,g) 

+ 2gA(l-A)|x,- 
X r,{r^,,,z)^ljg{r^,g,z) . (227) 
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Some details of the calculation can be found in | p.66| ]. The unit vector in z-direction 
is denoted by Cz- The polarization vector eV corresponds to transverse states of the 
7*, while the last term in the curly brackets is due to longitudinal polarization. There 
are two remarkable aspect concerning this last equation. First, it does not contain 
any derivative with respect to z any more, because of the transverse nature of the 
polarization vector e^ and of CzXex- Second, all dependence on the transverse momenta 
in the spinor part has canceled. 

From ( |219| ) and ( |220[ ) and from the definition of ip, ( |223D and ( ^24| ), an equation 



for ip is obtained. Assuming that these functions are only slowly varying with z, we 
omit the longitudinal part of the Laplacian and get 

^§-^%ir^,z) = (^-^ + U{f^,z)^^^{f^,z), (228) 

i^^5(r±,z) = (^_+U{f^,z)\ijg{r^,z). (229) 

Our ansatz yields two dimensional Schrodinger equations, where the 2;-coordinate plays 
the role of time and the mass is given by the energy. The Laplacian A^ acts on 
the transverse coordinates only. The functions ipq{r±,z) and ipq{r±,z) become two 
dimensional plane waves for 2; ^ oo, up to a phase factor that cancels in the square of 
the matrix element. It should be mentioned that the kinetic energy for the antiquark 
has a negative sign. This results from the fact that solutions of the Dirac equation 
with negative energy propagate backwards in time. The Laplacian in (|228[) and 



account for the transverse motion of the pair in which we are especially interested. The 
functions tpq (rj_, z) and ipq {f±, z) in the matrix element (|227|) may now be expressed 
in terms of the Green functions for ( |228| ) and p29| ) and their asymptotic behavior, 

V',(r-1,2,^2) = jdh^,^Gq{r^,2,Z2\f±,l,z^)e'P^^-■'^^■^e-'^\^^^\-P^"'^'-', (230) 
^5(^-1,2,^2) = JdWiGq{r^,2,Z2\r±,i,z^)e-'P^-r^^^'^e'^\P^-\-P-^^^^ (231) 



Now one uses the expression for rfcr'^*^, ( |208|) , and with the matrix element ( |227|) 
and the last two relations, we obtain 



dz / ciz'e^'?^'"^^-^') / rfVx Srq d\ / dVl d\ d 

00 J —00 J J 

X O {f±^q, f^^q) Gg {f±^g, Z I fq, Z^c) Gq {f±^q, Z \ fq, Z, 

X O (^_L,g,T_|_gj Gg (r_|_ g, 2; \Tg,Zoo) Ggyrj_g,z \Tg,ZQo) 



r±=r±„=rj_q 



X e'l^q~'^9j-PX,9 gH^?~'^9j'P-L.5 



X dY,qd^P^,dX 

(27r)'' A2(l-A)'4i/2 



For convenience the operator 



O {r_L^g, r_L,q) = xl \ fnfff ■ e^ + i (1 - A) <? ■ e^e^ ■ V {f±^q) + iXa ■ e^eV ■ V (rl,, 
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+ (1 - A) (e, X er) ■ V (r=l,g-) - A (e, x e^) ■ V (f^,,) (233) 

+ 2Q\{l-\)a-eAxg (234) 



is introduced. In order to obtain the total cross section, we integrate over p±^q and 
p±,q- The exponential factors give a (5-function which makes it possible to perform the 
integrations over all the rs. Note, that O does not depend on r. The result is 

1 oo oo 

at^ = Z) ^ 2Re | ^^^^. J dz J dz'e---^^'^ J d\, J £r', 

— oo z 

X 0(r-l,„fx,,)0*(rl,,,rl,,) 

X Gq{r'^ z'\r^^q,z) Gl{r'^^q,z'\rx,q,z) ^ ^ ^ ^,(235) 

Instead of four propagators one is left with only two, because of the convolution relation 

G (rl,2, ^2 I n_,i, ^i) = / c?V_L G (rl,2, ^2 I n_, z:) G (rl, z | fj_,i, 2:1) . (236) 

In order to derive an expression that is convenient for numerical calculations, one 
makes use of the path-integral representations of the propagators in (|235|). They read 



Gir'^,z'\r^,z)=Jvrexphjd^(±^f'-U{f,o] 
^ jVr exp <^ ij ±\r^di - i J^ X (f {z,) - f,,^) Q {z' - z,) Q 



(237) 



Zj - z] 



(238) 



where the upper sign corresponds to the quark and the lower to the antiquark. In 
this expression, r is a function of ^. The derivative with respect to ^ is denoted by r. 
Furthermore, the condition 

G(rl,z'|rl,^) =8^^\r'^-r^) (239) 

z'^z 

has to be fulfilled together with r (z) = rl and 7^(2;') = r^. In ( |23^ ), the phase shift 
function X(r — f^-^x) is introduced. The step function 6(2;) is for z < and 1 for 
2; > 0. As mentioned before, U is the superposition of all potentials of the nucleons, see 
eq. (|211|) . Let be the potential of a single nucleon in the target. The position of the 



nucleon with number j is denoted by the transverse vector f^-^x and the longitudinal 
coordinate Zj. If the range of interaction is much smaller than the distance z' — z^ 
the potential is practically zero outside the domain of integration over ^ and the phase 
shift function is given by 

X (rl - r,, x) = / di^ (fx - r-, X, i-z,). (240) 
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Note that t(^) was replaced by t{zj). This means, only an average value of the 
transverse coordinate is used for calculating the phase shift for scattering on a single 
nucleon. 

The two path integrals ( |238| ) for the quark (+) and the antiquark (— ) sum over all 
possible trajectories of the two particles. In order to calculate the cross section for pair 
production in the nuclear medium, one has to average over all nucleons. One obtains 
with the path-integral (|238|) 

{Gg {fig, z' |rl,g, z) Gl (f| ,j, z' |fx,g-, z)) 

z' 

r /' rP ■ rP- . ')\ 

(241) 




+i 1^ e (z' - z,) e {z, - z) (X {f, {z,) - rV) - X (f, (z,) - r,,^)) I \ (242) 

with boundary conditions Tq {z) = f±^q, fg {z') = f^^, fq {z) = fj_^q and fq{z') = f'^^. 
The averaging procedure is similar to the one described in [^ . All correlations between 
the nucleons are neglected and the average nuclear density pA is introduced, which is 
normalized to A, ( P43|) . Then, the whole expression may be written as an exponential, 
if A is large enough. 



X {Tq {zj)- r,- . 



i=i 



' A 



exp ( i ^ e (2' - Zj) e {zj - z) (X {Tq {zj) - r,- i 

z' 
^-jdhjdi'pA{s,i') 

z 

(l-exp(i(x(r,(e')-^1-X(r,(r)-^"'))))} (243) 

X |rf2.(^l-exp(^i(x(r,(r)-^1-X(r,(r)-5'))))}. (244) 



When pA is varying fairly smoothly inside the nucleus, one can replace its dependence 
of s, the transverse distance between the pair and the scattering nucleon, by the impact 
parameter 6, eq. (|244|) . Since we have a short ranged interaction, it is reasonable to 



choose as impact parameter h = {f± — f[_) /2. This way one finds the forward scattering 
amplitude for a dipole scattering on a single nucleon, see eq. (|244|) . The corresponding 
cross section, 



aqq (p) = 2Re f (fs (l - exp A (x (p - s) - X {s)) \ \ 



(245) 
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appears as imaginary potential in the propagator. At this point one notices that the 
forward scattering amphtude in the abehan model is predominantly real, while in a 
nonabelian theory it has to be imaginary. Nevertheless, we believe that the final result 
holds also for a nonabelian theory, because all dependence on the potential of the 
nucleons, (r) , is absorbed into a^g, for which we will use a parametrization. One 
finally arrives at 

{Gg (f| ,j, z' \f±^g, z) G\ {rlg,z' \r^^g, z) ) 
= |pf,|pf,exp|i|rfe(|f^' + fn' + ^P^(^0^,.-(|r,-r5|)^ i .(246) 

It is convenient to introduce center of mass coordinates, fVe/ = Tg — -fg and Tcm = 
{1 — X) Tg + \Tg, and to express the sum of the two kinetic energies as the sum of the 
kinetic energy of the relative motion and the center of mass kinetic energy, 

{Gq {r±,q, z' \r±,q, z) G*g {f'^^g, z' |fx,5, z)) 

= / T>Tcm / VTrel exp ll d^ ( -f^,; + -f^^ + -pA (&, ^99 i'^rel) j > (247) 



exp 1— 



• ^ (A Kg - r±,q) + (1 - A) {rig - f^,g)y 



2Tri{z' - z) \ 2 z' -z 

X J VTrel exp < i y d^ (^^rel + \pA (6, ^^qq {TrelYj \ ■ (248) 



We have introduced the reduced mass 

1 = 1 + 1 = I (249) 

/i P'q P'q ^A(l-A)- ^ '^ 

Since the imaginary potential depends only on the relative coordinate, the center of 
mass propagates freely and what remains is the effective propagator 



W (p', Z'\p,z) = / VTrel GXp I 1 d^ ( 2 ^™' " ^"P* ^^' ^^'^'' ^V 

with p' = r'j_ ^ — f'j_ and p = f±^g — r^ ^ and the optical potential 



(250) 



Vopt (b, p, z) = -^PA (b, z) agg (p) . (251) 



It fulfills the equation 

W {p',z' \p,z) = i6{z' - z) 5^2) (p" _ p) . (252) 



.^aJ^+ 2.A(l-A) -^°^^^^-^'^ 
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The propagator for the center of mass coordinate produces a (5-function and thus, one 



obtains from (|235|) after averaging over the medium 



a- 



7*A 
tot 



Zjaer 



(fh 2Re 



dX 



oo oo 



A2(l-A)^ 



dz / rf^'e^<'"(^-^ ) 



with 



X 0,,{p)Ol,{p') W{p',z'\p,z) 



O,, (P) = Ol (p) + O]-, (p) 



p'=p=0 



where the transverse part of the operator is 



O^q (p) = xj 1^/'? ■ ^T + i (1 - 2A) a ■ e^er ■ V (p) + (e, x eV) ■ V (p) W,- 



and the longitudinal part 



Ot- (p) = x;2 Q A (1 - A) a ■ e^x, = 2 Q A (1 - A) 6.x,. 



(253) 
(254) 

(255) 
(256) 



These operators are already known from the definition of the LC wavefunctions in 
section p73| , see ( [70|J7T| ). Equation ( |253| ) is the central result of this section. It is the 
total cross section for production of a gg-pair from a virtual photon scattering on a 
nucleus. We have not summed over the spins of the quark and the antiquark and not 
averaged over the polarizations of the photon. We have not summed over the different 
flavors, either. The expression for the operator O, ( |254| ) , depends on the spin vector of 
the quark and the antiquark. The directions of these vectors may be fixed arbitrarily. 
In order to represent the result in the form ( |192| ), W and it's complex conjugate is 
rewritten in an expansion. The results can be combined in the following way: 

W{p\z'\p,z) = Wo{p',z'\p,z) 

z' 

+ i dzi d'^piWo (p', z'\pi, zi) Vopt {b, pi, zi) Wq (pi, zi \ p, z) 

z 
z' Z\ 

- dzi d^pi / dz2 I d^p2Wo (p', z'\pi, zi) V*pt (6, pi, zi) 

Z z 

X W (pi, zx I p2, z-i) Vopt {b, p2, Z2) Wo (p2, Z2 \ p, z) . (257) 



Here, Wq is the propagator corresponding to ( p5CI| ), when the potential is absent. The 
first term gives a divergent contribution, which is the wave-function renormalization 
for the photon. The second term leads to the first contribution in ( |192| ). The operators 
Oqq applied to the free propagator Wq give the light-cone wave functions \E'^^ (p. A) , up 
to a constant overall factor. The third term in the above expansion is the interference 
term. Since this term contains the full propagator, there are no higher terms in this 
expansion. 
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As an example for further calculation, consider the integral 



oo oo 



I 



dz I dz'e"^"^""^'-''^ I dz, I d^p, / dZ2 / d^P2 



-oo z 



X Kpt iP^ Pi) ^l) Vopt {b, p2, Z2) 

X Wo (p", z' I Pi, zi) W (pi, zi I P2, Z2) Wo (p2, 2:2 I p, 2;) 



which is needed to calculate the interference part 



cr. 



int / 



zW. 



tot 



4z/2 



d^b 2Re 



dA 



AMl-A)- 



rO (p-) O* (p") X 



p'=p=0 



(258) 



(259) 



With the new variable e^ = A (1 — A) Q^ + rri^f, which is also defined in ^, one finds 



„rmn 



2z/A(l-A)' 



(260) 



and because of the relation 

'i^ + ^^[frf^] [Wo (p", .' I p, .) e--r"(^'-^)) = 15 (,' - ,) 5(2) (p- - p) , 

(261) 
one can write the propagator in the form 



Wo{p',z'\p,z)e 



-iq^^"[z'-z) 



(27r)' J 27r 



a; 



P,+e' 



2z/A(l- A) 



iO_ 



(262) 



The 0+-prescription for the pole in the complex tu-plane ensures that Wo (p', z' \p,z) 
for z > z' . Putting ( p62| ) into ( 258 ) yields after a short calculation 



J 



4Z/2A2 (1 - A)^ 



00 21 



dz, / d^2e'^2""(^^"^^) / rfVi / rf'p2 



-00 —00 



(2vr)^ 

X K*pt (^) Pi^^i) Vopt (b, p2 , 2:2) 

X Ko (5 \p' - pi\)W (pi, zi I p2, 22) Ko (5 |p2 - p1) . 

Kq is the MacDonald function of zeroth order. We have used the relation 



(263) 



(264) 



The next step is to insert (|263|) into (p59D and to calculate the contribution from the 
second term in the expansion (|257|) in a similar way. One obtains for the total cross 
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section 



7* A 



AZ 



2 "em 



dX / (fpiagg{pi) 



0,,{p)Ko{e\p-pi\) 



oo oo 1 



Zj -^ 2Re / d% / dz, / dz2 I dX / d'pi / d''p2 e 

— oo 2l 



p=0 



-ig™"(22-Zl) 



p=0 



X v;;, (6, pi, 2;i) Kp, {b, P2, ^2) (^o;- (pO Ko {e \p' - pil) 

X W{p2,Z2\pi,Zi) (Oqg{p)KQ{e\p-p2\) 

With help of the relation 

V ip) Ko (ep) = -e^K, (ep) 
P 

one finds the light-cone wave functions 



p'=0 



^?,Kp1 = ^/^c?,^,(p)KoM 



27r 



^/ ^^^ "^ Ko {ep) 6x„x, 6x„x, 



+ HAg (2A - 1) er ■ ep + (eV x e;) ■ Cpj eKi (ep) 5a, _a,- [ 



and 



^i-(«,p) 



gg 



^/^0i,(plKoM 

zvr 



Z^V^2QXil-X)Ko{ep)5x„x, 



2ti 



q,Aq 7 



(265) 
(266) 
(267) 

(268) 

(269) 
(270) 



with the Kronecker-5. The unit vector in p-direction is denoted by Cp. As spin vector 
we have chosen the unit vector in ^-direction and the parameter A^ takes the value +1 
for spin in positive z-direction and the value —1 otherwise. For the antiquark it is vice 
versa. For the photon, A-y = +1 for positive helicity and A-y = — 1 for negative helicity. 
The transverse light cone wave function has one part that depends on Ko and another 
part dependent on Ki, the MacDonald function of first order. Note the different spin 
structures of these parts. In the Ko-part, the spins of the quarks add up to the spin 
of the photon, but in the Ki-part of the transverse LC wave function, the spins of the 
quarks add to and the pair gets an orbital angular momentum. We finally sum over 
all fiavors, colors, helicities and spin states and get the following expression: 



a 



7* A 
tot 



A dX d^p agg (p) 



vl/^ 



qq 



(a,p)| 



1^-^ 



qq 



{a,p)\ 
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N 

3a 



^^ZJRe d% / dzi / dz2 dX d^pi j d 



'2p2 e-igr"(^2-zi) 



(2^) /=! 

•' — OO 21 

X Pa (6, zi) Pa {b, Z2) a^q (pa) <yqq (pi) 

X /(l-2A(l-A))£2^i:faKi(£pi)Ki(£p2) (271) 

I P1P2 

+ (mj + 4g2 a2 (1 - A)') Ko (epi) Kq (eps) W (p2, ^2 | Pi, ^1) • 



■,r,L/ N 12 



Here, 1^^,^ (ep) p are the absolute squares of the transverse and the longitudinal light- 
cone wavefunctions summed over all flavors, see (^51) and (^61) on page 0. This form 
is more convenient for numerical calculations than ( p53| ) and was used in [^ for a 



calculation of nuclear shadowing. ( P71| ) was for the first time suggested in a paper by 
Zakharov p7|. 

In the rest of this subsection, the assumptions and approximations entering this 
derivation are summarized. The starting point is the Dirac equation with an abelian 
potential. This simplification is used in anticipation that all dependence of this po- 
tential will be absorbed into the dipole cross section. Because the final result contains 



known approximations as limiting cases, see section |3.2| , we believe that it holds also 
for the case of a nonabelian interaction. 

Furthermore, no interaction between the quark and the antiquark is taken into 
account and therefore, the two Dirac equations decouple. The Furry- approximation 



| 164| |165|| of the wavefunctions is employed which is known to be a good approximation 



to the continuous spectrum of the Dirac equation for high energies. In the next step a 
two dimensional Schrodinger equation for a scalar function which may be regarded as 
an effective wavefunction is derived. The 2;-coordinate plays the role of time, since the 
particles move almost with the velocity of light. This Schrodinger equation is solved 
in terms of it's Green function. Averaging over all scattering centers in the nucleus 
yields an optical potential that is proportional to the total cross section for scattering 
a gg-pair off a nucleus. The real part of the forward scattering amplitude is omitted. 
All dependence on the potential is absorbed into the dipole cross section. 

The averaging procedure and the summation of the multiple scattering series is 
similar to the one in Glauber theory and most of the approximations come in 
at this point. First, all correlations between the nucleons are neglected. Then, the 
influence of the potential is described by a phase shift function. It is also assumed, 
that the interaction is short ranged and the pair interacts only with one nucleon at a 
given time. The phase shift for scattering a particle in the pair off a single nucleon 
is calculated for an average value of the transverse coordinate of the particle. This 
means, the transverse coordinates should not vary too rapidly within a longitudinal 
distance of the order of the interaction range. In order to obtain an exponential from 
the averaging procedure, the nuclear mass number A has to be large enough. Further 
approximations are, that both particles in the pair see the same nuclear density. The 
value of the density in the middle between the quark and the antiquark is used in the 
calculation. Furthermore, the motion of the center of mass of the pair is approximated 
by a free motion, since the pair is scattered predominantly in forward direction. 
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One finally arrives at the result ( f253| ). This formula allows to calculate the cross 
section cr7ot for arbitrary polarization of the photon and the pair. However, this 
equation is not convenient for numerical calculations and is modified by introducing 
the light-cone wavefunctions ( |265| ). Finally, one sums over all helicity and spin states, 
arriving at (P71|). 



3.4 The mean coherence length 

We emphasize again that nuclear shadowing is controlled by the interplay between two 
fundamental quantities. 

• The lifetime of photon fluctuations, or coherence time. Namely, shadowing is 
possible only if the coherence time exceeds the mean internucleon spacing in 
nuclei, and shadowing saturates (for a given Fock component) if the coherence 
time substantially exceeds the nuclear radius. 

• Equally important for shadowing is the transverse size of the hadronic fluctuation 
of the virtual photon. In order to be shadowed the fluctuation has to interact 
with a large cross section. As a result of color transparency [^, ^, ^, small size 
configurations interact only weakly and are therefore less shadowed. For the qq 
Fock component, the dominant contribution to shadowing comes from the large 



aligned jet configurations [B4, 140| of the pair. 



In this section a definition of the mean coherence length or the mean fluctuation 
lifetime of a Fock state, relevant for shadowing, is proposed. The mean coherence time 



for the qq Fock state is evaluated using the perturbative (^, |46| ) and nonperturbative 
{ \n\ , |78D wavefunctions. It is observed that the coherence length is substantially longer 
for longitudinal than for transverse photons. At the same time, both are different from 
the usual prescription /^ = {2m]\fXBj)~^- At high Q^ one approximately has 

(272) 



5m]sfXB 
4 



J 



'^ bniMXBj ' 

for transverse and longitudinal photons respectively. The coherence length is found to 
vary steeply with Q^ at fixed xbj and small Q^. The coherence length for the \qqG) 
Fock component, which controlles nuclear shadowing for gluons, is also calculated. The 
latter turns out to be much shorter than for \qq) components, see fig. |2^. Therefore, 
the onset of gluon shadowing is expected at smaller xbj than for quarks. 

A photon of virtuality Q^ and energy v can develop a hadronic fluctuation for a 
lifetime, 

^^ = n^rV = -^- ' ^274) 

Q-^ + M^g XBj rriN 

where Mgg is the effective mass of the fluctuation, and the factor P^^ = (1 + M^g/Q"^). 
The usual approximation is to assume that Mgg ^ Q^ since Q^ is the only large 
dimensional scale available. In this case P = 1/2. 
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The effective mass of a noninteracting gg-pair is well defined, M^^ = {'mj+p'^)/a{l — 
a), where Pt and a are the transverse momentum and fraction of the light-cone mo- 
mentum of the photon carried by the quark, respectively. Therefore, P has a simple 
form. 



P{kT,a) 



Pt + £'^ 



(275) 



where 



£^ = a(l - a)Q^ + mj . (276) 

To find the mean value of the fluctuation lifetime in vacuum one should average ( |275| ) 
over Pt and a weighted with the wavefunction squared of the fluctuation. 



{P\ 



^ 



qq 



P{kT,a) 



* 



qq 



^ 



qq 



* 



qq 



(277) 



The normalization integral in the denominator in the r.h.s. of ( p771 ) diverges at 
large pt for transversely polarized photons, therefore one arrives at the unexpected 
result {P'^)vac = 0. The divergence is related to the wavefunction renormalization 
constant 

I-Z3 (278) 



^ 



qq 



^ 



qq 



and can be interpreted as a result of overwhelming the fluctuations of a transverse 
photon by heavy qq pairs with very large Pt- Such heavy fluctuations indeed have 
a very short lifetime. However, they also have a vanishing transverse size p ~ 1/pt 
and interaction cross section. Therefore, such fluctuation cannot be resolved by the 
interaction and do not contribute to the DIS cross section. To get a sensible result 
one should properly define the averaging procedure. One is interested in the fluctua- 
tions which contribute to nuclear shadowing, i.e. they have to interact at least twice. 
Correspondingly, the averaging procedure has to be redefined as, 



(P) 



shad 



(/(7* -^ qq) P{pT,(y) f{-f* -^ qq)^ 



f{Y -^ qq) f{Y -^ qq) 



(279) 



where /(7* -^ qq) is the amphtude of diffractive dissociation of the virtual photon on 
a nucleon 7* p ^ qqp, 



Therefore, P has to be weighted with the interaction cross section squared cr^^(s, p) 



in the averaging procedure. Then, the mean value of factor P{a,pT) reads. 



-pTS^ 



T,L 



j^ da j d^ pid^ P2 ^qq {p2,a) a^g- (s,p2) ^ (P2 - Pi, «) %^ {pi,a) agg{s, p 



with 



P{P2 - Pi,a) 



j^dajd^p ^lq^{pa)aqq{s,p) 
d^PT 



{2ny 



exp {-ipT • {p2 - pi))P (a, p)- 



(280) 
(281) 
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Using expression (|275D one obtains for a non interacting gg-pair, 



P (p2 - Pi, a) = — Ko [e \p2 - pi\ 



(282) 



As a simple estimate for the mean value (|28CI|) one can use the small-p approximation 



for the dipole cross section aqq{s,p) = C{s)p . The Factor C{s) does not enter the 
result since it cancels in ( P80| ). One obtains for transverse and longitudinal photons 
with perturbative wavefunctions 



[P-) 



2Q2 J da {1- a) a 




/t/a(l-a)2a2 /^6 



(283) 



(284) 



respectively. 

The factor {P'^'^) is calculated as function of Q^ at xbj = 0.01 from (|283|) and 
( P84|) . The results depicted in fig. |2^ by dotted lines are quite different from the naive 
estimate P^'^ = 1/2. Besides, P^ turns out to be substantially longer than P^. This 
indicates that a longitudinally polarized photon develops lighter fluctuations than a 
transverse one. Indeed, the effective mass Mqq is maximal for asymmetric pairs, i. e. 
when a or 1 — a are small. However, such fluctuations are suppressed in longitudinal 
photons by the wavefunction (^61). 

The dependence of (P^'-^) on xbj depicted in fig. ^ for Q"^ = 4 GeV^ and Q"^ = 40 
GeV^ is rather smooth. Therefore, the coherence length varies approximately as Ic oc 

1/xBj- 

The simple approximation aqg oc p^ is not realistic since nonperturbative effects 
affect the large-p behavior. A more realistic parametrization of the dipole cross section 
( p.57| ) is introduced in section p?7[ 

Note that ( P80| ) can be represented in the form. 



(p.,.) 



N 



T,L 



DT,' 



(285) 



The angular integrations in ( |280|) for the denominators D^'^ are trivial and for the 
numerators N'^'^ one uses the relation P7[ , 



Ko {e |pi - P2I) = Ko (£p>) lo {ep<) + 2 ^ e^^^K^ (£p>) U {ep 



<) ' 



(286) 



m=l 



where p> = max(pi,p2), P< = min(pi,p2), cos0 = pi ■ p2/(piP2) and \m{z) are the 
modified Bessel functions of first kind (Bessel function of imaginary variable). It is 
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Figure 29: The dependence of (P) = IcIlT"" ^ defined in (|7|), on Q^ at xbj = 0.01. 
Here, /^"^ = l/rriNXBj is the largest possible coherence length. The curves are for qq 
fluctuations of transverse and longitudinal photons, and for qqG fluctuation, from the 
top to bottom, respectively. Dotted curves correspond to calculations with perturbative 
wavefunctions (HHHBI) and an approximate dipole cross section oc p^. Dashed curves are 



the same, except the realistic parameterization ( |157|) is employed. The solid curves are 
calculated with the nonperturbative wavefunctions (7^,7^). Note that as Q^ becomes 
to low at fixed Xbj, no high energy approximation can be applied any more. 



clear from this relation that after angular integration only one term in the sum gives 
a non-vanishing contribution. One finally obtains for transverse photons 



K 



T 



2Q'^ daa (l-a) dp2 p2 / dpi pi {mJKl (epa) Kq (spi) Iq (epi) 



[a^ + (1 - af] e^ Kj (erpa) Ki (epi) Ii (epi) } Uqq (s, pi) cr^^ (s, pa) , (287) 



1 oo 



Dl= I da I dp p {m) Kl {ep) + [a^ + (1 - af] e" Y.\ (ep) } a% (s, p) , (288) 
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Figure 30: xbj dependence of (P) = lc/l'^°'^ defined in ( |274|) , corresponding to the 
cofierence length for shadowing of transverse and longitudinal photons and gluon shad- 
owing, respectively. Here, 1^°"^ = l/niNXBj is the largest possible coherence length. 
Solid and dashed curves correspond to Q^ = 4 and 40 GeV^ and are calculated with 
the nonperturbative wavefunctions (|77|,|7Sp. The realistic parametrization of the dipole 
cross section (|157|) is used. The dotted curve shows the factor (P) for the qqG Fock- 
state. In this case the perturbative LC wavefunction (^) and an approximate dipole 
cross section oc p^ are used. The dotted curve is independent of Q^. 



and for longitudinal photons. 



m 



= 2 Q'^ da a^ {1 — a) / dp2P2 / dpipi 



X Kg {ep2) Ko (epi) Iq (epi) agg {s, pi) agg {s, P2) 
00 
Di = jdaj dppa^ (1 - af K^ (ep) a% (s, p) . 



(289) 



(290) 



The factor {P^'^{x, Q"^)) calculated in this way is depicted by dashed lines in fig. ^ 
as function of Q^ at xbj = 0.01. It is not much different from the previous simplified 
estimate demonstrating low sensitivity to the form of the dipole cross section. 

Although the quarks should be treated perturbatively as nearly massless, at the 
endpoints a or 1 — a —> the mean qq transverse separation p ~ I/5 becomes huge 
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~ l/fTLf and nonperturbative effects are important. A more sophisticated treatment of 
these effects is developed in [|12|, where an interaction between the particles in the pair 
is introduced, as explained in section |2.3| . In the following, we will perform the same 
calculations again, but now with the nonperturbative LC wavef unctions. 

Note that in px representation the free Green function G°^(2;i,pi; Z2,P2) integrated 
over longitudinal coordinate is simply related to the coherence length ( [^74| ), if one 
performs an analytic continuation to imaginary time, z — >• —iz 



dZ2G°Jzi,pi;Z2,p2) 



d^PT_ 

(2 7r)2 



exp 



-ipT ■ (p2 - Pi) lc{pT,a 



(291) 



This is easily generalized to include the nonperturbative interaction. Then, making use 
of this relation one can switch in ( [^8U| ) to p representation and express the mean coher- 
ence length via the Green function. One arrives at new expressions for the functions 
A^^'-^ and D'^^^ in 



N 



T,L 



da / d'^pid'^p2 \l/!'^^(p2,a 



(^qq[S,P2) 



dZ2 Gq|^ [Pi, Zi, P2,Z2) 



,T,L 



X ^gg- (Pl,a) (^qqis,pi] 



D 



T,L 



da / d'^r 



*9g- {ra)(^qq{r,s) 



(292) 
(293) 



where the nonperturbative qq wave functions are given by ([77| ) and ([T^). 

For a harmonic oscillator potential the Green function is known analytically. 



G^qT (P2,2i;pi, 2:2) 



a^ (a) 



27rsinh {lu Az) 



a? (a) 



exp 



e^ Az 



2v ail — a) 



(294) 



X exp 



(Pi + P2) "^oth {uj A2 



2pi • P2 
sinh iui Az) 



where Az = Z2 — Zi and 



u 



a{a) 



V ail — a) 



(295) 



is the oscillator frequency, cf. section | 

Now we have all ingredients which are necessary to calculate ( |285|) . Two from the 
eight remaining integrations, over the angles, can be performed analytically. The final 
result is 



XBj I da j dpipidp2p2 / dAz [^J^^ (e, A,p2)]* *Jg (e, A,pi) 
Jo Jo Jo 

(]? (a) 



X ^qq{p2,s) aqq{pi,s) 



X Ii 



a^ (a) P1P2 
sinh iu Az) 



exp 



sinh {u Az 
G? ia) 



exp 



e^ Az 



2 V ail — a) 
(Pi +P2) coth(cjA2;) 



(296) 
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N' 



rriN 



p1 /"OO /'OO 

xbj da dpipidp2p2 / dAz [^^^ (e, A, pa)]* * 
Jo Jo Jo 



e, A,pi) 



X 



X 



0"™(P2, 



'qq 



S) aan{pi,s) 



a (a) 



'qq 



a^ (a) pi p2 
sinh {ijj IS.Z) 



sinh [uj Az) 

2 



exp 



e^Az 



2 i/afl — a) 



D 



LT 



exp 



(ia 



a a 



(Pi + P2) coth (cu A^) 



drr 



^\t (^' ^> ^. 



99 



r, s 



(297) 



(29J 



'0 ^0 

For a dipole cross section that levels off at large separations like the older parametriza- 
tion ( |152| ) the integrations over p\ and p2 can also be done analytically. However, we 
prefer to work with the more general expressions that hold for arbitrary aqq{s,p), as 
long as it depends only on the modulus of p. The remaining integrations are performed 
numerically. The results for l'^'^{x, Q"^) are shown by solid curves in figs. ^ and |3y. 

It is instructive to compare our calculations with the VDM which is usually sup- 
posed to dominate at small Q^ < m^. The corresponding coherence length l^^^^ is 
given by (|274|) with Mqq = nip. The ratio of l^ calculated with the nonperturba- 
tive wave function to /^^^^ as function of Q^ is shown by solid curve in fig. ^. It 
demonstrates an unexpectedly strong deviation from the VDM expectation at quite 
low Q^. We also calculated /^ with the perturbative wavefunctions and a quark mass 
of ruf = 200 MeV. This choice mimics the nonperturbative effects quite well, as one 
can see from fig. RI. 



Shadowing in the nuclear gluon distributing function at small xbj which looks like 
gluon fusion GG ^ G in the infinite momentum frame of the nucleus, should be treated 
in the rest frame of the nucleus as shadowing for the Fock components of the photon 
containing gluons. 

The lowest gluonic Fock component is the \qqG) state. The coherence length rele- 
vant to shadowing depends according to (|274|) on the effective mass of the \qqG) which 
should be expected to be higher than that for the \qq) state. Correspondingly, the 
coherence length (If) should be shorter and a onset of gluon shadowing is expected to 
start at smaller Xbj- 

For this coherence length one can use the same formula as in the qq case ( |274| ), but 
with the effective mass. 



M'qqG 



Pt 



acil -ac, 



+ 



M2_ 



ac 



(299) 



where ac is the fraction of the photon momentum carried by the gluon, and Mqq is the 
effective mass of the qq pair. However, ( p99| ) is valid only in the perturbative limit. It 
is modified by the nonperturbative interaction which is much stronger for gluons than 
for quarks |jl2|. The nonperturbative interaction in the qqG state is discussed in detail 
in section T^. Therefore, we switch to the Green function formalism which recovers 
( p99|) in the limit of high Q^. 

Gluons are treated as massless and transverse. For the factor P defined in ( p74| ) 
one can write. 



P 



G 



N^ 



(300) 
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Figure 31: Q^ dependence of ratio of (Ij) calculated with (|280|) and nif = 200 MeV 



to l^^^^ calculated with ( ^74]) and M = nip. 



where 

JY^ _ _..^_ / ^2^ „.72^ ^2^_„.72^_ ^^, ^l^,'^,_^,T,t 



TTiN XBj / d Tied riqq d r2G d r2qq duq d In(aG) ^Jgc (^2G, ^2gg-, Og, og) 



dZ2GqqG{r2G,r2qq,Z2]riG,riqq,Zi)\ ^gqG (^IG, ^qq, "q, "g) (301) 



D 



G 



= / d^riG d^riqq rfVac d^r2qq duq d In(aG) ^Jg-^ (^2G, ^2qq, "g, «g) 

X 5^^^ (r2G - Tig) ^^^^ (rggg- - fi,^) ^^^^ (^g, ^gg, "q, «g) (302) 



Here we have introduced the Jacobi variables, Tqq = Rq — Rq and tg = Rg ~ {ciqRq + 
aqRq)/{aq + ttg) . RG,q,q a^G the position vectors of the gluon, the quark and the 
antiquark in the transverse plane and aG,q,q are the longitudinal momentum fractions. 
The Green function describing the propagation of the qqC system satisfies the two 
dimensional Schrodinger equation [^, 



d Q'^ aq + Uq 

dzo 



2v 2vaqaq 



mi 



A± (r2G) 
IvaG (1 - ac) 



y {f2G,r2qq,Olq,aG,Z2) 



X GqqG (^2G, r2qq, Z2\ T^G-, Aqq-, Zx) = 5 (^2 " ^l) 5^^^ (^20 " ^g) 5^"^^ {r2qq " ^qq) 



(303) 
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This equation accounts for the relative transverse motion of the three particles and for 
the interaction between them. 

In order to calculate the coherence length relevant to shadowing, we employ the 
amplitude for the diffractive dissociation 7* —* QqG, which is the qqC wavefunction 
weighted by the cross section |12], in analogy to ( ]96|) 






^S^ {rqq, aq) 



^qG [ ,rG^ Tqg 



an 



an 



a. 



m 



s,rG- 



aq 

(jIq \ (JIq 



qq 



9 



cr, 



Ola 



■Tr, 



qq \ ^1 ^G ~r ' qq 

CXq ~Y~ aq 



^qq I (^qq \S^ Tqq) 



(304) 



As different from the case of the \qq) Fock state, where perturbative QCD can be 
safely applied at high Q^, the nonperturbative effects remain important for the \qqG) 
component even for highly virtual photons. Since calculations for the full three body 
problem are obviously quite involved, one has to introduce appropriate simplifications. 
At high Q^ the size of a qq pair from a longitudinal photon is always of order ~ 1/Q- 
However, the mean quark-gluon separation at ac -C 1 depends on the strength of gluon 



interaction which is characterized in this limit by the parameter 60 ~ 0.65 GeV |[T2 



For Q"^ ^ 6q the qq is small, r^^ <^ r|., and one can treat the qqG system as a color 



octet-octet dipole, cf. section |2^ , 

GqqG {r2G, r2qq, Z2] TlG-, ^iqq, Zi) -^ G qq {r2qq, ^2\ riqq, Zi) G GG {^2G, ^2'-, ^IG, ^l) 



(305) 



Such a Green function Ggg satisfies the much simpler equation [jl2 

Ax (r2G) 



d Q' 



/,4 2 
"0 ' 2G 



GGG{r2G,Z2]riG,Zi) 



dz2 2u 2uaG (1 — Qg) ^uog (1 — Qg) 

= S{z2- zi) 5^^) (f2G - Tig) 



Correspondingly, the modified qqG wave function simplifies too, 

^qqG (^G, rqq, Og, ttc) => - "^qq (^gg, "g) ^W " V ^ gG (Tg) (TgG («, T 

where the nonperturbative quark-gluon wave function has the form [IT 



2 /aj e ■ tg / ^o 2 

^qG [tg) = hm ^gG (ttG, re) = - \ hr —^ exp -- r^ 



• G 



(306) 



(307) 



(30J 



and the color-octet dipole cross section reads. 



9 



(^GG{s,rG) = -^crgqis^rG) 



(309) 



as it was already discussed in section p^. 



3 NUCLEAR SHADOWING IN DIS AND FOR DY 88 

With the approximations given above, the factor {P'^) in ( p^74| ) for the gluon co- 
herence length is calculated in Appendix A with the result 



[p") 



31n(a^'^^/a 
6 



rmn\ 
G ) 



d6 
T 



S) 



+ 



[1 + 35] 



52 



iv*'p)-nl+i 



where 



'?/' (x) 



d\nT 



X) 



6 = 



2hl 



(310) 
(311) 

(312) 



dx ' Q2 ac 

Both the numerator and denominator in (|3 1 1|) diverge logarithmically for a^*" — ^ 0, 
as it is characteristic for radiation of vector bosons. To find an appropriate lower cut 
off, note that the mass of the qqG system is approximately given by 



M. 



qqG 



2fc2 



(313) 



where ( p99| ) was used with {p\) ^ 6q. Demanding M'^^q < 0.2s leads to a^*"" = 
26q/(0.2s — Q"^). Furthermore we work in the approximation of ac <^ 1 and also have 
to choose an upper cut off. We use 



n 



2hl 
< an < — - 
0.2s - Q2 - ^ - g2 ' 



(314) 



which means that we take only masses 2Q^ < M^^q < 0.2 s into account. The two 
limits become equal at xbj ~ 0.1. 



Our results for (P*^) = {l^) /IT°'^ ^^^ depicted in fig. ^ The approximations made 
above break down at low Q^. Calculations are therefore done only at high virtualities 
Q^ > \ GeV^. The coherence length for the qqG state is much shorter than both /J 
and l^ for \qq) fluctuations. This conclusion corresponds to delayed onset of gluon 
shadowing shifted to smaller xbj predicted in [|12] . 



3.5 The bremsstrahlungs contribution to low mass dileptons 



An advantage of the light-cone approach to the DY process |13, ffSl O, which is 



explained in section ^]^, is that the dipole cross section contains the effects of higher 
order corrections as well as nonperturbative effects. Therefore, the light-cone approach 
can be applied in the low mass region, M < 1 GeV, where pQCD is questionable. 
Although the term DY is normally only used for high mass dileptons, we will also refer 
to low mass dileptons, which are produced via the bremsstrahlungs mechanism, as DY 
dilepton since they are produced by the same mechanism. 

Low mass dileptons have been measured recently in proton nucleus {pA) ||1 70|| and 
in nucleus-nucleus {AB) ||171| , |173|j collisions at the CERN SpS collider. At a laboratory 
energy of 158 A GeV approximately three times more dileptons were detected, than one 
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(a) 





Figure 32: Diagrams for the two mechanisms for /+/^ production considered in this 
paper, (a) Dilepton production in AW colhsions which can be viewed as bremsstrahlung 
in qq scattering (hght cone approach) or as gluonic Compton scattering {Gq -^ 7*g) 
from the gluon cloud of the second nucleon (parton approach), (b) Dileptons from 
prompt gluons Gq — >■ 7*g, where the on-mass-shell gluon G (crossed by a dashed line) 
is produced in another NN collision. 



expects from the calculation of hadron decays into e^e~ in the final state. Furthermore, 
the observed shape of the spectrum is different from the calculated. Most theories 
locate the origin for the observed dileptons in the hot and dense phase of hadrons [|174| - 
185 1 , which is rather late in the time evolution of a heavy ion collision. The dilepton 
enhancement is a highly interesting effect, because it might point at a restauration of 
chiral symmetry. 

In this section, we estimate the contribution of the DY process to the dilepton 
cocktail, using the light-cone approach. These dileptons arise from the very early 
stages of the heavy ion reaction in which partons are the relevant degrees of freedom. 
We calculate (i) the direct production NN -^ l^l~X in the light cone approach and in 
a parton model and (ii), following an idea of J. Hiifner ||186|| , the lepton production via 
a gluonic Compton process GN —^ l'^l~X from prompt gluons. These prompt gluons 
are radiated in a nucleon-nucleon collision and may lead to an enhancement of the 



dilepton yield via the gluonic Compton process, see fig. p^b, where another nucleon 
scatters off the previously released gluon and radiates a virtual photon. Prompt gluons 
have been recently identified as as an important source for charmonium suppression in 
heavy ion collisions ||187|| and therefore the question arises whether they might also be 
responsible for a part of the dilepton enhancement. However, in ||186|| this mechanism 
could be excluded. 

First we discuss the direct production of lepton pairs within the light cone approach, 
fig. |32|a. We employed the older parametrization ( p.52| ) from \V2 in our calculations, 
which works very well at the low M^ values we are interested in. The improvement 
( |157|) modifies the DCS only small p, corresponding to large virtuality, and will not 
alter the result. The partonic cross section (EST) has to be embedded into the hadronic 
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process. The cross section of direct (D) dilepton production in nucleon-nucleon (AW) 
takes the form 



ttf. 






dMdy SnM 



E- 



da 



XI 



Xi 



(If 



If 



dajgN -> 7*giV) 
(i(lna) 



y 



-y 



(315) 



where y is the rapidity of the lepton pair in the cm. frame and Xi = {\/^x'jp + AM'^/s^ 
Xf)/2. The first term corresponds to radiation from the projectile quarks and the 
second term to radiation from the target quarks (see fig. |32|a), which corresponds 
to replacement y -^ —y in the first term. Since the dominant contribution to the 
cross section comes from large values of xi/a, we neglect antiquarks in the projectile, 
qj{xi/a.) = 0. We parameterize the valence quark distribution in the form 



qu[x) 



(X 



-^(l-x)^ 

\/x 



:i-x) 



c. 



X 



(316) 
(317) 



where Cu,d are defined by the normalization to the number of valence quarks iV^ ^ in a 
nucleon 



dxqu,d{x) = Nu,d- 



(318) 



In the case of the proton-nucleus (pA) and nucleus-nucleus (AB) collisions the cross 
section of the direct production is given by the integration of the nuclear densities pa,b 
over the impact parameter 



— ^ = d'^b jd^S I dZApA{s,ZA) I dzB PB{b-s, Zb) 



^2 D j2 D 



dMdy 



dMdy 



(319) 



In this expression we neglect small shadowing effects. To obtain the rate of the dilepton 
events per interaction one has to divide expression (|319|) by the total cross section, 
which can be calculated using the Glauber expressions 



CpA 



C^AB 



d% 



d% 



1-exp -a]^^TA(6) 



1 - exp -aX?v / dhTA{s) Teih - s, 



(320) 
(321) 



where cr]^^ ~ 30 mb denotes the inelastic AW cross section and we used the standard 



Woods-Saxon parameterization [|189|| for the nuclear density pa^b in 

7a,b(&) = / dzpA,B{b,z). 



(322) 



There are still two sources of uncertainty: 



The dipole cross section aqq{s,p) cannot be calculated perturbatively at large 
transverse separations, which are important in the small mass region, however. 
We employ the phenomenological expression ( 152 ). 
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Figure 33: Data ||170|| on dilepton events in p-Be and p-Au scattering together with 
predictions from different mechanisms. The various curves correspond to meson decay 
channels. The shaded areas represent our results for direct production (in the light 
cone and the parton model approaches) and via prompt gluons: The upper limits 
correspond to the quark mass ruf = 150 MeV and the lower ones to rrif = 300 MeV. 
The curves for the parton model and for prompt gluons were calculated by Yu. Ivanov. 



• The light-cone wavefunction \i/^.q(a, p) depends on the mass rrif of the quark. To 
illustrate its influence we calculate the cross sections for the two extreme cases 
of the constituent quark mass rrif = 150 MeV and ruf = 300 MeV. 

Our results for proton (p-Be, p-Au) and heavy ion (Pb-Au) scattering are shown 
in figs. |33| , |3^ under the label "light cone". The border lines of the shadowed areas 
correspond to different choices of the constituent mass ruf (150 MeV and 300 MeV for 
the upper and lower lines, respectively). We believe, that also a different dipole cross 
section cannot produce values, which are significantly out of the shaded area, because 
aqq is well constrained to describe the structure function F2 and hadronic cross sections. 

We compare these results with predictions of a parton model, which treats direct 
radiation of the lepton pair as Compton scattering of a target gluon on a beam quark 
(or vice versa) in pQCD. This formulation is in principle identical to the light-cone 
approach, but all nonperturbative and higher order effects which are parametrized in 
the DCS are not taken into account. The amplitude for the subprocess qG -^ l^l~X 
is well known and be found e.g. in [Q. The main uncertainties in this case originate 
from: 

• Possible failure of pQCD calculations for the soft Compton amplitude. 

• Poor knowledge of the gluon distribution function, especially at a soft scale. 
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Figure 34: Data [|T7l|, [171 on 
dilepton events in Pb-Au scatter- 
ing. The different curves corre- 
spond to meson decay channels. 
The shaded areas are our results 
for direct production in the light 
cone approach and via prompt 
gluons: the upper limits corre- 
spond to the quark mass nif = 
150 MeV and the lower ones to 
rrif = 300 MeV. Both processes 
contribute separately to the dilep- 
ton production and therefore their 
contributions have to be added. 
The curves for the parton model 
and for prompt gluons were calcu- 
lated by Yu. Ivanov. 



M (GeV/c^) 



We use the following gluon distribution function, 
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where Cg is defined from the condition that gluons carry half of the nucleon momentum 



dxxfcix) 



(324) 



The results for proton (p-Be, p-Au) and heavy ion (Pb-Au) scattering are shown 
4[ At low M the results of pQCD (labeled as "parton model" 



in figs. ^, ^ At low M the results of pQCD (labeled as "parton model" ) are higher 
than those from the light cone calculations but start to agree for M>1 GeV. The 
reason for this is that the saturation of the dipole cross section becomes important 
at small M, and the cross section decays weaker than 1/M^. Although the parton 
model results are supposed to coincide with the predictions of the light-cone approach 
different approximations are used and the amount of disagreement between the results 
may serve as a measure of theoretical uncertainty. Nevertheless, we believe that light- 
cone prediction is more reliable. 

The data in figs. ^ and |3^ have been measured in a pseudorapidity interval of 
2.1 < rj < 2.65, while our calculations have been performed for rapidity y = 2.4. Note 
that the data are also subject to pr cuts, which exclude very small transverse momenta, 
Pt < 50 MeV. These cuts are not included in our calculation. In the region of interest 
(M = 0.2-0.8 GeV) however both approaches underestimate the experimental data by 
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at least a factor of 10. We therefore believe, that also a more careful calculation would 
not change our results. 

The pQCD calculation of the Compton amplitude is similar to that for an interac- 
tion with on mass shell gluons. The second process, labeled "prompt gluons", considers 
gluons which are radiated in an elementary AW collision and which convert into a vir- 
tual photon via gluonic Compton scattering on another nucleon. To estimate their 
contribution to the dilepton pair production we start with the elementary Compton 
subprocess G q ^ '-/* q (see the upper part of Fig. P^ where the on-mass-shell gluon 
is marked by a dashed line). To calculate the prompt gluon spectrum we use the 



perturbative evaluation for the cross section of gluon radiation [jl5| 



d^a{q^qG) _ 3a,C 2mja^k' + [1 + (1 - aY] (F + a^mj) 



dadk"^ 



TX 
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where /c^ is the transverse momentum of the gluon, G is the factor of the dipole approx- 
imation for the cross section of a qq pair with a nucleon (in this energy range C ~ 3 
||160|| ) and a is the fraction of the quark light cone momentum carried by the gluon, as 
usual. 

Finally, we combine the cross section for the elementary process with the distribu- 
tion function qf{xi) for the quark and the distribution dnc/dxc of the gluon to the 
cross section for the prompt 



"P") process on the nucleon 
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Here, dnojdxQ is the gluon distribution over the longitudinal momentum fraction, the 
on mass shell gluon carries away from its parent quark, fig. R^. The calculations for 
the parton model and the prompt gluons were done by Yu. Ivanov and we refer to [|186|| 
for an expression for dnc/dxc and for details of the calculation. 

For nucleus-nucleus collisions the geometric factor has to take into account that 
points of gluon creation and interaction are different 
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(329) 



dzpA{s, z 
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Q^P 
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dMdy 



Again, expression ( |329| ) has to be divided by the total inelastic cross section to compare 
with experiment. 
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The results of our calculations are summarized in Figs. 531 and 133: the mecha- 



nisms considered contribute to the observed spectrum less than 10% and therefore are 
unimportant on the present level of discussion. On this level also the theoretical un- 
certainties, for instance in the choice of the value for constituent quark mass or in the 
evaluation via two mechanisms (light cone vs. parton model) are not yet of importance. 

3.6 DY shadowing in proton- nucleus collisions 



In section p.6| we introduced the light-cone approach to the DY process at low X2, 
in which DY dilepton production is viewed in the rest frame of the target and ap- 
pears as bremsstrahlung from an incident quark. In the case of a nuclear target, this 
quark undergoes multiple rescatterings and in the simplest case radiates electromag- 
netic bremsstrahlung every time it is hit. As it is well known, this naive expectation 
is wrong at high energies. Landau and Pomeranchuk [|T6i |T7|1 were the first who no- 



ticed that the cross section for bremsstrahlung of a high energetic charge in an amor- 
phous medium is suppressed compared to the Bethe-Heitler cross section. Somewhat 
later, a quantum mechanical treatment of this effect, today known as the Landau- 
Pomeranchuk-Migdal (LPM) effect, was given by Migdal [|1^. The LPM effect was 
measured recently at SLAG pO, pTI, |22|. 



At low X2, multiple rescattering of the quark inside the nucleus leads to an LPM 
suppression of bremsstrahlung from that quark. This suppression is the nuclear shad- 
owing effect for the DY process, which was measured for the first time by the E772 
collaboration p|. 

What is the physical mechanism behind this suppression? The answer is that the 
quark needs a rather long time to recreate its electromagnetic field. This coherence 
time is related to the longitudinal momentum transfer by the uncertainty relation 



lc = -= . ±:"y%, . , (330) 



ql {l-a)M^ + a'^m'l+p^' 



where Eg is the energy of the projectile quark. If the quark is hit again before it has 
recreated its field, it cannot radiate another photon. Note however that the Fourier 
modes of the field with high transverse momentum px are recreated earlier as the 
one with low px- A large transverse momentum corresponds to a small size of the 
7*g-fluctuation. The small size fluctuation are therefore less shadowed. The low pt 
components of the field correspond to large 7*g fluctuations which are strongly shad- 
owed. 

Since the photon carries the momentum fraction Xi away from the proton, the 
coherence length for DY can be written as 

1 (l-a)M2 

771^X2 (1 — a)M2 + a'^m'l + p"^ 

The role of xbj in DIS is played by X2 in DY. One can define a mean coherence length 
for DY in analogy to the one for DIS, see section |3.4|. We will however not pursue the 



concept of the mean coherence length for DY in this work. This has to be postponed 



3 NUCLEAR SHADOWING IN DIS AND FOR DY 



95 



to future studies. Like in low xbj DIS, the coherence length for DY can become 
substantially longer than a nuclear radius. 

In the case of an infinitely long coherence length Ic —^ oo, the whole target acts like 
a single scattering center and the partonic DY cross section can again be written, cf. 
(|129|), in factorized light-cone form 



din a 



d^p\'^^*^qia,p)\'^j:gq{s,ap), 



(332) 



where 



Sgq(s, ap) = 2 d b ll — exp 



cr, 



qg 



[s,ap) 



Tib) 



(333) 



is the eikonalized dipole cross section. One immediately sees that (|332|) is the analog 
of the eikonal formula for DIS ( |189| ). The eikonal term obviously leads to a reduction 
of the total cross section. It contains the LPM suppression. As already pointed out 
in section p.l| for the case of shadowing in DIS, also ( p32| ) is different from the usual 
Glauber formula, where (Jqq{s, ap) is averaged in the exponent. 

In the case of a finite coherence length, which is relevant for all present and most 
of the future experiments, the DY cross section on the partonic level can be calculated 



with the same Green function technique as the DIS cross section p5[] 
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(334) 



cf. section |3.2| . The Green function G fulfills a two dimensional Schrodinger equation, 
similar to ( |197D , 
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(335) 



The two dimensional Laplacian in ( p35| ) acts on the transverse coordinate and takes into 
account that the quark is deflected due to the radiation process by a small, but non-zero 
angle. Note, that for convenience we included the phase factor exp(— ig™"(z2 — -^i)), de- 
scribing the longitudinal motion, into the Green function G, (g£**" = iff /2Eqa (1 — a)). 
This phase factor is not contained in W and appears explicitly in ( |194| ). For a calcu- 
lation that can be compared with experimental data, one has to embed the partonic 
cross section into the hadronic environment. For a proton as projectile, this is done in 
the same way as in ( |137| ), 
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dW^dxF 67rM2 xi + X2 



XX 
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a ^ \ a J din a 



(336) 
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Figure 35: The dependence of shadowing in proton-nucleus colhsions for the DY 
process on xp for FNAL (E772) (v^ = 38.8 GeV) and RHIC (v^ = 200 GeV) 
energies. The calculations are for tungsten. Both curves are calculated for M = 5 GeV 



and show the same points as in fig. ^ At RHIC, the whole Xi?-range is shadowed. 



Given a parametrization of the proton structure function Fl", like the one from 
|124|| , one is in the position to calculate shadowing for the DY process. A comparison 



to experimental data will be presented in section |3.7| . First we investigate the behavior 
of shadowing and find out, how much it will influence dilepton production at RHIC. 

An analytical expression for the Green function G can be obtained if one employs 
the approximation agg{a, ap) = C{s)a^p'^ for the dipole cross section and a uniform 
nuclear density, cf. section p.2| . Such an approximation works remarkably, especially 
for heavy nuclei. Nevertheless, it can improved if one makes use of the fact that the 
asymptotic expression ( P32| ) is easily calculated with the realistic parametrization of 
the dipole cross section (|157|) and realistic nuclear density. One needs to use the full 
Green function only in the transition region from no-shadowing to a fully developed 



shadowing given by (|332|) . Therefore, the value of the uniform nuclear density is fixed 
by demanding that the asymptotic shadowing ratio is the same for the realistic Woods- 
Saxon form of the nuclear density ||189|| and for the effective uniform density po- We 
have checked that the found value of po is practically independent of the value of the 
cross section in the interval 1 — 50 mb. The factor C{s) is fixed in the same way, but 
now separately for transverse and longitudinal photons and for each value of a. All 
further analytical calculations are similar to the one performed for shadowing in DIS, 
see appendix |A.4| . 
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Figure 36: The dependence of shadowing in proton-nucleus colhsions for the DY 



process on X2- The curves show the same points as fig. q^. The calculations are for 
tungsten. Much lower values of X2 are reached at RHIC than at FNAL. 



Our predictions for RHIC are depicted in figs. 36 and 35 in comparison with a 
curve calculated at the highest energy that can be reached today at FNAL. ISR at 
CERN can actually reach higher energies, \/S = 62 GeV, but did never include nuclei 
nor performed measurements at low X2- The curves in the two figures show the same 
points, but in one figure, |3l 

Xp. 



displayed vs. X2 and in the other figure, ^, vs. Feynman 
Shadowing in DY is more complicated than in DIS, because there a two structure 
function, of the projectile and of the target, involved in DY. Calculations are performed 
only for values of Xi;' > 0.1, because at smaller values quarks from the target will also 
radiate in the same direction, which is not taken into account in 



One immediately recognizes that RHIC can reach much lower values of X2 than 
FNAL. Even more striking, the whole xp range at RHIC shows strong shadowing. 
This illustrates the close connection between heavy ion collisions at high energies and 
low X physics. The physics at RHIC is governed by coherence effects. The situation 
is different at SpS, where the energy is not high enough to allow for strong coherence 
effects. For a reliable interpretation of the data from RHIC and LHC it is essential to 
have good theoretical understanding of shadowing effects, because DY dilepton serve as 
a reference for the study of heavy quarkonium suppression, esp. J/ip, which is believed 
to be a signal for quark gluon plasma formation [lE^ . 

It is also interesting to investigate the influence of the nuclear medium on the 
transverse momentum distribution of the DY pairs. The formula for the differential 
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cross section was derived in |T3. It reads 
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(337) 



where kx is the transverse momentum of the quark in the final state. It is convenient to 
introduce the transverse momenta relative to the direction of the virtual photon. The 
transverse momentum of the incoming quark relative to this direction is then given by 



Pi 



Pt 



a 



and the transverse momentum of the final quark by 



P2 = kT- 



a 



a 



-Pt- 



(338) 



(339) 



The interaction with the nuclear medium is contained in the Green function G, which 
fulfills the two dimensional Schrodinger equation ( p35| ) and in the optical potential 



Voptib,p,z) 



-pA{b,z)agg{s,p). 



(340) 



This potential describes the shadowing corrections and is also present in ( P35| ). The 
operators O, acting on the variable pi or p2 in the Green function G in ( p37|) , are given 
by ( |132| ) and (|133|) . They are the same operators which appear in the definition of 
the LC wavefunctions for the transition q —>■ 7*g, cf. section p.6| . It has been shown 
in [jl5|, that in the limit of infinitely high energy of the incoming quark, Eg — > 00, the 
differential cross section ( |337|) takes an eikonal form, similar to (|146|) , 



da{qp -^ qi*p) 
d In ad^pT 



(2vr)2 
1 



d^pid^p2 exp[ipT ■ (Pi - P2)]^**q(a,Pi)^7*g(«,P2 



X 2 {^M-(«Pi) + Sgg-(ap2) - Sq5-(a(pi - P2))} 



(341) 



where S^q is given by ( 13331) . 

The transverse momentum dependence of the LPM effect, calculated from (|33 
is depicted in fig 



The calculations are for tungsten. The mass of the lepton 
pair is M = 5 GeV and Feynman xp = 0.7. Since ( |337| ) is very complicated, the 
following simplifying approximations are employed. First, a uniform nuclear density of 
Pa = 0.16 fm~^ is used, which is well justified for a heavy nucleus like tungsten, and 
second, the approximation aqq{s,p) = C{s)p'^ is applied. The latter approximation 
works remarkably well, especially for heavy nuclei [p.52|] . With these approximations. 
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Figure 37: The transverse momentum dependence of the LPM effect for different 
energies. The curves are for tungsten and M = 5 GeV, xp = 0.7. an approximate 
dipole cross section aqq{s,p) = C{s)p'^ is used. 



some of the 12 integrations in ( |337| ) can be performed analytically. The result of the 
analytical calculations is given in appendix |A.3| . One further integration has to be 
performed, in order to account for the hadronic environment, cf. ( |149| ). To make the 
calculations even more realistic, an energy dependence of the factor C{s) is introduced. 
At FNAL energies, a value of C{y/s = 38.8 GeV) = 3 is appropriate, while at RHIC 
energies one rather has C{^/s = 200 GeV) = 6. 

Note that pairs with low transverse momenta are suppressed at small Pt < 1 GeV, 
but enhanced at p^ > 1 GeV. This enhancement in the intermediate pt region is already 
known for a long time as Cronin-effect [ |190| |. It is usually assumed that the Cronin- 
effect is caused by multiple scatterings and the light-cone approach reproduces this 
phenomenon. At very large Pt, all nuclear effects vanish. Furthermore, the shadowing 
at small pt increases with energy, since the dipole cros section becomes larger. 

Although ( p371) looks rather complicated, it contains intuitively understandable 
physics. The nucleus acts as a color filter [Q. Small size dipoles, having large px, can 
pass the nucleus nearly undisturbed while dileptons with small transverse momentum 
are shadowed, because they correspond to large arguments p of the dipole cross section. 



Note that ( p371 ) contains also the phenomenon of nuclear broadening. Due to multiple 
rescatterings, the projectile quark performs a random walk in impact parameter space, 
while propagating through the nucleus. Therefore, not the whole suppression at small 
Pt is due to nuclear shadowing. Part of the dileptons reappear at larger transverse mo- 
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mentum, which leads to the enhancement in the intermediate pt region. Indeed, it was 
found in ||191|] , that in the approximation aqq{s,p) = Cp^, the transverse momentum 
broadening is given by 

HptV^ = {Prr^ - {Prr = C{T{b)), (342) 

where (T(6)) = j d%T'^{h)/A, is the mean nuclear thickness. 

First numerical calculation with the eikonal approximation ( |341|) are performed in 



|T5|, but only for quark- nucleus scattering. Our results are consistent with the one from 



15|. Note that the transverse momentum dependence of the LPM effect in QED is also 



investigated in [^. The case of gluon radiation off a quark is studied in [|192| , |193 



Recently, also the approach by Luo, Qiu and Sterman |68, ^] has received much 



attention, see e. g. [ |194|| - [T97| . This approach extends the QCD factorization theorem 



198|| to higher twist operators. The double scattering contribution to the DY process is 
described in terms of a matrix element, which contains all the nonperturbative physics, 
and a hard part that can be calculated in perturbation theory. However, it was already 
mentioned in section |2.5| that the calculation of the transverse momentum distribution 
is very complicated in pQCD, since it requires the resummation of large logarithms 
ln(p|,/M^) at p\ <^ M^. It is known, how to perform this resummation for the leading 
twist ||120| , |121| , |122|| , but it still has to be investigated, how these techniques can be 



generalized to higher twists. 



3.7 Comparison with data 



The mean coherence length introduced in section ^]J alone is not sufficient to predict 
any nuclear effect in the structure function F2. It is merely intended to be a tool for 
qualitative considerations. The Green function technique from section p^ does not need 
the mean coherence length as input. Instead, all coherence effects are automatically 
taken into account by the Green function ( p.97| ). Together with a realistic phenomeno- 
logical dipole cross section we are now able to perform parameter free calculations for 
nuclear shadowing in the region which is most difficult for theory: The transition region 
between no shadowing at Xbj ~ 0.1 and saturated (for the \qq) component) shadowing 
at very small Xbj- 

In the rest frame of the nucleus shadowing in the total virtual photoabsorption 
cross section aj^^ (or in the structure function F^) can be decomposed over different 
Fock components of the photon, 

at/ = Aafj - Aa,,t{qq) - Aa^qqC) - Aatot{qq2G) - ... (343) 



It was found in section |3]J that the coherence length for the \qqG) fluctuation is rather 
short at xbj ~ 0.01, compared to the mean internucleon spacing in a nucleus, because 
of it's large mass, see fig. 29. Although the calculations were done only at large Q^ 
and some data are at quite small Q^, we neglect the \qqG) component of the virtual 
photon for the present comparison with data. It is one of the most important tasks for 
the future, to include also the gluon, since it is important for shadowing of the nuclear 
gluon distribution. The treatment of the propagation of the three body system qqG 
through the nucleus is however numerically quite involved. 
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Holding only the first two terms on the r.h.s. of ( ^43| ) , the total cross section for 
transverse and longitudinal photons on a nucleus can be represented as 
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(344) 



where p^(&, z) is the nuclear density dependent on impact parameter b and longitudinal 
coordinate z. In (E 



|) we have generalized our result from section p.2| to include also the 

nonperturbative interaction proposed in |jl2| , see also section |2]^. The nonperturbative 
wave functions for the qq component of the photon are defined in (|77D-([78D. The Green 
function G{p2,Z2 \pi,zi) describes propagation of a nonperturbatively interacting qq 
pair in an absorptive medium. It fulfills the two dimensional Schrodinger-equation, 



d , A^{p2)-e' , i 
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a^{a)pl ^ 



2z/a (1 — a) 
= iS {z2 - 



G{p2,Z2\pl,Zi) 

zi) S^'Hp2 - Pi) , (345) 



where a{a) and A are given by (|76|) and ( |8TD respectively. The third term on the l.h.s. 
of (|345|) describes the absorption of a qq pair with cross section agg (s, r) in the medium 
of density pa {b, z) . Note that the Green function G includes now also the motion in 
longitudinal direction, in contrast to the Green function W in ( |197|) . The longitudinal 
motion is contained in the term e^ jlva (1 — a) in ( |345| ) . 

An explicit analytical expression for W can be found only for a dipole cross section 
Ogg^s, p) = C{s)p'^ and a constant nuclear density pa(&, z) = Pq. Such an approximation 
has a reasonable accuracy, especially for heavy nuclei. Nevertheless, it can be even 
more precise if one makes use of the fact that the asymptotic expression (|189|) is easily 
calculated with the realistic parametrization of the dipole cross section ( |157| ) and a 
realistic nuclear density. One needs the full Green function only in the transition 
region from no-shadowing to a fully developed shadowing given by (in 



Therefore, 

we fix the value of the uniform nuclear density by demanding that the asymptotic 
shadowing ratio is the same for the realistic Woods-Saxon form of the nuclear density 
|189|| and for the effective uniform density po. We checked that the found value of po is 
practically independent of the value of the cross section in the interval 1 — 50 mb. The 
factor C is fixed in the same way, but now separately for transverse and longitudinal 



photons and for each value of a. More details can be found in appendix [A. 4 
We perform calculations with the parameter v in 



fixed at w = 0.5, because the 
parameters in the dipole cross section are fixed for this value of f , see section |2.7| . The 
dependence on v is very weak. All calculations are also performed with perturbative 
wavefunctions (^, ^) and the quark masses given on page |^. The results are shown in 

figs, ps^ 



Curves obtained with nonperturbative wavefunctions are drawn with solid 
lines and the results from perturbative wavefunctions are depicted by dashes lines. The 
calculations are for the same kinematics as the datapoints. We neglect shadowing in 
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the deuterium structure function F.p, which is of order < 5% ||199||. The antishadowing 



effect which occurs at xbj ~ 0.1 is neglected as welL Shadowing is calculated only for 
nuclei with A > 12, since the basic formula is derived in the approximation A ^ 1. 
Although the calculations are parameter free, agreement is reasonably good. 



Fig. 38 shows a comparison of the calculations with experimental data for the 



structure function ratios of carbon and calcium over deuterium from NMC ||200|, 1201 



Unlike other NMC data, these data reach down to quite small xbj , where the error bars 
for carbon are quite large. At the last point in the carbon plot, the solid curve turns 
unexpectedly upwards. No such behavior is visible in the perturbative calculation. 
Note, that the data are at different values of Q^. The leftmost point has Q^ = 0.035 
GeV^ and the rightmost point Q^ = 9.6 GeV^. Because of the very small Q^ of the 
point at the lowest value of xbj, we attribute the decrease of shadowing to the model 
of the nonperturbative interaction. The calcium data and part of the carbon data 



shown in this figure are the result of a reevaluation of the data from ||202| , p03|| . Among 
other reasons, the reevaluation was necessary, because the radiative corrections require 
knowledge of F:P over a large range of xbj and Q^. When new information about 
F^ became available from SLAC, the radiative corrections were performed again ||200|| . 
This lead to an increase of the structure function ratio by at most 2.5% for Ca/D at 
the lowest value of xbj [pO0|] . 



In fig. |39|, we compare our calculation to NMC data from ||161| , pOO|| for different 
nuclei, ranging from aluminium to lead measured relative to carbon. Shadowing is al- 
ways slightly larger in the nonperturbative calculation. This is in accordance with the 
observation in section |2.7| that the choice of the quark mass squeezes the wavepacket 
somewhat more than the nonperturbative interaction. Since shadowing is determined 
by diffraction, which is an almost entirely soft process, one is more sensitive to the large 
p behavior of the dipole cross section and to details of the nonperturbative interaction. 
The perturbative calculation seems to reproduce the data better than the nonpertur- 
bative. Note however that taking into account the antishadowing effect would shift the 
curves upwards. Whether this will happen at all values of Xbj or only around Xbj ~ 0.1 
is model dependent. We prefer not to introduce new parameters. Note that the data 
in fig. ^ correspond to different muon beam energies, 90 GeV for the calcium data in 
the upper row and 200 GeV for all other data. The structure function ratio is however 



independent of the beam energy. Indeed one sees from the equations in section BTT] that 



pZ? 7*1) I 7*D' 

while for the measured lepton-target cross section one has 



A ^■y'A I ^■y'A 

(346) 



F,^ _al^ + a. 



a^A F^ {1 + eRA)il + R_ 



D) 



where 



4(1 -y)-^^ 
4(1 -y) + 2y2 + '-^ 



(347) 



(348) 
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is the photon polarization parameter and 
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The lepton nucleon cross section ratio is only equal to the structure function ratio, if 
there are no nuclear effects in R. This is assumed in all NMC data and supported by 
experimental observations. In the kinematical region relevant for NMC, R data show 
practically no A dependence [p04| . 



The best data available today for nuclear shadowing in DIS are for the structure 
function ratio Sn/C from NMC ||205|| . These are shown in figs. [4^ and ^ Although 
our calculation describes the existing data quite well, the dip in the Q^ dependence of 
the solid curve in fig. |T looks rather suspicious. This dip is practically invisible if one 
employs the older parametrization (|152|) of the dipole cross section ||206|| . It is therefore 
caused by the suppression factor in the improved parametrization ( |157| ). Since the dip 
is not present in the perturbative calculation, it is also related to the nonperturbative 
interaction. We do not believe that this indicates a physical effect. One rather has to 
improve the dipole cross section and the nonperturbative interaction, in order to get a 
better description of the Q^ dependence of shadowing. We also hope to overcome the 
numerical prescription we use for the dipole cross section in the future and to perform 
the calculations with the exact parametrization. 

We also compare our parameter free calculations to data from E665 |p07| , |208|| , where 
much lower values of xbj are covered, fig. ^. The decrease of shadowing in the solid 
curve at very low xbj and Q^, which was already visible in fig. |3^, is also seen in the 



Xe/D plot. The agreement with data is less good for the E665 data than for the NMC 
data. Indeed, the two datasets contradict each other ||161| , p04|| and no calculation can 
reproduce both. The disagreement between NMC and E665 vanishes however, in the F2 
ratio relative to carbon. Note also that the qqG component will become important at 
very low Xbj and calculations, taking only the qq state into account, will be unreliable 
at xbj ^ 0.01. 

Besides describing all data reasonably well, the success of the Green function tech- 
nique is to reproduce the correct shape of the curves vs. xbj- Without treating the 
coherence length properly, shadowing would set in quite abruptly at Xbj ~ 0.1 and 
quickly saturate as one goes to lower xbj-, cf. fig. ^. 

Nuclear shadowing is very challenging from the experimental point of view, because 
of large radiative corrections. These corrections occur, e.g. when the incident lepton, see 
fig. ^, radiates bremsstrahlung and does not produce a deep inelastic event [^ . In the 



older publications ||202| , |203|| , the NM collaboration calculated the radiative corrections 
employing the computer program FERRAD, which relies on the theoretical analysis 



For the reevaluation of the shadowing data [ 200 | three different 
FERRAD, an improved version of FERRAD and TERAD. The 



by Mo and Tsai [gl] 

codes were tested 

latter relies on the work of Akhundov, Bardin, and Shumeiko pOQl |210|, pil|] and was 



used for the published NMC result. A more detailed discussion on radiative corrections 
for NMC and HERMES can be found in ||212|| . We point out that the whole shadowing 



effect in the data shown in figs. 38-41 is calculated. Without radiative corrections, the 
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Figure 38: Comparison between calculations for shadowing in DIS and experimental 
data from NMC [pOOj , |201|| for the structure functions of carbon and calcium relative 
to deuterium. The x^j-dependence of shadowing is shown. The leftmost point of the 
carbon data is measured at Q^ = 0.035 GeV^ and the rightmost point at Q^ = 9.6 
GeV^. The Q^ values for calcium are 0.6 GeV^ < Q^ < 9.7 GeV^. The inner error 
bars show the statistical error and the outer error bars statistical and systematic error 
added in quadrature. The full curves are calculated with the Green function method 
( P44|) including the nonperturbative interaction. The dashed curves are calculated with 
the perturbative wavefunctions 



cross section ratio would be around unity ||212| , |215|| . However, the correctness of the 



calculation was checked by comparing to the "hadron tagged" data, making sure that 
really a deep inelastic event was measured. 

A different way to identify the deep inelastic events was chosen by the E665 col- 
laboration ||207| , |208|| , where cuts in the electromagnetic calorimeter were applied. In 
| 207| the xenon data obtained in this way were compared to an evaluation with hadron 
tagging, see fig. 42. The two methods give consistent results. In [208|, the E665 col- 



laboration also applied the FERRAD code for a comparison with NMC. The results 



are depicted by open circles in fig. |4^. The code wass assumed to be reliable only in 
the region y < 0.7, see (|]), and xbj > 0.002 ||208|| . Therefore it was not applied to 



all data points, but only to those, where it was supposed to work. One recognizes a 
systematic discrepancy between the two evaluation methods. The radiative corrected 
cross section ratios are lower by ~5% than the points from the calorimeter analysis. We 
argue however from the date of the publication that the improved input for FERRAD 
used in |^0U|] was not taken into account by E665 in ||208||. It is of course impossible 



for the author to conclude which experiment is correct. 
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Figure 39: Comparison between calculations for shadowing in DIS and experimental 
data from NMC ||161| , |200|| for the structure functions of different nuclei relative to 
carbon. The aj^j-dependence of shadowing is shown. The Q^ range covered by the 
data is approximately 3 GeV^ < Q^ < 17 GeV^ from the lowest to the highest xbj bin. 



The different curves and error bars have the same meaning as in fig. 38 
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Figure 40: The xbj dependence of nuclear shadowing in DIS for the structure function 



ratio of tin relative to carbon. The data are from NMC ||205|| . The different curves and 



error bars have the same meaning as in fig. ^. 
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Figure 41: The Q^ dependence of nuclear shadowing in DIS for the structure function 
ratio of tin relative to carbon. The data points [POH] are the same as shown in fig. EC. 



The different curves and error bars have the same meaning as in fig. 38. 
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Figure 42: The Xbj dependence of nuclear shadowing in DIS for various nuclei. The 
data are from E665 [p07 , 



niql - Full circles show data taken with electromagnetic 
calorimeter cuts, while the open circles result from the FERRAD radiative correction 
code. Hadronic cuts were applied for the triangles in the lower left plot. For better 



visibility, triangles and open circles are displaced to slightly lower xbj- The Q range 
of the data is 0.15 GeV^ < Q^ < 22.5 GeV^ except for the xenon data, 0.03 GeV^ < 
Q^ < 17.9 GeV^, from the lowest to the highest value of xbj- The meaning of the 



different curves and error bars is explained below fig. ^. 
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Figure 43: The left figure sfiows tfie leptonic cross section ratio measured by HERMES 
for helium and nitrogen. The data are not at variance with other experiments, since 
the points correspond to different values of Q^. The large enhancement in the i?-ratio 
is shown on the right, Ra = ^1 /cr'^ The figures are from ||213 . 



Much attention has been drawn to a possible A dependence of R and the importance 
of radiative corrections by the data released recently by the HERMES collaboration 
| 213|| . HERMES measured nuclear shadowing in DIS at the same values oixsj as NMC, 
but at smaller virtuality, 0.5 GeV^ < Q^ < 2.018 GeV^ and larger y < 0.85 (|). The 
27.5 GeV positron beam at HERA was used. The results shown in fig. |4^, expose several 
striking features. The cross section ratio for nitrogen to deuterium is suppressed much 
stronger than it would have been expected extrapolating previous NMC measurements. 
This ratio is related to the F2 ratio by ( p47| ). The two quantities are equal, only when 
e = 1 or when there are no nuclear effects in R (|349|) . Since NMC could not find any 
nuclear effects in R, all NMC data were evaluated assuming Ra = Rd- Note that 
at HERMES e ~ 0.4 at low xbj- The discrepancy in the left part of fig. ^, can be 
interpreted in terms of different shadowing for transverse and longitudinal photons. 
It was concluded in ||213|| that while aL is enhanced, ay is suppressed on nitrogen by 
at least factor of two compared to a deuterium target. Then one obtains a structure 



function ratio of the same order of magnitude as observed before, fig. P 

Radiative corrections are huge at HERMES, because the positron mass is much 
smaller than the muon mass and because radiative corrections become very large at y ^ 
1 1^. Furthermore, F2 and R are needed as input to calculate the radiative corrections 
and have to be known over a wide range oixBj and Q^. Therefore an iterative procedure 
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Figure 44: The F2 ratio for nitrogen to deuterium measured by HERMES [213 



com- 



pared calculations with the Green function method (|344|) . The full curve is calculated 
with nonperturbative and the dashed curves with perturbative wavefunctions. The Q^ 
range of the data is 0.5 GeV^ < Q^ < 2.018 GeVl 



was applied by HERMES pi2| , |213| |, using existing data as input to calculate the 
radiatively corrected cross sections in a first approximation. The new structure function 
and the new R obtained in this way were again used as input to the computer program. 
This was repeated until convergence is obtained. To make the result reliable, three 
different codes were compared, FERRAD [|T|], TERAD [|0|, |T0], |2lTl and POLRAD 
| 216| |. In spite of the statement in | |208| | that FERRAD is not applicable for y > 0.7, 
all codes gave the same result within 2%. More details can be found in ||212 . 



The HERMES experiment drew attention to the fact that only very few data are 
available in this kinematical region. Moreover, no reliable theoretical calculations are 
done yet. Approaches based on QCD evolution equations need parton distributions 
at a semi hard scale as input and their applicability is doubtful at the low Q^ values 
of HERMES. More promising is the intuitive approach treating nuclear effects in the 
spirit of vector dominance model (VDM) ||147|| as shadowing for the total cross section 
of hadronic fluctuations of the virtual photon (see e. g. ||140|| ). However, the VDM is 
sensible only at small Q^ ^ not so easily generalized to longitudinal photons. 

Eventually we are able to provide theoretical predictions from the light-cone ap- 
proach, since nonperturbative effects are assumed to be parametrized in the dipole 
cross section and we also included the nonperturbative interaction introduced in |1^ 



between the q and the q for a more realistic description of the soft physics. The results 
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Figure 45: The result of calculations in the light-cone approach compared to the 



HERMES measurement [pl3|] of Rn/Rd- The data points are for the same kinematics 
as those shown in fig. 0. On this scale, the lines corresponding to perturbative (dashed) 
and nonperturbative (solid) wavefunctions lie on top of each other and predict an R- 
ratio of practically unity. The large enhancement for o"2 measured by HERMES is 
clearly not reproduced in the light-cone approach. 



of our calculations for the nitrogen F2-ratio is shown in fig. ^ in comparison with 

HERMES data. We expect somewhat larger shadowing for F2 than seen by HERMES. 

The most striking feature of the HERMES data for shadowing is a dramatically 

rising R = a^/aT ratio on nitrogen compared to proton target at Q^ < 1 GeV^ 



213|| . Our predictions for R^/Rp are plotted in fig. ^ versus Q^ at the same values 
of xbj as the experimental points. Like in all other calculations, we neglect nuclear 
effects in deuterium. Apparently, we do not expect any remarkable effect. There is 
a priori no reason to think that the light-cone approach leads to equal shadowing for 
longitudinal and transverse photons. It is calculated in section |3.4| that the coherence 
length for longitudinal photons is significantly longer than for transverse. Apparently 
this is compensated by the smaller size of gg-fluctuations from longitudinal photons, cf. 
section p.3| . One has however to bear in mind, that we can calculate shadowing for a^ 
only at small Q^. In the approximation, when only the qq pair is taken into account, 
shadowing for longitudinal photons vanishes oc l/Q"^ at large Q^. Taking also the qqG 
component into account, shadowing for ai will persist at high Q^, cf. section |2.4|. We 



also calculated R for a proton target and compared to NMC data [Q in the low xbj 
region, fig. ^ Good agreement with the data is achieved over the whole Q^ range. 
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The inner error bars represent the statistical 
error and the outer error bars statistical and systematic error added in quadrature. 
The full curve is calculated using nonperturbative and the dashed curves using per- 
turbative wavefunctions. From the left to the right, the data are for the Xbj values 
0.0080, 0.0045, 0.0125, 0.0175, 0.025, 0.035, 0.05, 0.07, and 0.09. 



without introducing any new parameters. One might argue that the beginning of an 
increase of R is visible from the low Q^ points in fig. |^. There are however not enough 
data to allow for definite conclusions. Note that nuclear effects for R were estimated 
previously in a different approach 



We suppose that our results provide a reliable base line for nuclear effects in this 
region. The dramatic effects revealed by the HERMES experiment probably cannot 

[HI 



While a 



be explained without involving new nonstandard dynamics, see e.g. 

nuclear enhancement in the longitudinal cross section might be qualitatively explained 

by higher twist contributions, the enormous suppression for O'y which is implied by 

the o"L enhancement and fig. ^ still remains. Indeed, for the two leftmost points in 

fig- 



one finds 



g a^"" lair" 



0.45. For nitrogen, this implies that the virtual 
photon splits up into something that interacts with an effective cross section of ~ 100 
mb, compared to the usual ~ 25 mb for real photons and for Q^ > \ GeV^. 

Experimental data for shadowing in DY are available from the collaborations E772 
and E866 |^17 . We compare only to the E772 data, although even smaller values of 
X2 could be reached by the E866/NuSea collaboration. The latter measured however 
relative to a beryllium target. The large A approximation and the uniform nuclear 
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Figure 47: Comparison between calculations in the light-cone approach and E772 data 
P] for shadowing in DY. Only perturbative wavefunctions are applied, since there is 
no real dipole in DY. The calculations are performed at the mean values of the lepton 
pair mass. From the left to the right, these values are 5 GeV^, 5.7 GeV^, and 6.5 GeV^. 



density we have to apply is not justified for such a light nucleus. These simplifications 
are necessary, because of merely technical difficulties and are not demanded by the 



nature of the light-cone approach, section |3.6| . The result of the calculation, shown in 
fig. 0, agrees well with the data. 

Practically all DY data are taken at masses M higher than the J/tp mass, because 
of the large underground of dileptons from other sources at small M^. It is therefore 



impossible to reach values of X2 as low as the values of xbj in the E665 data, fig. |4^, at 
present. This situation will however change in the nearest future, when RHIC starts 
operating. With PHENIX, values of Xa ~ 0.001 and M ~ 2 GeV can be reached |218|. 
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4 Summary and outlook 

In this work a phenomenological approach to nuclear shadowing in DIS and DY in 
the target rest frame is developed. As it is well known, the partonic interpretation of 
high energy scattering processes depends on the reference frame. In a frame, where the 
target is fast moving, DIS looks like the scattering of virtual photons off sea quarks 
in the target (at low xbj)- In contrast to this, in the rest frame of the target the 
virtual photon splits up into a gg-pair which then scatters off the target gluon field, 
see fig. |l|. Note that in the infinite momentum frame the gg-pair belongs to the target 
and is generated from gluon splitting, while in the target rest frame the incoming 
virtual photon is decomposed into a superposition of Fock-states. The gg-pair has to 
be regarded as the lowest Fock-state of the photon. 

The DY process may be seen in a way similar to DIS: In the target rest frame 
DY dilepton production looks like bremsstrahlung [|T3[, see fig. 0. A quark from the 



projectile scatters off the target gluon field and emits a photon of mass M, which 
decays into a dilepton pair. Like in DIS, the incoming quark is decomposed into a 
superposition of Fock-states, but now the lowest nontrivial state is the g7*-fiuctuation. 
The single quark and the g7*-state scatter with different amplitudes off the target. 
This disturbes the coherence between the Fock-states and the 7* is freed. 

It is a special virtue of the target rest frame point of view that coherence effects, 
which are so important at high energies, can be understood in a very easy and intuitive 
way. 

It is convenient to work in a mixed representation, where the two transverse di- 
mensions are treated in coordinate space while the longitudinal direction is described 
in momentum representation. The advantage of this representation is that Fock-states 
with fixed transverse separations between the particles are eigenstates of the interac- 
tion at infinitely high energy P, |^, §. The total cross section for DIS and DY can then 
be written in factorized light-cone form, namely as product of a LC wavefunction and 
the dipole cross section, integrated over the two transverse directions and the longi- 
tudinal momentum. The differential DY cross section can also be expressed in terms 
of the dipole cross section |T^. The LC wavefunctions describe the probabilities for 
the transitions 7* — > gg (DIS) and g -^ 7*g (DY), respectively. In the case of DIS, 
a nonperturbative interaction [|12| between the particles in the pair is introduced to 
model confinement. 

Working in a representation, in which the scattering matrix is diagonal, makes it 
particularly simple to resum multiple scattering terms. In both cases, DIS and DY, 
shadowing is caused by multiple scattering within the coherence length. The coherence 
length is the lifetime of the fiuctuation of the projectile and can substantially exceed the 
nuclear radius at high energies. In DIS, large size dipoles are absorbed at the surface of 
the nucleus, casting a shadow on the inner nucleons. For the DY process, shadowing is 
due to the LPM suppression of bremsstrahlung, as explained in the introduction. The 
conditions which have to be fulfilled in order to observe shadowing are 

• The fiuctuation of the projectile has to interact with a cross section that is large 
enough to make the mean free path of the fiuctuation shorter than the nuclear 
radius. Only for these fiuctuations, multiple scattering is possible. 
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• The coherence length has to be longer than the mean internucleon spacing. Oth- 
erwise, the projectile would scatter incoherently off the different nucleons and no 
shadowing is observed. 

Since shadowing is dominated by large size ffuctuations, it is an essentially nonper- 
turbative phenomenon, even at high virtualities. Shadowing does not vanish at high 
virtuality, because the large cross section of the aligned jet configurations compensates 
for the small probability to create a very asymmetric pair. Note that there are no large 
gg-pairs from longitudinal photons. Nevertheless, shadowing for longitudinal photons 
does not vanish with Q^ either, when higher Fock- component of the photon, containing 
gluons, are taken into account. Since longitudinal photons measure the gluon density 
in the target, one expects that gluon shadowing is a leading twist effect. 

The most important nonperturbative input to all equations in this work is the 
dipole cross section which is essentially unknown and has to be parametrized. The 



parametrization presented in this work is an improvement of the one given in []T2| . The 
dipole cross section is known to rise quadratically with the transverse separation p at 
small p. It is expected to become constant as function of p at large separations. The 
energy dependence of the dipole cross section correlates with p. At small separations, 
the dipole cross section grows with a hard pomeron intercept while at large separations 
it grows much more slowly with a soft pomeron intercept. The dipole cross section is 
written as function of the cm. energy squared s of the colliding objects, rather than 
as function of xbj, in this work. It can therefore also be applied to calculate total 
hadronic cross sections. 

The propagation of the gg-pair through the nucleus is described in the light-cone 
approach by a Green function which fulfills a two dimensional Schrodinger equation. 
The kinetic part of the Hamiltonian accounts for the transverse motion of the particles 
in the pair which is important at finite energies. The dipole cross section multiplied 
with the nuclear density enters as imaginary potential in the Schrodinger equation. The 
optical potential accounts for the shadowing corrections. The Green function contains 
all multiple scattering terms and treats the nuclear form factor correctly. Furthermore, 
the nonperturbative interaction from [|T^ is also implemented into the formalism. It 
leads to a real part of the potential that describes the confining interaction between 
quark and antiquark. The total nuclear cross section is represented as the difference 
of the unshadowed cross section and a suppression term that represents the shadowing 
correction. This interference part contains also the longitudinal momentum transfer. 
Several known approximations can be obtained as limiting cases. For infinitely high 
energy, an eikonal formula similar to the Glauber formula [Q is recovered. For light 
nuclei, where only the double scattering term is important, one reproduces the formula 
of Karmanov and Kondratyuk [p.58|| . Shadowing vanishes as xbj — »■ 0.1, because the 



longitudinal momentum transfer leads to an oscillating phase factor which makes the 
shadowing term vanish. The same formalism can also be applied to shadowing in DY. 
The two dimensional Laplacian in the Schrodinger equation describes the deflection of 
the quark after radiation of the photon. This term becomes negligible at infinitely high 
energies and one obtains a similar eikonal expression as in DIS. The Green function 
method is however less intuitive in this case since there is no physical dipole in the DY 
process. For the same reason, the nonperturbative interaction is not applied to DY. 



4 SUMMARY AND OUTLOOK 116 



The main achievements of this work can be summarized as follows: 

• The Glauber-Gribov theory of multiple scattering is improved to include the 
correct nuclear form factor in all higher order scattering terms. Calculations for 
DIS and DY are in good agreement with data, without introducing additional 
free parameters. 

• The mean coherence length of a Fock state is defined. This quantity is intended 
to be a tool for qualitative considerations. We estimate that shadowing for the 
qqG Fock state of the photon is negligible at Xbj > 0.01. 

• For proton-nucleus collisions at RHIC energies, considerable shadowing of DY 
dileptons is predicted for the whole xj^-range. The transverse momentum de- 
pendence of the LPM effect is studied, too. While shadowing is predicted for 
dileptons with low transverse momenta, an enhancement at larger transverse 
momentum is expected. 

• The contribution of low mass dileptons, which are produced via the brems- 
strahlungs-mechanism, to the dilepton cocktail is estimated to be of order 10% 
only. 

• The transverse momentum distribution of DY pairs in proton-proton collisions is 
calculated. The result is finite at zero transverse momentum due to the saturation 
of the dipole cross section. Note that first order pQCD leads to a divergent result. 

• A new parametrization of the dipole cross section is presented. It is written as a 
function of the cm. energy squared s and the dipole size p. The parametrization 
describes well all HERA data up to Q^ < 35 GeV^, as well as DY data at large 
xp and hadronic cross sections. No K factor is introduced to describe the DY 
data. 

In spite of the successes of the light-cone approach, there are also some open issues 
which have to be considered in the future: 

• The transverse momentum spectrum of DY pairs in proton-proton collisions, 
calculated in the light-cone approach, looks different from what one would expect 
from data at higher X2- It will be interesting to compare to data at low X2, once 
they are available. 

• The strong enhancement of the longitudinal cross section o"2 , recently observed 
by HERMES in DIS off nitrogen, is not reproduced in in the light-cone approach. 

• There are still unsolved technical problems, how to do calculations with the Green 
function method, without the approximation of a uniform nuclear density and a 
quadratically rising dipole cross section. 

This work is far from being complete. There are many interesting applications and 
challenges which have to be investigated in the future. The two most natural continu- 
ations of this work are, 
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• to develop a framework that consistently takes into account higher Fock states of 
the projectile, containing gluons. This is especially important in view of nuclear 
gluon shadowing. First calculations are already done in [|l^] . 

• to investigate the suppression of gluon radiation in nucleus-nucleus collisions and 
it's influence on the physics at RHIC and LHC 

• to extend the Green function method to vector-meson production and to study 
the influence of the nuclear medium on the total cross section and on the trans- 
verse momentum distribution. 
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A Appendix 

A.l Calculations for DY cross sections in pp scattering 

In section |2.6| , the differential DY cross section is expressed as a four-fold Fourier 
integral (|146|) 



da{qp -^ (lYv) 



dlnad'^PT (27r) 



/ d'^pid'^p2 eXp[ipT ■ (pi - P2)]***q(a,Pl)*7*9(«,P2) 



X - {agg{api) + (Tgg{ap2) " ^^^(^(pi " ^ 2 ) ) } , (350) 

where 

vl>;T,(a,Pi)^r.("'P2) = ^|mXKo(r/Pi)Ko(r/p2) 

+ [l + {l-af]v'^^^^K,{7^p,)K,{r^pA, (351) 

P1P2 J 

vl/;{'^(a,pi)vI/:^.^(a,P2) = ^M^ (1 - a)^ Kq (r/pi) Ko (77P2) . (352) 

The Fourier-integral is inconvenient for numerical calculations, but one can perform 
three of the integrations analytically for arbitrary (jqq{ap). 
Consider the Kp-part first. With help of the relation 

Ko(r/p) = ^ [d'lj^-^, (353) 

one finds 

da{qp —>■ q'j*p) 



d In ad'^pT 



Kq— part 



-^ [rrifa + 2M 1 - a -—^ / d pid p2 



X 



giPT-(pi-P2) g-i^l-pi gi'2-P2 

~(^F+WTf+7r 



X - {(Tqqiapi) + (Xqq{ap2) - (rqq{a{pi - Ps))} . (354) 

Note that the term in the curly brackets consists of three contributions, which depend 
either only on pi or on p2 or on the difference pi — p2- Thus, the integral (|354|) 
can be split into three terms. In the integral which arises from the (Jqq{api) part, 
the p2-integration is trivially performed and leads to a two dimensional delta-function 
5^'^\pT — l2)- This makes it possible two perform also the I2 integration. The integration 
over /i gives just the MacDonald function Kq ( |353| ). Thus one is left with a two fold 
integration over pi. Provided the DCS depends only on the modulus of p, one can use 
the relation 

h = ^ I dct>^''' (355) 
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to perform one more integration. Here, Jo is a Bessel function of first kind. The 
contribution arising from (Jqq{ap2) is calculated in exactly the same way. For the 
o"gg(a(pi — p2))-part one has to introduce the auxiliary variable d = pi — p2, before the 
procedure described above can be applied. 

The Ki part is calculated in a similar way. Note that 



KiM 



1 d 
7] dp 



Ko(r7p). 



(356) 



The Ki-part reads 

da{qp -^ qYv) 



d In ad^pT 



a. 



Ki— part 



Tl" 



[l + (l-«)'](^/^VrfV2 



X 



giPT-(pl-p2) g-ih-pi gH2-p2 



{II + r^^) {II + r/2 



(27r)2 
k ■ h 



d% d% 

271 27r 



X 2 {'^^^("^i) + ^g?("P2) - (yqM{pi - P2))} • (357) 



Like the Ko-part, the complete integral (|357| ) is split into three pieces, corresponding 
to the three terms in the curly brackets. Again, one integration over p is immediately 
performed, leading to a 5-functions, which allows to do one integration over /. With 
the second /-integration, one recovers the MacDonald function Ki via ( p56|) . For the 
Ki-part, one also needs the relation 



dz 



(358) 



Although the calculation is slightly more cumbersome for the aqq{a{pi — p2))-part, 
which leads to an additional Ko-term in the final result, all calculations are easily 
performed. 



Finally, one finds 

dcr{qp — > qYp) 
d In ad^pT 



Up 



-n" 



+ [i + (i-«r] 



[m)a^ + 2M^ (1 
2r]pT 



a] 



_p\ + 772 



Xi 



4q' 



-Xo 



2 , 2-^3 — Xi + -X2 
pf^ + 7]^ 2 



with 



Xi 
X2 
X3 



drrJo{pTr)Ko{r]r)agg{ar) 
drr'^3o{pTr)Ki{r]r)aqq{ar) 
drrJi{pTr)Ki {rjr)aqq{ar). 



(359) 

(360) 
(361) 
(362) 



The remaining integrals are evaluated numerically with the Numerical Recipes [p2 
routines. 
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A. 2 Calculations for the mean coherence length 

In this appendix, some details of the calculation of the mean coherence length for the 
ggG-component, section |3.4| , are presented. We consider the case Q^ ^ b^. With the 



approximations (|3(J5D to (|3(J7D we obtain for the denominator of 



D^ = d^cd^rgg dag / -^ |^^ (r^^-a^^-)!^ [a^g- (s, rc)]^ r^g- ■ W^g (^g)' 

(363) 
This integral diverges logarithmically for a^™ ^ 0. To find an appropriate lower cut 
off, note that the mass of the qqG system is approximately given by 

<a - ^ + Q^. (364) 

We demand that M^^^ < 0.2s which leads to a^™ = 26q/(0.2s — Q'^). Furthermore we 

work in the approximation of ac <C 1 and we also have to choose an upper cut off. We 

use 

2bl _ 2b? 



The two limits become equal at xbj ~ 0.1. In our further calculation of D'^ we do the 
replacement rqqiVqqj — > r^^dij and perform the derivative. This yields 

2 G 



a'/ 



X K, is, tg)]' e-'o^o f\ + ^ + b{]. (366) 

For the integration over r^g we use the integral representation 

K„w4^°°fexp(-«-^) (367) 

of the modified Bessel function and obtain 



OO Q 

drq,r'^-Kl{erq,) = —. (368) 



Thus we have for the denominator 

Now we restrict ourselves to the a dipole cross section of the form 

aqq{s,rG)=C{s)rl (370) 
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and perform the last integration with the result 



^a ^ 32aemasC\s) ^^ a^ 



n^QX 



a 



(371) 



G 



Note that the factor C{s) will drop out, when one takes the ratio < P^ >= N'^ /D^. 
Next we calculate the numerator 



iV^ = rriNXRi i -J^j — ^4Q^ j (friG(fr2Gd'^raa I da, 






qq I ^^q 




a 1 - aq) 

ac ^ 



-..mm 
*G 



X Kl{er„.)aqAs,riG)(T„Js,r2G 



^ \7 ^' ^^^ -^^ 

^qq ■ Vr-iG 2 ^ ^ 

'IG 



'2G 



dAzGcG (^2G, ^G, A^;) 



(372) 



where 



Ggg(^2g, riG, A^) 



&^e- 



27rsinh(cjA^)^^^| 2 



h^ 
% 



{rlc + ^Ig) cth(u;A; 



2riG ■ fsG 



sinh(u;Az) 



is the solution of (|306|) with u = 6q/(z/q;g)- 



(373) 



Again we can do the replacement rggiVggj — > r^^Sij, perform the derivatives, sum 
over gluon polarizations and use ( p68| ) to obtain 



A^' 



G 



rriNXEj ■ 






da, 



hle- 



d^riGd'^r2GdAz / agq (s, tig) agg (s, r2G) ^ 



2v 



OtG 



sinh (cuA^l 






4 (ru^^Gf ^ ^4(nG-r;2G)^ ^ 26^ (^"^^ ■ ""^^^ 



4 4 
' 1G'2G 

"0 "0 



^2 ^2 
1 



'^1G'^2G 



^0 2 „4 

' IG' 2G 



26; 



2 I^IG • r2Gj 
^4 „2 
' IG' 2G 



2 2 2 2 

'IG 'IG '1G'2G 



exp (-/3 (r^G- + Ag) + 27^ig ■ ^2g) 



(374) 



where 



ft2 
/? = ^(l + coth(cuA^)) 



7 



ft2 



2 sinh ('(x'A^) 



(375) 
(376) 



With a cross section like agg (s, r) = C{s)r'^ the integrations over tg are easily per- 
formed with the result 



A^' 



G 



rriNXEj - 



37T^Q% 



dAz 



da 



_QfA2 

G e 2. 



10 



og sinh(t^A2) i (1 + coth (cuA^))' 
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12 



8 sinh (uAz) In 



1 + coth(a;Az)Y, 



■377) 



sinh^ (uAz) (1 + cth {ujAz)f 

^max 

= Sn^Q^l I ^^y ^r' (l-y)ln(l-,) 

In the last step we have introduced the new variables y = e~'^^^ and 5 = 2h'^/{Q'^aG)- 
According to ( p65| ) the limits for 5 are 



1<*<0,2^-1, 



(379) 



For the integral containing the logarithm it is convenient to do one more substitution, 
X = y^ . Then one finds 



/ (iy2/^2 (1 _ y^~) In (1 - t/2) 

Now only one integration is left in A^*^ 

5 



(380) 



- lim -— / dxx'^^ 2 (1 — x) "*"'' 



Q 64:aemasC^{s) f dS 



37T^Q%1 



6 



7 ^ / , .ON , M 1 



(382) 



and we end up with the result (|3TID for the factor (P^) = N'^/D'^ at Q^ > 6g, 



Smax 

"g 



(P-) 



31n'^ 



d6 
T 



7 



6 



(1 + 5) ' 8(1 + 35) 52-1 V"^^^^ ^[2^26 



(383) 



A. 3 Calculations for DY cross sections in pA scattering 

The calculation of the transverse momentum distribution of DY pairs in pA collisions 
requieres rather involved analytical and numerical calculations. In order to be able to 
perform some of the integrations in ( p37|) analytically we employ the approximations of 
a uniform nuclear density pA and a dipole cross section agg{ap) = Ca^p^ . With these 
approximations, it is possible to write down an analytical expression for the Green 
function, 



G{p2,Z2\pl,Zi) 



27r sinh (coAz) 



X exp<;-- 



(Pi + P2) "^o^h (ujA: 



2pi • P2 
sinh (uAz) 



. (384) 
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Here, 



a={-l + i)^PAEga^l-a)C/2 (385) 

and 

LJ = -{1+ i)^pACa/{2Eg (1 - a)) (386) 

is the complex oscillator frequency. 

We decompose (|337D in six integrals which correspond to the integration regions 



1. Xi : z,<Z2<-L = -^Ra - 62, 

2. J2 : zi < -L < Z2 < L, 

3. X3 : 2^1 < — L and L < Z2, 

4. J4 : -L< Zi<Z2< L, 

5. X5 : —L < Zi < L < Z2, 
Q. Tq : L < Zi < Z2. 

The differential cross section is then given by 



d{\na)d'^PT ._ 



For convenience we introduce the notation 

\2 



J2^r (387) 



Hi = l + {l-a)\ (388) 

Ho = 2M^ (1 - af (389) 

and employ the auxiliary functions 

2 

Giz,t) = I (390) 

?7^ cosh^cjzj + 2at smh^cjzj 

= 1 2atcosh(^2)+r,2sinh(a..) 

2a ?7^ cosh(cijz) + 2at smh(ci;z) 

F{L,z) = 2Ca^pA{2L-z). (392) 

The calculations are too long, to be shown in detail. The result is 

I, = ^ iHo + H,pl) 

Ra 2L 

I2 = -^-^2Re^^^^^J LdL I dz[Ho + H,p^^G{z,0)]G{z,0) 



X exp{-iq'^'''z-p'^Tp{z,0)), (394) 
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X. 



2a. 



-2Re 



n 



^2 (^^2 +p2^) 



f LdL f dt [Ho + HipI G{2L, t)] G{2L, t) 



exp (-igr'^2L - t - p^/?(2L, t)) 



X4 



2ap 



-2Re 



(C 



TT 



LdL f dz f dz'G{z',0) 



Ho + H,G{z',0) 





FiL,z) 



Pt 



1 + F(L, 2;)/3(2;', 0) (1 + F(L, 2)/3(^', 0))^ 



X 



«^P \^ i^L ^ l+F{L,z)l3{z',0) 

l + F{L,z)P{z',0) ' 

Ra 2L oo 



I, 



=^^2Re^ 



TT 



Tj" 



LdL I dz I dtG{z,t) 



Ho + H^G{z,t) 





F{L,z) 



Pt 



1 + F(L, 2;)/3(^, t) (1 + F(L, ^)/?(;2, t))^ 



X 



«^P ^ iVL -^ '^ l+F(L,2)/3(2,t) 



l+F(L,2)/?(^,t) 



X« 



2a„ 



iJA 



TT 



L(iL / dt 



E^t + H,il-t) 



exp I —t 



Pt 



ri^+F{Lfi)t 



r/2 + F(L,0)t 



(395) 



(396) 



(397) 
(398) 



Note that in Xi and Xg, the integration over the longitudinal coordinates Zi and Z2 
does not extend over the nucleus, which is located between —L and L. Nevertheless, Xi 
and Xg are different from zero. One can however check, that the cross section vanishes 
in the limit p^ — »■ 0. In particular one finds Xi + X3 + Xg = — X2 = X4 = — X5. 



A. 4 Calculations for shadowing in DIS 

In the approximation where only the gg-Fock-component of the virtual photon is taken 
into account, it was found that nuclear shadowing in DIS is given by the formula 



(T. 



7*A 
tot 



Nf 

A(T^*P--ReV 
2 ^^ 



d% 



/•oo /»1 

dzi / dz2 da I d^pi / d'^p2 

Jzi Jo 



X \1/ 



T,L 



', A, P2) Pa (^ 2^2)0^59 (S,P2) 



X G (p2, 2:2 1 Pi, zi) Pa {b, zi) aqq (s, pi) '^t,l {e, A, pi) 



(399) 



where summation over the polarisations (T, L) of the virtual photon is understood. The 
Greenfunction G describes the propagation of the gg-pair through the nucleus. With 
the nonperturbative interaction which was introduced to model confinement effects, G 



A APPENDIX 125 



fuUfills the equation 



d Ai (02) — e^ % ,, . , , a^ [a] pi 
0Z2 2ua[l — a) 2 2ua[l — a) 



G{p2,Z2\pl,Zi) 



i5{z2-zi)5^^\p2-pi), (400) 



where 



e^ = a (1 - a) Q2 + mj (401) 

is the so-called extension parameter which depends on the flavor / via the quark mass. 
The strength of the potential between the particles is characterized by the parameter 

a^ (a) = v^-^\ll2MeVy + (1 - vf'^^ (leSMel/)^^ (1 - a) . (402) 

Since ( |399|) is too complicated for a direct numerical evaluation, we introduce the 
following simplifying approximations. First, we chose a dipole cross-section of the 
functional shape 

a,,-(s,p) = C(s)p2, (403) 

where the coefficient C{s) does not depend on p. Second, we work with a uniform 
nuclear density 

PA {b, z) = pAoQ {R\ -b^- z^) . (404) 

These two approximations allow us to solve ( |400| ) analytically, because ( [40 0| ) reduces 
to the Schrodinger-equation of a two dimensional harmonic oscillator, where z plays 
the role of time. The solution is known as 

G{p2,Z2\pi,zi) -- 



27r sinh {ljAz] 
B(a) 



X exp 



2 



{pl + pl)coth{u;Az) '^^''^' 



sinh (uAz) 



.(405) 



Note however that the frequency of the oscillator has a complex value. 



u = 1^^. (406) 

Furthermore we introduced the notation 

Az = Z2- zi (407) 

and 

^2 {a) = -ipAMC{s) + a^ (a) . (408) 

The quantity 

M = ua{l-a) (409) 

plays the role of the reduced mass of the pair and the minimal longitudinal momentum 
transfer when the photon sphts is given by 

Ql = • (410) 

^^ 2M ^ ^ 
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With these simphfications ( |399D reads 

1 ^ r n/^ 



/=i 



dzi 

1/2 Jo 



l/2~zi 



d{Az) I da I d^pi 



d^p2 



X e-'^-(^-^^)p^oC'(^)P?P2 



$0 (^, A, pi) $0 (e, A, P2) + <^i (£, A, pi) ■ $1 (e, A, P2) 



(411) 



5 



27rsinh('u;A^)^^^ I 2 



B 



(Pi + P2) coth (cuA; 



2pi • P2 
sinh iuj/S.z) 



with / = 2a/-R^ — &^. For simphcity we do write out the dependence of -B on a 
exphcitly. From the physical point of view it is reasonable to distinguish between 
shadowing for transverse and longitudinal photons, as indicated by the subscripts T 
and L in ( |399| ). For calculational purposes it is however more convenient to use the 
functions 



$o(e,A,pi)<l>o(e,A,p2) 

00 

A2 
dudt 



QaemZf 



(mj + AQ'a' (1 - af) 



(412) 







sinh (At) sinh (An) 



exp 



;^g-2 2 Xe'^o'^ 

^^' coth (At) - t -p- coth (An) - n 



4 



$1 (e, A, pi) ■ $1 (5, A, p2) = ^^^ (1 - 2« (1 - «)) ^ 

(27rj P1P2 



(413) 



X / dudt exp 





— -p- coth (At) - t — coth (An) 



4 



n 



instead. The strength of the nonperturbative interaction is now characterized by the 
dimensionless parameter 

A = ^. (414) 

Note that in the limit A ^ one obtains the expressions for the perturbative wave- 
functions. We will not write out the a dependence of A explicitly in the following. 
In particular, $0 goes to the Ko-pait of the light-cone wave- functions and $1 to the 
iiTi-part. 

The total cross section '~f*A can now be decomposed into three terms. 



a 



r^ = Aa^'P - a{^ - aj 



7* A 



(415) 



The shadowing correction is given by the last two terms, where the indices and 1 
refer to the $0 and $i-part of the wavefunctions, respectivly. 

First we consider the $o-part. Since the integrand in ([41 1| ) depends on zi only via 
the upper limit of the integration over Z2, we are able to integrate two times by part. 
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with the result 



cr, 



■y*A 



-| f pD fl /• /• 1*00 

- Re y^ / dl da d^pi I d^p2 / dudt 



TX 



12 



if - ?,DH + 2D') e-'^^'C\s)pl,^^^ {m} + 4Q'a' (1 - a)') 



2„--u-t 



A^e 



B 



{Any 
-pIpI exp (-7p2 _ /3p2 ^ 25pi ■ pi) . (416) 



sinh (At) sinh (Am) 27rsmh [uol] 

The polynomial in / results from making the boundary terms due to the integration by 
parts vanish. We also introduced the nuclear diameter D = 2Ra and the notation 

B coth {uol) Ae^ coth (At) 



7 
5 



2 4 

B coth {ul) \e^ coth (Au) 



2 



4 



(417) 
(418) 
(419) 



2 sinh (a;/) 

The exponential in ( 416 ) is obtained when one inserts the expressions for the Green 
function and the nonperturbative wave-functions into ([41 1|) . Now one notes that the 
integrations over pi and p2 are simply Gaussian and can immediately be performed, 



Nf 



(T, 



■y'A 



Oif., 
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C\s)p\oReJ2z] / dl da 



/=i 



D 



dudt 



e-'"^' (m^f + AQ^a^ (1 - af) - 



^0 



B 



sinh (ujI) 



sinh (At) sinh (Am) 
252 



(420) 



Thus we have reduced the number of integrations to four. Substituting ([417|) in ( |420|) 
one obtains the somewhat lengthy expression 



Nf 



af^ = 32aemC^{s)plf^R\Re^ZJ dL da 



X e 






dudt 



(421) 



X < 



+ 



{m} + 4QW (1 - af) {L' - 3L' + 2) e'^-^A^i? 

sinh {uoDL) sinh (At) sinh (Am) 
[sinh {uDL) {AB^ sinh (At) sinh (An) + A^e^ cosh (At) cosh (An)) 

+25 Ae^ cosh (cuDL) sinh (A(m + t))f 

8B^ sinh^ (At) sinh^ (Am) 

[sinh (uDL) (45^ sinh (At) sinh (An) + X^e^ cosh (At) cosh (An)) 

+2B\e^ cosh (uDL) sinh (A(n + t))f 
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The $i-part can be treated in the same way. Once again we integrate two times by 
part and obtain 



^f rD rl 






- Re y^ / dl da d^pi / d'^p2 / dudt 



X 



2 ^<^emZ f 



^ (/^ - ^DH + 2D') e--^^C\s)p\,^^ (1 - 2a (1 - a)) 



X e"""*^— ^^Pi ■ P2 exp (-7P? - PpI + 25pi ■ P2) (422) 

with the same notations as in ( |417| ). The Gaussian integral over pi and p2 is performed 
with the result 



^2/ \ 2 '^-f pD pi 



32 



Re 2^ ^f / dl I da I dudt 
r_i Jo Jo Jo 



X 



(/^ - 3D^/ + 2D'') e-'""^^ (1 - 2a (1 - a)) 



B 5 



X e-"-* ■ w n 2' (423) 

smh (cu/) (7/? - ^2)2 

and after substitution of the expressions in ( |417|) one finally obtains 

TV, 



7*A 

a.' 



J pi pi poo 

QAaemC'^is)p\oR\Re^ZJ / dL da dudt (424) 

r ^ Jo Jo Jo 



f = 

u~t d2 



^-ig^DL ^^ _2a{l- a)) {L' - 3L^ + 2) e-"-*5 

sinh^ (At) sinh^ (Am 



[sinh {cuDL) (45^ sinh (At) sinh (Am) + A^e^ cosh (At) cosh (An)) 

+2BXe^ cosh (ooDL) sinh (A(n + t))]^ 

Taking the limit A — *> 0, one obtains the corresponding expressions for the pertur- 
bative wave-functions, 

f pi pi POO 

a^'*^ = 32aemC'^{s)p\QR\Re^Z] / dL da dudt (425) 

r^ Jo Jo Jo 

X e-'"^^^ (mj + 4Q2«2 (i _ ^f) ^l' - ?>L^ + 2) e'^-'A^E 

sinh (uDL) ut 



X 



[sinh (ujDL) {AB^ut + e^) + 25^2 cosh (luDL) {u + t)f 



8B^uH^ 
+ 



[sinh (cjL)L) {AB^ut + e^) + 25e2 cosh (ooDL) {u + t)] 

^f /•! .1 

/=1 



3 I ' 



f p1 p1 poo 

a7*^ = 6AaemC'^{s)plQR\Rey2zj / dL da dudt (426) 

r_i io io io 
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^-ig^DL (^ _ 2a (1 _ a)) {L' - 3L^ + 2) e' 



u-t Ty2 



uH^ 



[sinh {ujDL) {AB^ut + e^) + 2Be^ cosh (uDL) {u + t)Y 



To make our simplifications ( [40 3D and ( [404| ) more realistic, we determine C(s) by 
demanding, 



2jd%d^r\<I/T,L{er) 



exp 



CT.Lr-' 



T(6) 



Jd\\^T,Lier)\'CT,Lr' 



2jd%d^r\^T,L{er) 



1 — exp 



o-<;g(s/'-) 



^T(6) 



/rfV|^T,L(£r)| Crgg-(s, 



where 



Tib) 



dz Pa {b, z) 



(427) 



(428) 



is the nuclear thickness. This means that C{s) is chosen such that in the limit of 
infinite coherence length, the same value of the cross section as obtained with the 
realistic parametrization ( |152| ) is reproduced. This procedure is performed separately 
for transverse and longitudinal photons and for each value of a. We use realistic 
parametrizations for pa(^; z) with parameter values from | ]189 |. 

The value for the mean nuclear density p^o is determined in a similar way from 



d:'b 



l-exp(-^2^/i?i-6V 



AO 



d% 



l-exp(-|T(6) 



(429) 



We tried values between 1 mb and 50 mb for ctq and found the result for p^o practically 
independent of this value. 
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